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1. Synopsis 



This report summarizes the activities of the working group on the "Minimal Supersym- 
metric Standard Model" or MSSM for the GDR-Supersymetrie. It is divided into five 
main parts: a first part dealing with the general features of the MSSM spectrum, a sec- 
ond part with the SUSY and Higgs particle physical properties, and then two sections 
dealing with the production and decay properties of the MSSM Higgs bosons and of the 
supersymmetric particles, and a last part summarizing the experimental limits on these 
states. A brief discussion of the prospects of our working group are given at the end. In 
each section, we first briefly give some known material for completeness and to establish 
the used conventions/notations, and then summarize the original contributions to the 
GDR-Supersymetrie workshop[[ For more details, we refer to the original GDR notes 
[1-18]. The report contains the following material: 

Section 2 gives a survey of the MSSM spectrum. We start the discussion by defining 
what we mean by the MSSM in section 2.1, thus setting the framework in which our 
working group is acting [||. We first introduce the unconstrained MSSM (uMSSM) defined 
by the requirements of minimal gauge group and particle content, R-parity conservation 
and all soft-SUSY breaking terms allowed by the symmetry. We then briefly summarize 
the various phenomenological problems of the uMSSM which bring us to introduce what 
we call the "phenomenological MSSM" (pMSSM) which has only a reasonable set of free 
parameters. We then discuss the celebrated minimal Supergravity model (mSUGRA), in 
which unification conditions at the GUT scale are imposed on the various parameters, 
leading only to four continuous and one discrete free parameters. We finally discuss the 
possibility of relaxing some of the assumptions of the mSUGRA model, and define several 
intermediate MSSM scenarii (MSSMi), between mSUGRA and the pMSSM. 

In section 2.2, we focus on the [radiative] Electroweak Symmetry Breaking mechanism, 
starting by a general discussion of the two requirements which make it taking place in 
the MSSM. We then summarize some analytical results for theoretical bounds on the 
parameter tan (3, obtained from the study of electroweak symmetry breaking conditions 
at the one-loop order ]2], [|. The key-point is to use, on top of the usual conditions, the 
positivity of the Higgs boson squared masses which are needed to ensure a [at least local] 
minimum in the Higgs sector. We then show analytically how these constraints translate 
into bounds on tan f3 which are then necessary and fully model-independent bounds. The 
generic form of the bounds are given both in the Supertrace approximation and in a more 
realistic one, namely the top/stop-bottom/sbottom approximation. 

In section 2.3, we first write down for completeness the Renormalization Group Equa- 
tions for the MSSM parameters in the one-loop approximation, that we will need for later 
investigations. In a next step, we discuss some exact analytical solutions of the one-loop 
renormalization group evolution equations of the Yukawa couplings |3J]. Solutions of the 
one-loop equations for the top and bottom Yukawa couplings exist in the literature only 

1 Since the spectrum of the various contributions is rather wide, this leads to the almost unavoidable 
situation of "un rapport quelque peu decousu"... 
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for some limiting cases. Here, we give the exact solutions with no approximation: they are 
valid for any value of top and bottom Yukawa couplings and immediately generalizable 
to the full set of Yukawa couplings. These solutions allow to treat the problematic large 
tan /3 region exactly, and thus are of some relevance in the case of bottom-tau Yukawa 
coupling unification scenarii. 

Section 3 is devoted to the physical parameters of the MSSM. In section 3.1. we 
first briefly summarize the general features of the chargino/neutralino, sfermion mass 
matrices to set the conventions and the notations for what will follow, and discuss our 
parameterization of the Higgs sector exhibiting the Higgs boson masses and couplings that 
we will need later. We then discuss a method which allows to derive the MSSM Lagrangian 
parameters of the gaugino sector from physical parameters such as the chargino and 
neutralino masses ||. This "inversion" of the MSSM spectrum is non-trivial, especially 
in the neutralino sector which involves a 4 x 4 matrix to de-diagonalize. The algorithm 
gives for a given tan/3, the values of the supersymmetric Higgs-higgsino parameter /i 
and the bino and wino soft-SUSY breaking mass parameters Mi and M2, in terms of 
three arbitrary input masses, namely either two chargino and one neutralino masses, or 
two neutralino and one chargino masses. Some subtleties like the occurrence of discrete 
ambiguities are illustrated and a few remarks on inversion in the sfermion and Higgs 
sectors are made. 

Section 3.2 describes the Fortran code SUSPECT which calculates the masses and 
couplings of the SUSY particles and the MSSM Higgs bosons ||. The specific aim of 
this program is to fix a GDR-common set of parameter definitions and conventions, and 
to have as much as possible flexibility on the input/output parameter choice; it is hoped 
that it may be readily usable even with not much prior knowledge on the MSSM. In the 
present version, SUSPECT1.1, only the two extreme "models" pMSSM on the one side 
and mSUGRA on the other side are available. We describe the most important subrou- 
tines in the case of the pMSSM and then discuss the mSUGRA case with the various 
choices of approximations and refinements, paying a special attention on the renormaliza- 
tion group evolution and the radiative electroweak symmetry breaking. We then briefly 
describe the input and output files and finally collect a list of various available options 
and/or model assumptions, with a clear mention of eventual limitations of the present ver- 
sion of the code, and list the improvements which are planed to be made in the near future. 

In section 4, we discuss some aspects of the production of the MSSM Higgs bosons at 
future hadron and e + e~ colliders, and their possible standard and SUSY decays. We start 
in section 4.1 with Higgs boson production at the LHC and first make a brief summary 
in section 4.1.1 of the main production mechanisms of the lightest Higgs boson h [in 
particular in the decoupling limit where it is SM-like] and the heavier particles H, A and 
H , as well as the most interesting detection signals. 

In section 4.1.2 we analyze the effects of stop loops on the main production mechanism 
of the lightest h boson at the LHC, the gluon fusion mechanism gg — > h, and on the 
important decay channel h —> 77 |7|]. We show that if the off-diagonal entry in the t 
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mass matrix is large, the lightest stop can be rather light and its couplings to the h boson 
strongly enhanced; its contributions would then interfere destructively with the ones of 
the top quark, leading to a cross section times branching ratio cr(gg — > h) x BR(/i — > 77) 
much smaller than in the SM, even in the decoupling regime. This would make the search 
for the h boson at the LHC much more difficult than expected. Far from the decoupling 
limit, the cross section times branching ratio is further reduced due to the additional 
suppression of the Higgs couplings to SM fermions and gauge bosons. In the case of the 
heavy H boson, squark loop contributions to the cross section <r{gg —>■ H) can be also 
large, while they are absent for the A boson because of CP-invariance. 

In section 4.1.3, we discuss the production of the light h particle in association with 
light top-squark pairs at proton colliders, pp — > gg + qq — ► ith ||. The cross section for 
this process can substantially exceed the rate for the SM-like associated production with 
top quarks, especially for large values of the off-diagonal entry of the t mass matrix which, 
as mentioned previously, make the lightest stop much lighter than the other squarks and 
increase its coupling to the h boson. This process can strongly enhance the potential 
of the LHC to discover the h boson in the 77+ lepton channel. It would also allow for 
the possibility of the direct determination of the tth coupling, the largest electroweak 
coupling in the MSSM, opening thus a window to probe directly the trilinear part of the 
soft-SUSY breaking scalar potential. Finally, this reaction could be a new channel to 
search for relatively light top squarks at hadron colliders. 

In section 4.2 we analyze the possible decays into SUSY particles of the neutral 
and charged [10 1 Higgs particles of the MSSM. For the light h boson the only SUSY de- 
cay allowed by present experimental data are the invisible decay into a pair of lightest 
neutralinos or sneutrinos. The decays are possible only in small areas of the parameter 
space in the constrained MSSM; however, relaxing for instance the assumption of a uni- 
versal gaugino mass at the grand unification scale, leads to possibly very light neutralinos 
and the decays into the latter states occurs in a much larger area. Decays of the heavy 
neutral H, A bosons into chargino/neutralino pairs and H boson decays into stop pairs 
can be also dominant in some areas of the parameter space. We then show that the de- 
cays of the particles into chargino/neutralino and slepton pairs are also still allowed 
and can be dominant in some areas of the parameter space; we also briefly discuss some 
additional charged Higgs boson decay modes present in non-supersymmetric two-Higgs 
doublet models. The SUSY decays should not be overlooked as they can strongly suppress 
the branching ratios of the Higgs boson detection modes, and therefore might jeopardize 
the search for these particles at the LHC. 

In section 4.3 we briefly summarize the main production mechanisms of the MSSM 
Higgs bosons at e + e~ colliders and describe the Fortran program HPROD [IT] which cal- 
culates the production cross sections for SM and MSSM Higgs particles in e + e~ collisions. 
In the SM, it includes the bremsstrahlung off the Z-boson line and the WW/ZZ fusion 
processes; some higher order production processes, such the production in association with 
tt pairs and the Higgs boson pair production in the bremsstrahlung and the WW/ZZ fu- 
sion processes, are also included. For the MSSM CP-even Higgs bosons, it includes the 
Higgs-strahlung, the associated production with the pseudoscalar Higgs boson A, and 
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the WW/ZZ fusion processes; for the H ± boson it includes the pair production in e + e~ 
collisions as well as the top quark decay process. The complete radiative corrections in 
the renormalization group improved effective potential approach are incorporated in the 
program, which computes both the running and pole Higgs boson masses. The possibili- 
ties of having off-shell Z or Higgs boson production in the bremsstrahlung and in the pair 
production processes, as well as initial state radiation, are allowed. Future improvements 
will be listed. 

Finally, in section 4.3.2, we discuss the associated production of the lightest h boson 



with stop pairs in e + e~ collisions [12]. The final state t\tih can be generated in three 
ways: (i) the production of e + e~ — > t 2 t\ through Z-boson exchange and the subsequent 
decay £2 — ► t\h; (ii) the production in the continuum in e + e~ collisions e + e~ — > t\tih with 
the main contribution coming from h emission from the t\ lines; and (Hi) the production 
in the 77 mode of the e + e _ collider, 77 — ► t\t-\h. Due to the clean environment of e + e~ 
colliders, this final state might be easier to be detected than at the LHC if kinematically 
allowed, and would provide a more precise determination of the t\t\h coupling. 

In Section 5, we discuss some of the production and decay properties of the SUSY 
particles as well as their virtual effects. We start in section 5.1 by discussing some of 
the virtual effects of charginos, neutralinos and sleptons of the first generation at LEP2 
energies [|H|. In the production of lepton pairs in e + e~ collisions, box diagrams involving 
neutralino/selectron or charginos/sneutrino pairs occur and alter the production cross 
sections and asymmetries; at LEP2 energies the effects can be sizeable and experimentally 
measurable in the process e + e~ — > fi + n~ if the masses of the neutralinos and/or charginos 
are close to the beam energy due to threshold effects. 

In section 5.2, we analyze the correlated production and decay of a pair of the lightest 
charginos |141P|. We show that the chargino polarization and the spin-spin correlations give 
two observables which do not depend on the final decays of the charginos, and hence on 
the neutralino sector. Combined with the production cross section and with the chargino 
mass which can be measured via a threshold scan, these two observables allow a complete 
determination of the SUSY parameters /1, M 2 and tan /3 [and the sneutrino mass] in a 
completely model-independent way. With the knowledge of the lightest neutralino mass 
from the energy distribution of the final particles in the chargino decays, the parameter 
Mi can be also determined, leading to a full reconstruction of the gaugino sector. 

Section 5.3 deals with the production of a chargino/neutralino pair at proton colliders, 
pp — > qq — ► Xi %2- We first discuss the analytical expression of the tree-level cross section 



15|], that is calculated with the help of the FeynMSSM package that we briefly describe, 
and compare with the results available in the literature. We then discuss this process 
at the LHC, and propose to use the charge asymmetry in the three-leptons plus missing 
energy channel to determine some MSSM parameters [flq| . In particular, it will be shown 



that this asymmetry depends only on the mass of the final state, and can be used to 



2 This contribution has not been, strictly speaking, entirely made in the framework of the GDR- 
Supesrymetrie. However, one of the authors could not resist to the temptation of including it in this 
report since it is one of the hot topics of the GDR, especially in the working group "Outils" . 
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measure the sum of the chargino and neutralino masses in a model-independent way. 



Some important three-body decay modes of sfermions |17| are discussed in section 
5.4. In particular, we discuss the decays of the lightest stop into a 6-quark, the lightest 
neutralino and a W or H + boson which can compete with the loop-mediated decay into 
charm+neutralino in the case of light stops. We also discuss decays of the heavier stop 
(sbottom) into the lighter one and a fermion pair through off-shell gauge or Higgs bosons, 
as well as decays of squarks into third generation sleptons+leptons (squarks+quarks) 
through a virtual exchange of a neutralino/chargino (gluino). These decays can have 
sizeable branching fractions in some areas of the MSSM parameter space. 

In section 5.5, we discuss stop and sbottom squark searches at LEP200 These 
squarks have a special place in the SUSY spectra due to the strong Yukawa couplings of 
their partners, the t and b quarks, and could have masses accessible at LEP200. After a 
brief analysis of the cross section including all important radiative corrections, we discuss 
the two main decay modes which are relevant for stop squarks with masses accessible 
at LEP200, namely the loop induced flavor changing decay into a charm quark and the 
lightest neutralino, t± —>■ cxI, and the three-body decay t\ — > bl + v through the exchange 
of an off-shell chargino. Some remarks will be given on the decays of a light sbottom. 



Finally, section 6 deals with the limits and constraints on the SUSY particle and MSSM 
Higgs boson masses from present experiments fl9j , mainly from LEP and the Tevatron. 
At LEP, the searches for SUSY particles concern sleptons, stops, sbottoms, charginos 
and neutralinos, while at the Tevatron the main focus is on squarks and gluinos. The 
lightest CP-even h and the CP-odd A bosons are also searched for at LEP, while the 
charged Higgs bosons H ± are searched for at the Tevatron. Here, we will summarize the 
experimental limits on the production cross sections and masses of these particles at LEP2 
and the Tevatron, paying attention to all the decay channels. The most recent results 
have been used, including preliminary results reported at the last summer conference in 
Vancouver and the last LEPC meeting. 
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2. The MSSM Spectrum 



2.1. Definitions and properties 

The Minimal Supersymmetric Standard Model (MSSM) is the most economical low- 
energy supersymmetric (SUSY) |H| extension of the Standard Model (SM); for reviews 



see Refs. pT| , [22] , p3|] . In this section, we discuss the basic assumptions which define the 
model and the various constraints which can be imposed on it. This will allow us to 
set our notations and conventions for the rest of the discussion. We will mainly focus 
on the unconstrained MSSM (uMSSM), what we will call the phenomenological MSSM 
(pMSSM) and the constrained minimal Supergravity (mSUGRA) model. 

2.1.1 The uMSSM: unconstrained MSSM 

The unconstrained MSSM is defined by the following four basic assumptions: 

(a) Minimal gauge group: 

The MSSM is based on the gauge symmetry SU(3) C X SU(2) L x U(1) Y , i.e. the SM 
symmetry. SUSY implies then, that the spin-1 gauge bosons, and their spin-1/2 super- 
partners the gauginos [bino B, winos Wis and gluinos Gi_s] are i n vector supermultiplets. 

(b ) Minimal particle content: 

In the MSSM, there are only three generations of spin-1/2 quarks and leptons [no 
right-handed neutrino] as in the SM. The left- and right-handed chiral fields belong to 
chiral superfields together with their spin-0 SUSY partners the squarks and sleptons: 

Q , u R , d R , L , l R (2.1) 

In addition, two chiral superfields Hd, H u with respective hypercharges —1 and +1 for 
the cancellation of chiral anomalies, are needed [[M, |25|. Their scalar components: 




give separately masses to the isospin +1/2 and —1/2 fermions. Their spin-1/2 superpart- 
ners, the higgsinos, will mix with the winos and the bino, to give the mass eigenstates, 
the charginos Xi,2 an d neutralinos Xi2 3 4- 
(c) R-parity conservation: 

To enforce lepton and baryon number conservation, a discrete and multiplicative sym- 
metry called R-parity is imposed PEfl . It is defined by: 

R = (_1)2*+3B+L (2.3) 

where L and B are the lepton and baryon numbers, and s is the spin quantum number. 
The R-parity quantum numbers are then R = +1 for the ordinary particles [fermions, 
gauge and Higgs bosons], and R = —1 for their supersymmetric partners. In practice the 
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conservation of i?-parity has important consequences: the SUSY particles are always pro- 
duced in pairs, in their decay products there is always an odd number of SUSY particles, 
and the lightest SUSY particle (LSP) is absolutely stable. 

The three conditions listed above are sufficient to completely determine a globally su- 
persymmetric Lagrangian. The kinetic part of the Lagrangian is obtained by generalizing 
the notion of covariant derivative to the SUSY case. The most general superpotential [i.e. 
a globally supersymmetric potential] compatible with gauge invariance, renormalizability 
and R-parity conserving is written as: 



W 



E -Y% u R H u .& + Y* d R H d .Q j + l R H u .U + ^H u .H d 



(2.4) 



t,]=gen 



The product between SU(2) L doublets reads H.Q = e a bH a Q b where a, b are SU(2) L indices 
and 6i2 = 1 = —£21, and Y^' d ' 1 denotes the Yukawa couplings among generations. The 
first three terms in the previous expression are nothing else than a generalization of the 
Yukawa interaction in the SM, while the last term is a globally supersymmetric Higgs 
mass term. The supersymmetric part of the tree-level potential Vt ree is the sum of the 
so-called F- and D-terms [27], where the F-terms come from the superpotential through 
derivatives with respect to all scalar fields (fi a 



v F = El^l 



W a = dW/d</> a 



(2.5) 



and the D-terms corresponding to respectively U(1)y, SU(2)l, and SU(3)c gauge sym- 
metries are given by 

(2.6) 



V D = -{D X D X + D 2 D 2 + D 3 D 3 



with 



9i 



92 



9s 



i=gen 



1 T 



-L\U - -u Ri u Ri + -d R 4 Rz + VjJr.) + -HtH u - -H d H d 



i=gen " — — 

~ A ~ A* ~ A* ~ 

E Q\lfii - u Ri —u Rt - d ] R —d Ri 

i=gen " ~ 



(2.7) 



Here the tildes denote the scalar quark and lepton fields, and (cxfc)fc = i_3 and (Xk)k=i-8 the 
Pauli and Gell-Mann matrices; #1,2,3 are the three gauge couplings. 

For completeness we also write down the fermion-scalar sector of the supersymmetric 
part of the Lagrangian, which is needed for later purpose. Following the conventions of 
Ref. [2(J , with two-component spinors, this part of the Lagrangian contains on the one 
hand a purely chiral contribution 



£ch ir . = --Er k tt + h.c. 

a, b 



W ab = d A W/d<p a d<p b 



(2.8) 
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where ijj a is the supersymmetric fermionic partner of (j) a , and on the other hand a mixed 
chiral-vector contribution coming from the gauged matter field kinetic term yielding in 
component form 

£ mix = -iV2£,(D A ) a if> a \ A + h.c. , (D A ) a = dD A /d<p a (2.9) 

A, a 

The index A denotes the gauge group with which the gaugino X A is associated [there is 
of course also a contribution from the usual gauged kinetic term for the Higgs which we 
do not write here]. 

(d) Soft-S US Y breaking: 

To break Supersymmetry, while preventing the reappearance of the quadratic diver- 
gences [soft-breaking] we add to the supersymmetric Lagrangian a set of terms which 
explicitly but softly break SUSY |28| : 



Mass terms for the gluinos, winos and binos: 
1 



- C - - 

^gaugino 2 



3 8 

M X BB + M 2 W a W a + M 3 J2 G a G a + h.c. 

a=l a=l 



(2.10) 



Mass terms for the scalar fermions: 

,2 r\\r\ i _2 ft 



- ^sfcrmions = m Q,iQiQi + m L,i L i L i + m l,i\ U ^\ 2 + m S,i\ d Ri\ 2 + m l,i\ l ^\ 2 

i=gen 

(2.11) 

• Mass and bilinear terms for the Higgs bosons: 

- £ffi ggs = rn 2 Hu HlH u + m 2 Hd H j d H d + B/i(H u .H d + h.c.) (2.12) 

• Trilinear couplings between sfermions and Higgs bosons 

i,j=gen 

(2.13) 

The soft-SUSY breaking scalar potential, which will be discussed later in detail, is the 
sum of the three last terms: 

^soft -^sfermions -^Higgs ^tril. (2-14) 

Up to now, no constraint is applied to this Lagrangian, although for generic values of 
the parameters, it might lead to severe phenomenological problems, such as flavor chang- 
ing neutral currents [FCNC] , unacceptable amount of additional CP-violation, color and 
charge breaking minima, etc... The MSSM defined by the four hypotheses (a)-(d) above, 
will be called the unconstrained MSSM or in short the uMSSM. 
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2.1.2 The pMSSM: phenomenological MSSM 



The uMSSM contains a huge number of free parameters, which are mainly coming from 
the scalar potential V so ft. Indeed, the sfermions masses are in principle 3x3 hermitian 
matrices in generation space, with complex matrix elements, leading to 5 x 6 x 2 arbitrary 
parameters. The matrices A u and A d for the trilinear couplings on the other hand are 
arbitrary 3x3 complex matrices in generation space leading to 2 x 9 x 2 = 36 parameters 
[the A 1 matrices are diagonal and the matrix elements are real, since in the SM the neu- 
trinos are masslessQ and there is a separate conservation of the e, \x and r lepton numbers, 
leading to a much smaller number of parameters] . Thus, if we allow for intergenerational 
mixing and complex phases, the soft-SUSY breaking terms will introduce a huge number 
of unknown parameters, 105 parameters |$0| in addition to the 19 parameters of the SM! 
This feature of course will make any phenomenological analysis a daunting task, if possi- 
ble at all, and in addition induces severe phenomenological problems as mentioned above. 
One definitely needs to reduce the number of free parameters to be able to use the model 
in a reasonable and somewhat predictive way. 

There are, fortunately, several phenomenological constraints which make some assump- 
tions reasonably justified to constrain the uMSSM. These assumptions will be discussed 
in the report of the "Saveurs" group [31] to which we refer for details. Here we will simply 
and briefly mention them. 



(a) No new source of CP-violation 

New sources of CP-violations are constrained by the experimental limits on the elec- 
tron and neutron electric moments and in the K system [e.g. e'/e] which are extremely 
tight. [Of course, since the phases at hand in the uMSSM are numerous, a kind of fine 
tuning can be made which will allow for cancelling contributions in the various quanti- 
ties]. Assuming that all phases in the soft-SUSY breaking potential are zero eliminates all 
new sources of CP-violation and leads to a drastic reduction of the number of parameters. 

(b) No Flavor Changing neutral currents 

The non-diagonal terms in the sfermion mass matrices and in the trilinear coupling 
matrices, can induce large violations of FCNC which are severely constrained by present 
experimental data. Constraints have then to be imposed to suppress operators which 
lead to these large effects. These constraints amount to a severe limitation of the pattern 
of the sfermion mass matrices: either they are close to the unit matrix in flavor space, 
or they are almost proportional to the corresponding fermion masses [flavor universality 
and flavor alignment respectively]. We will assume here that both the matrices for the 
sfermion masses and for the trilinear couplings are diagonal, which also leads to a drastic 
reduction of the number of new parameters. 



3 A signal for neutrino oscillations, thus implying neutrino masses, has been very recently reported 
by Super-Kamiokande |p9| . However, these neutrino masses are so tiny that they have no effect on the 
discussion given here. 
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(c) First and Second Generation Universality 

Experimental data, e.g. from K°-K° mixing, severely limit the splitting between the 
masses of the first and second generation squarks, unless squarks are significantly heavier 
than 1 TeV. One can assume therefore that the soft-SUSY breaking scalar masses are the 
same for the first and second generations. There is no experimental constraint on the third 
generation masses [note in addition that in this sector significant mass splitting between 
the mass eigenstates can be generated by the off-diagonal matrix elements, as will be 
discussed later]. Furthermore, one can assume also that A u ,A d and A 1 are the same for 
the two generations. In fact, since they are always proportional to the fermion masses, 
these trilinear couplings are important only in the case of the third generation; one can 
therefore set the ones of the first two generations to zero without any phenomenological 
consequence in this context. 

In addition, some parameters in the Higgs sector can be related to SM parameters, 
see below. Thus, making the three assumptions (a)-(c) will lead to the following 19 input 
parameters only: 

tan/5: the ratio of the vev of the two-Higgs doublet fields. 

Ma- the mass of the pseudoscalar Higgs boson 

//: the Higgs-higgsino mass parameter 

Mi, M2, M3: the bino, wino and gluino mass parameters. 

mq,mu R ,m^ R ,rni,me R : first /second generation sfermion masses 

mQ,mt R ,m~ bR ,mi,m fR : third generation sfermion masses 

A t , Ab, A T : third generation trilinear couplings. 

Note that the remaining three parameters m 2 H , m 2 H and B are determined through 
the electroweak symmetry breaking conditions and the value of M^', alternatively, one can 
use directly the Higgs mass relations, which are equivalent to the electroweak symmetry 
breaking conditions, only when supplemented with an extra relation [see section 2.2 for a 
further discussion of these issues.] 

Such a model, with this relatively moderate number of parameters [especially that, in 
general, only a small subset appears when one looks at a given sector of the model] has 
much more predictability and is much easier to be discussed phenomenologically, compared 
to the uMSSM. We will refer to the MSSM with the set of 19 free input parameters given 
above as the "phenomenological" MSSM or pMSSM^. 



2.1.3 mSUGRA: the constrained MSSM 

All the phenomenological problems of the unconstrained MSSM discussed previously are 
solved at once if one assumes that the MSSM parameters obey a set of boundary conditions 
at the Unification scale. These assumptions are natural [but not compulsory, see Ref. [3~2~| 



e.g.] in scenarii where the SUSY-breaking occurs in a hidden sector which communicates 



4 Note, however, that at the time being the program SUSPECT which will be discussed later does not 
use in the option pMSSM, exactly the same set of input parameters as proposed above. The underlying 
physical assumptions are nevertheless identical [see section 3.2]. 
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with the visible sector only through gravitational interactions. These unification and 



universality hypotheses follows 
• Gauge coupling unification: 

a^Mu) = a 2 (Mu) = a 3 (Mu) = a v (2.15) 
with oti — 9i/^ 1T - Strictly speaking, this is not an assumption since these relations 



are verified given the experimental results from LEP1 ||33|| . In fact, one can view 
these relations as fixing the Grand Unification scale Mjj. 

• Unification of the gaugino masses: 

M^Mu) = M 2 {M V ) = M 3 (M V ) = m 1/2 (2.16) 

Since the gaugino masses and the gauge couplings are governed by the same renor- 
malization group equations, the former at the electroweak scale are given by: 

M, = <M0l mi/2 - M 3 (M Z ) = ^AMM Z ) = ^\M 1( M Z ) (2.17) 
ctu ol 2 [M z ) Oii{M z ) 

For instance, one has the well-known relation: Mi = ^tg 2 9wM 2 ~ \M 2 . 

• Universal scalar [sfermion and Higgs boson] masses 

M Q = Mz R = M 3r = M- l = M Ir = M Hu = M Hd = m (2.18) 

• Universal trilinear couplings: 

A u (Mu) = A d (M v ) = Ai(Mu) = A (2.19) 

Besides the three parameters m^jWo an d the supersymmetric sector is described 
at the GUT scale by the bilinear coupling B and the Higgs-higgsino mass parameter 
//. However, one has to require that electroweak symmetry breaking takes place. This 
results in two necessary minimization conditions of the Higgs potential [see next section 
for details]. The first minimization equation can be solved for the second equation 
can then be solved for B. Therefore, in this model, we will have only four continuous and 
one discrete free parameters: 

tan/3 , mi/2 , m , A , sign(^) (2.20) 

This model is clearly appealing and suitable for thorough phenomenological and experi- 
mental scrutinity. This constrained model, is usually referred to as the minimal Super- 
gravity model, or mSUGRA. In addition, one can also require the unification of the top, 
bottom and tau Yukawa couplings at the GUT scale p4]| . This would lead to a further 



constraint if the "fixed point" solutions are chosen. Depending on whether one includes or 
not the top Yukawa couplings, the parameter tan/5 should be either small tan/3 ~ 1.5, or 
large tan/3 ~ 50 |35|]. The values taken by the parameter A happen to be also constrained 
in this situation which further reduces the number of free parameters. 
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2.1.4 The MSSMi: the intermediate MSSMs 



mSUGRA is a well defined model of which the possible phenomenological consequences 
and experimental signatures have been widely studied in the literature. However, it should 
not be considered as THE definite model, in the absence of a truly fundamental description 
of SUSY-breaking. Indeed, some of the assumptions inherent to the model might turn 
out not to be correct. In fact, in many models some of the universality conditions of 
mSUGRA are naturally violated; see e.g. Ref. [|32], [57| for a discussion. 

To be on the safe side from the experimental point of view it is therefore wiser to 
depart from this model, and to study the phenomenological implications of relaxing some 
defining assumptions. However, it is desirable to limit the number of extra free param- 
eters, in order to retain a reasonable amount of predictability when attempting detailed 
investigations of possible signals of SUSY. Therefore, it is more interesting to relax only 
one [or a few] assumption at a time and study the phenomenological implications. Of 
course, since there are many possible directions, this would lead to several intermediate 
MSSM's between mSUGRA and pMSSM, denoted here by MSSMi's [with i an integer 
and finite, although possibly large, number]. Some of these MSSMi's are similar to the 
Minimal Reasonable Model discussed recently []36|| . A partial list of possible MSSMi's can 



be as follows [many other possibilities are of course possible including the relaxation of 
two assumptions at a time and the introduction of an amount of CP-violation]: 

(1) MSSMI: mSUGRA with no sfermion and Higgs boson mass unification: 

The Higgs sector of the MSSM can be studied in a relatively model independent way, 
since at the tree-level only two input parameters are needed: tan/3 and one of the Higgs 
boson masses. Although the squark mass parameters and the trilinear couplings of the 
third generations, as well as the parameter \x enter through radiative corrections, one can 
study their impact in a thorough manner without invoking any strong assumption [c.f. 
the LEP analyses]. The mSUGRA model is therefore too restrictive, and to have more 
freedom, one can decouple the Higgs sector for the squark sector by relaxing the equality 
of the sfermion and Higgs boson masses in eq. (2.18): 

m g = m u R = rrid R = m L = m e R ^ M Hu = M Hi (2.21) 

Different sfermion and Higgs boson masses are in fact suggested by some SUSY-GUT 
theories e.g. based on SO(10) p7j. One would then have an additional parameter since 
one of the Higgs masses for instance will remain free. 

(2) MSSM2: mSUGRA without sfermion mass unification: 

Scenarii with light stops are appealing in several respects; for instance, a light stop 
with a mass of the order of 100 GeV might trigger Baryogenesis at the electroweak scale 
| Jj8[j . However, it is rather difficult, with the scalar mass unification assumption to have a 
light stop while the other squarks remain rather heavy [without e.g. large A t values] . One 
can then disconnect the third generation from the first two ones [where the experimental 
constraints on FCNC are most stringent] and allow for non universal scalar masses in the 
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third generation squarks [and similarly for sleptons] 

m Q ^ m R + m b R ^ m L^ m r R ^ rrif (2.22) 
with mj the (common) mass parameter of the first /second generation sfermions. 

(3) MSSM3: mSUGRA with no gaugino mass unification: 

The chargino/neutralino sector depends only on three parameters: /i, M 2 and tan/3 
if gaugino unification is assumed. This allows to make thorough experimental analyses 
which led to important results, such that the mass of the lightest neutralino should be 
larger than ~ 40 GeV [see section 6]. Furthermore, searches for charginos at LEP2 and 
gluinos at the Tevatron are connected since the gaugino masses are related. One can go 
one step downwards in model-dependence and relax the gaugino mass unification: 

M^Mu) £ M 2 (M V ) ^ M 3 (Mu) (2.23) 

This would make the connection between the chargino, neutralino and the gluino sectors 
less strong and will e.g. leave open the possibility of very light neutralinos. 

(4) MSSMO: pMSSM with partial mass and coupling unification 

This model, with 7 free parameters, is the most used in phenomenological analyses: 

tan/3 , Ma , // 
M^Mu) = M 2 (Mu) = M 3 (Mu) = m 1/2 , A t = A b = A T = A 
m Q = m i R = m b R = m <i = m u R = m iR , mi = m fR = mj = mj = m- eR (2.24) 

2.2 Electroweak Symmetry Breaking 
2.2.1 General features 

We turn now to the discussion of the electroweak symmetry breaking (EWSB). Using the 
notations introduced in the previous section, the Higgs potential takes the form 

^Higgs = Vtree + V\ 

Vtree = (m^ + /i 2 )if]# d + (m^ + // 2 )^F U + Bfx(H u .H d + h.c.) 

+ 9 j(HlH d - HlH u f + 9 f(HlH u )(HlH d ) (2.25) 

where g 2 = g\ + g\\ in some cases we will also use the shorthand notations: 

m i = m H d + V 2 > m l= m H u + V 2 > m l = - B V ( 2 - 26 ) 
V\ contains the higher order corrections to this potential. Here we consider solely the 



one-loop corrections which have the well-known form in the MS scheme 39 



2 



ft M' 
V, = Str[M 4 (Log— - 3/2)] (2.27) 
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where /i# denotes the renormalization scale, M 2 the field dependent squared mass matrix 
of the scalar or vector or fermion fields, and Str[...] = S S pin( — l) 2s (2s + 1) (...)„, where the 
sum runs over gauge boson, fermion and scalar contributions. 

Even if no model assumption apart from minimal Supersymmetry is made, i.e. no 
unification of the gauge couplings, no universality of the soft-SUSY breaking terms, no 
Yukawa coupling unification, etc... it is clearly important to still require the electroweak 
symmetry breaking to take place. The usual necessary condition for EWSB is obtained 
from eq. (2.25) as 



1 , ,9 m 2 - ml tan 2 3 . „ n —2ml 
2 z tan 2 3-1 m\ + m 



M| = L — - 7 i — ^— - , sin 2(3 _., (2.28) 

2 



with M§ = g 2 /4(t> 2 + v 2 ,) and tan/3 = v u /vd, where we assumed 

<H d >= ±( V *) <H U >= -L( ; 1 (2.29) 




Eqs. ( |2.28| ) contain two complementary and necessary requirements, for i) the breaking of 
the electroweak symmetry, ii) the Z mass value to be reproduced correctly through this 
breaking. These two conditions are however generally not sufficient to ensure electroweak 
symmetry breaking. For one thing, beyond tree-level they only express the existence 
of a stationary point not necessarily a global (nor even local) minimum of the effective 
potential. In a phenomenological analysis one then usually checks numerically for the 
globality of the minimum [see next section for an analytical approach]. The other reason 
is the possible existence of color or charge breaking minima which can be either lower than 
the electroweak minimum or sufficiently stable to become dangerous from a cosmological 
point of view [we will have, however, nothing to say about these minima in the present 
report] . 

Thus it should be clear that even in the unconstrained MSSM one should at least 
impose the constraints of eq. Q2.28 ). These equations correlate not only the parameters 



of the Higgs sector, but actually all the other parameters of the model when the radiative 
corrections to Vmggs are taken into account. Then one has typically 

m\ = m 2 Hd + /i 2 + rad. corr. (v u , v d ) 
m~l = m\ + /i 2 + rad. corr. (v u , Vd) 

ml = -Bfj, (2.30) 
In this case eqs. (2.28) are no more quadratic in tan 3 and linear in sin 23 [not even 



polynomial in these variables anymore] so that one usually resorts to numerical methods in 
solving these equations. At this level we should perhaps restate a question of terminology. 
What is usually called radiative electroweak symmetry breaking is the fact that eqs. (|2.28|) 
are satisfied at a given scale [presumably the electroweak scale] through the running of 
the quantities which enter these equations, down from a high scale where some initial 
conditions were assumed. Relaxing the radiative breaking requirement simply means 
that one no more insists on starting from a high scale and specifying initial conditions, 



17 



but requires directly the electroweak symmetry breaking. This would be a fully model- 
independent but still a physically consistent approach. Of course one could eventually run 
the masses and couplings up to a high scale to assess their consistency with any model 
assumption at that scale. 

Finally one can also implement eqs. (|2.28|) indirectly. For the sake of illustration we 



give here a tree-level example: the usual tree-level Higgs boson mass relations 

Kb = \i M z + M lT y/(M% + M\f - AM%M\ cos 2 2(3 ) 

M 2 H± = .U-; • .Uf r (2.31) 

which are a consequence of the special form of the tree-level part of the Higgs potential 
eq. ( |2.25| ) and of eqs. ( |2.28| ). The four free parameters m^m^m 2 and tan/3 reduce to 
two free parameters due to the constraints of eq. (|2.28|) , usually chosen as Ma and tan/? 



in eqs. ( |2.31D 



However one can show explicitly that eqs. ( [2.3 1|) do not necessarily imply electroweak 



symmetry breaking i.e. using expressions with non-zero masses is not sufficient to ensure 
vacuum stability. Such mass relations can be realized and still have eqs. ( [2.28| ) violated. 
Only if one imposes on top of these mass relations Ma to be given by 

M\ = m\ + m\ (2.32) 

that EWSB is assured. Also this result is established at tree-level and it is not clear 
whether it still holds when loop corrections are taken into account. Thus it is always safer 
to check explicitly eqs. ( [2.28| ) whenever possible. 

2.2.2 EWSB and model— independent tan/3 bounds 

In this section we report on some analytical results for model independent theoretical 
bounds on tan f3 obtained from the study of electroweak symmetry breaking conditions to 
one-loop order [fj, [|. The point is to use, on top of eqs. ( |2.28|) , the positivity of the Higgs 



boson squared masses which are needed to ensure a, at least local, minimum in the Higgs 
sector. We then study analytically how these constraints translate into bounds on tan /3 
which are then necessary and fully model-independent bounds on this parameter. We re- 
call here that the positivity of the squared Higgs boson masses is automatically satisfied 
as a consequence of eq. fl2.28|) at the tree-level [or tree-level renormalization group im- 



proved] effective potential. However, this property is not expected to be generic when the 
one-loop corrections are taken into account [beyond the ones included in the runnings]. 
This was explicitly shown in Ref. in a specific approximation. Since the full one-loop 
effective potential has a complicated form we relied on two different approximations and 
showed that they both lead to qualitatively similar results. The first of these approx- 
imations consists in absorbing all logarithms of the one-loop effective potential in the 
runnings of the tree-level quantities thus casting Vmggs in the form, see eqs. ( |2.25| , [2.27j) , 



Vkiggs = ^tree(/4) + ("| ) StrM 4 (2.33) 
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From now on we refer to this approximation as the Supertrace approximation. Here 
^tree ^ s obtained from V tree by replacing all the tree-level quantities by their running 
counterparts. This would be fully justified in a model with just one mass scale and would 
mean that we have resummed to all orders the leading logarithms in the MS scheme. 
Note that the residual one-loop correction in eq. (|2.33| ) is scheme dependent but should 
be consistently kept as a residual correction [and not reabsorbed in the runnings as it 
is sometimes suggested] since it would otherwise jeopardize the resummation procedure; 
see for instance Ref. PP[ . Of course the MSSM has many different mass scales and the 
above approximation is therefore very rough. It has however the merit of allowing a 
full analytical proof of the existence of bounds on tan j3 free from any phenomenological 
assumption. This is significant in the sense that our approximation with one mass scale 
tends to increase the symmetry of Vniggs so that if new bounds appear because of the 
difference in structure between tree-level and one-loop, then these would hardly disappear 
in more realistic, and less symmetric, approximations. 

The analysis will not be carried further here, the interested reader is referred to Ref. 
for full details. Hereafter we only give the generic form of the tan/5 bounds and then 
discuss briefly how some of these bounds can arise also in a more realistic approximation, 
namely the top/stop-bottom/sbottom approximation. 

In the Supertrace approximation, the bounds on tan (3 read: 

if tan/3 > 1: tan/3_ < tan/3 < Min(tan/3+, -X 2 m jX 2 mz ) (2.34) 



where 



tan /3_ = Min(T + , m t /m^) and tan (3 + = Max(T + , m t jm^) (2.35) 
if tan/3 < 1: Max(T_, -X^JX^J < tan/3 < 1 (2.36) 



where 



2X 2 

mi 



-X 2 mi - X 2 m2 =F J(X 2 mi + X 2 m2 f - 4X^ 3 (2.37) 



The X^. are generalizations of the m 2 's which include residual one-loop corrections from 
eq. ( |2.33| ). Note here that these bounds are slightly improved with respect to the ones in 



Ref. p| due to the presence of —X^/X^ and — X^/X^. Since we choose by convention 



tan/3 > 0, the positivity of M\ translates into X^ a < 0, while the signs of X^ and X^ 2 
are not fixed 



We come now to the top/stop-bottom/sbottom approximation and give here a little 
more details about our approach. In this approximation, the one-loop correction to the 
Higgs effective potential is approximated by 

rr 6/i r _ 4 / ml 3\ 4 / m 2 3 

Vi = ^ ^< L °^-2j- 2m *v L0g[ ^ ] -2 



64tt 2 
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^2.38) 
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Taking a linear combination of the two stationarity conditions of Vm ggs at the electroweak 
minimum, and defining u = v u Vd and t = tan f3, one finds 



T = K 2 {u)u 2 + ld(u)u + Ko(u)t = 



(2.39) 



where 



K 2 {u) = ^{t 2 - I) 2 
v ' 96 v ; 



with 
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(2.42) 



3 



x(rr& - m 2 3 ) - -[g 2 (t 2 - 1) - 8Y t V](A t * - - fit) 

2 Log[(m9)V(m9) 



K) 2 - K 2 ) 2 



6y b 2 [(A 6 t - aO(a 6 - /it) + 2t(A - fit) 2 } 
3 



- 9l)(t 2 ~ !)*("& - m}) + -g 2 (t 2 - l)(A h t - fi)(A b - fit) 



2 WZ J1/y ~ " /-v " Q3 
The last term in % takes the following form: 



/Co(u) = -(1 + t 2 ) A°H + t^ ( M ) 



(2.43) 



where: 
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and 
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[with k = ^/(647r' / 

In deriving the above formulae, the dependence on the Higgs fields in the logarithms 
was fully taken into account, together with the convention > m~ 2 ,m 2 > m 2 ^. The 
key point now is to note that on one hand 



and on the other hand 



K. 2 (u) > 



/Ci(«)t > 



(2.46) 



(2.47) 



The first inequality is obvious from eq. ( |2.4(J| ), while the second requires some mild phe- 
nomenological assumptions to overcome the analytical complexity of eq. ( |2.41| ) Let us give 
here just two examples of such assumptions, leaving a more detailed study to Ref. ||. (?) 
If the mixing between the stop eigenstates is weak, A t tan/3 -/i< 1, and neglecting the 
gauge contributions to JCi, one easily sees that /Ci is dominated by /C^ which behaves 
like m\. Since in this limit the stops are almost degenerate and heavier than the top 
[assuming all squared soft masses to be positive], then eq. (|2.47|) is readily verified, (ii) A 
second example of mild phenomenological assumption, is to take the heaviest stop mass 
larger than ~ 360 GeV, and tan/3 small ( > 1). Combined with the experimental lower 
bound on the lightest stop [~ 80 GeV, see section 6], one again obtains eq. ( [2.47| ). 



K t < (2.48) 

The positivity of M\ is not automatic beyond the tree-level, and should be imposed 
explicitly] one retrieves a generalization of the T± bounds of the previous approximation. 
With the convention tan /3 > 0, which we did not need to impose in the previous discussion, 
the generalized T± bounds read: 



for tan (3 > 1 , tan /3 > T+ (2.49) 
for 0<tan/?<l , tan/3 < T_ (2.50) 



with (T+ > T_) 



-/Co ± V(M 2 " (A^) 2 



T ± = v A U o (2.51) 
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In summary, we have determined analytically and in a model-independent context, calcu- 
lable bounds on tan/5. These bounds are nothing but partial necessary constraints coming 
from the requirement of electroweak symmetry breaking in the MSSM. Since these nec- 
essary bounds are calculable in terms of the parameters of the MSSM, they can be used 
to delineate domains which are theoretically inconsistent with EWSB. This would be a 
valuable guide in the standard procedure of the numerical check of electroweak symmetry 
breaking, where one can avoid from the start inconsistent domains in the input param- 
eter space. Furthermore, we emphasize that the above bounds have in principle a wider 
applicability, and are more quantitative, than the usual bounds 1 < tan/5 < mt/m^. 

We show in Fig. 1 the sensitivity of T + to the mixing parameter A t = A t tan j3 — /i 
in the stop sector. This illustrates the amount of exclusion of tan/5 depending on the 
approximation used. 




Figure 1: The T + model-independent lower bound on tan/5 in the top/stop- 
bottom/sbottom approximation and in the Supertrace approximation. Here: M\ = 
m\ + m\ = (300 GeV) 2 , m 2 3 = -(100 GeV) 2 , Q = m t = 175 GeV, M SU sy = 630 
GeV, u.± = ±100 GeV and Y t = 0.9 
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2.3 Renormalization Group Evolution 



In the minimal SUGRA model, the MSSM parameters [couplings and masses] are defined 
at the Unification scale with some unification conditions, and then are evolved down to the 
electroweak scale through Renormalization Group Equations (RGE). The RGE evolution 
is thus an important ingredient of mSUGRA, and more generally in any theory incorpo- 
rating GUT unification. In the RGE's of the MSSM, different levels of approximations 
are available as will be discussed later. In this section, we will write for completeness the 
RGE's for the masses and couplings in the one-loop approximation, that we will need 
later when we will discuss the program SUSPECT. 

In a next step, we will discuss some exact analytical solutions of the one-loop RG 
evolution equations of the Yukawa couplings. Solutions to these equations exist in the 
literature since many years [|2[ in the the limit where the top Yukawa coupling is assumed 
to dominate all the others. This limiting solution is relevant for small tan/5 values and 
remain numerically useful for tan/3 values up to 10 or so. Later on, various attempts 
were made to obtain general solutions, but still relying on some approximations, such as 
neglecting the U(l) coupling as compared to those of SU(3) and SU(2). However, even in 
this case, the solutions given where actually implicit |Q and not directly generalizable 
to include more than two Yukawa couplings (Y t , Y&). Below, we will give the exact solu- 
tions with no approximation whatsoever. They are valid for any value of Y t and Yj, and 
immediately generalizable to the full set of Yukawa couplings. 



2.3.1 The one-loop RGE's 



In the following, we list the RGE's for the MSSM parameters in the one-loop approxima- 
tion. Although two-loop evolution equations for some parameters [such as the gauge and 
the Yukawa couplings] are available, for many purposes it is a rather good approximation 
to use only the one-loop equations [especially that it is much faster when implemented in 
numerical programs]. The list that we give below is ordered as in the program SUSPECT 
which will be discussed in section 3.2. 



Gauge Couplings [<?j with i — 1, 2, 3 and n g the generation number]: 
dg t 1 



-b i g i with bi 



1 10 h 
~ ~3~ n9 ' 



dt 32?r 2 

Yukawa Couplings [i = 1,2,3 generations]: 
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32vr 2 [ 


dYj 


yi r 
1 d 


dt 


32tt 2 [ 


dY? 


Yi r 


dt 


327T 2 . 



13 2 o 2 16 2 



E{3(^ fc ) 2 + (Yf) 2 } + ra 2 + my 

k 

£{ W) 2 + 3(^ fc ) 2 } + m? - 3(g 2 + gl) 



16 



g9 2 + 3#2 + y al 



2.53) 
2.54) 
2.55) 
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The \i term and the vacuum expectation values vi, v 2 [k = 1, 2, 3]: 

djj, jj, 
~dt ~ ~32^ 
dvi v\ 



dt 
dv 2 
~dt 



327T 2 
V2 

32tt 2 



E{3(^) 2 + (^) 2 }-^(^ 2 + ^ 2 ) 



(2.56) 
(2.57) 
(2.58) 



The scalar Higgs masses and the parameter B [k — 1, 2, 3]: 



dm 2 Hi 

dt 
dm 2 H2 

dt 
dB 

~dt 



16tt 2 
1 

167T 2 

1 

lQ^ 2 



£{3(1?) 2 P| + (>f ) 2 Pft - - 9 lTr(Ym 2 ) - {g\M\ + 3g 2 2 M&.59) 



£3(Y U *) 2 P| + ^(Ym 2 ) - {g\M\ + 3g 2 2 M 2 ) 
k A 

£{ 3 ^) 2 + 3^ ) 2 + A^) 2 } - (gfMi + 3^ 2 M 2 ) 



(2.60) 
(2.61) 



The trilinear A couplings [i, k — 1, 2, 3]: 



d^ 
dt 



dAj 



32tt 2 



64(yj) 2 + 24(yi) 2 + 6E^(^) ; 

fc 

/26 32 



327T 2 



6A*(yj) 2 + 24«) 2 + 2^{A k (Y l k ) 2 + SA k d (Y d k ) 2 } 



14 



32 



9 giM 1 + 6g z 2 M 2 + —g z 3 M 3 



QA\(Y?f + 2 £{Af (if) 2 + 3^(F d fe ) 2 } - Q{g\M x + g 2 2 M 2 ) 
k 



(2.62) 

(2.63) 
(2.64) 



dt 32vr 2 

The scalar fermion masses [with Pt - r = m%„ „ „ n + m\ A f + m 2 ? r +(A A iw/) 2 l 

^ u,«M H 2 ,H 1 ,H 2 Q k ,Qk,L k u k ,d k ,l k V «,<W J 



dm 2 

l *i _ 


1 ■ 


dt 


16tt 2 [ 


dm 2 f 


1 ■ 


dt 


16tt 2 [ 


dm 2 -: 


1 ■ 


dt 


167T 2 [ 


dm l Rl 


1 ■ 


dt 


16tt 2 [ 


Am \ 


1 ■ 


dt 


16tt 2 [ 



2(Yi) 2 P{ + g 2 Tr(Ym 2 ) - 4gfM 2 



{Y?ypt - -^Tr(Fm 2 ) - (gfMf + 3^ 2 M 2 ) 



2(^) 2 P] + Wm a ) - ( ^ 2 M 2 + ^ 3 2 M| 



1 



! • 3 
2 



16 



9' 
/16 



16 



2(^) 2 ^ " §^ 2 Tr(Fm 2 ) - (y^M 2 + ^g 2 3 M 2 

ra 2 ^ + ra 2 ^]+^i 2 Tr(rm 2 ) 



(2.65) 
(2.66) 
(2.67) 
(2.68) 
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1 16 

~glM? + 3g 2 2 Ml + -g 2 3 M* 



The gaugino masses [i = 1,2,3 and the bi are given above] 



dMi 
~df 



16n'- 



■MAg- 



(2.69) 



^2.70) 



A few remarks are in order, here: 

(i) The evolution parameter t is defined by t = Log(M^/Q 2 ); this is different from 
the one used in the RGE's of the program SUSPECT, where t = l/2Log(Q 2 /M^) 



Tr(Ym 2 ) is the isospin pondered sum of the squared soft masses of the scalar fermions; 
in the case of universal soft masses, the trace vanishes at any scale due to anomaly 
cancellation, (in) The RGE's for the gaugino masses and the gauge couplings are related 
and from eqs. (|2.52| ) and (|2.70|) , one can easily see that d/dt(Mi/g 2 ) = 0. 



2.3.2 Exact solutions for the Yukawa coupling RGE's 



We are interested here in eqs. (|2.53| - ^755D . They have the nice feature of being completely 
decoupled from the rest of the system, especially from the gauge couplings whose running 



is determined a priori via eq. (|2.52|) . [This is no more true at two-loop order where the 
gauge and Yukawa equations become interwound.] When all Yukawa couplings except 
Y t are neglected, eqs. ( |2.54|) and (|2.55|) become trivial while eq. (|2.53|) becomes of the 
Bernoulli type in the variable y t = Y 2 



-r t Vt = h(t)Vt 



by 2 



where 



_L (-^ 2 + 3 g 2 +- g 2 ) 
16vr 2V 3 y3 y2 9 yi) 



16tt 2 



and is easily solved to give [pB| , ^2 | 

Vt{t) = - 
where 

E(t) = e Io^ dt ' and y° = Y 2 {t = 0) 



y°E(t) 

by°J*E(t>)dt> 



(2.71) 



(2.72) 



(2.73) 



(2.74) 



In the more general case where both Y t and Yj, are kept in the game, but neglecting all other 
Yukawa couplings, eqs. (|2.53| , |2.54]) become after the change of variable, y t = Y 2 ,y b = Y 2 



d 

dt Vt 
d 

dt Vb 



fi{t)y t + ay b y t + by 2 
h{t)Vb + aybyt + by\ 



(2.75) 
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where fi(t) and b are given in eqs. (|2.72|) and 



Ut) = J—(— g 2 + 3gl + -g2) , a = l — (2.76) 

J v ; 16tt 2V 3 ys y2 9 yiJ 16tt 2 v ; 

As far as we know, the system eqs. ( |2.75| ) is not treated in standard text books, and 
although it looks at first sight simple, we could not find a systematic way of relating it to 
a standard formQ. It is also relatively easy to solve the system up to first order in Yj, in the 
region Y t ^> Yj, Q. This is already an improvement of the known solutions with Yf, ~ 0. 
It extends the numerical validity much further than tan f3 ~ 10. More importantly, this 
approximate solution gave us a hint of the structure of the exact solution which was then 
found by sheer guess Hj: 



l-by?f*E 12 (t')dt> 

y° b E 21 (t) 

l-by° b foE2i(t>)dt> 



7/ - ™ (9 77\ 



- gtfglW (9 78^ 

Vb - ; — r. , r t n ^jjj \ 2 -™) 



where 



El2{t) ~ (i-^S(w (2,79) 
4(t) = (i-^S(w (2,80) 

E . = e f h{t')dt' i = 1?2 (2.81) 



and |/° = 3^ 2 (t = 0),y° = V 6 2 (t = 0) are arbitrary initial conditions. The solutions 
eqs. ( |2.77| , |2TT8"| ) are exact for any value of tan f3. They resemble formally eq. ( |2.73| ) of which 
they are a generalization. One should note however the important difference, namely that 
due to eqs. ( ^i.79| . |2.80"D , the general solutions for y t and ?/& are actually continued integrated 
fractions. Indeed eqs. fl2.79U2.80D give an implicit definition of Ei 2 (t) and E 2 i(t), the first 



being defined in terms of the second and vice- versa. One could still write E\ 2 (t) in an 
explicit though infinite series form: 

El2(t) = n hv Q[t ) a ~ h (2 - 82) 

^ °% Jo o r*2 g 2 (t 3 jdt 3 /b > 

and similarly for E 2 i(t) with the substitution 1 <-> 2, y° «-> We will see later on that 
both forms for E^- are useful. In any case, we should stress here that the solutions for 

5 The situation would be much simpler if f\(t) = f2{t), i.e. when neglecting g\. In this case the 
equations can be solved in quadrature after some change of variables, leading though only to implicit 
solutions involving some hypergeometric functions f43] 
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yt,yb are themselves explicit. 

What do we gain from these exact analytical solutions? 

(z) First of all one can prove rigorously the convergence of the infinite series 
and determine explicitly the convergence criteria [|J. This implies that, for practical 
purposes, keeping only the first iteration of the Eij(t) series is numerically a very good 
approximation [within the convergence region] . We give in Table 1 a numerical comparison 
versus a Runge-Kutta method. 

(ii) The large tan (3 region is treated exactly and one can follow precisely the various 
features of the running of the Yukawa couplings in this regime. 

(Hi) The fact that the coefficients of the quadratic parts of eqs. ( |2.75|) are equal and 



that there are only two non-zero Yukawa couplings is actually unessential in finding the 
general solutions in the present form. This form generalizes straightforwardly for an 
exact solution of eqs. (|2.53| - |2~55|) which will be given elsewhere , and would thus be of 



relevance in the case of bottom-tau Yukawa coupling unification scenarios B4 



(iv) Finally, one can control analytically the acceptable regions for the initial condi- 
tions y®. This is related to the necessity of avoiding Landau poles and more generally 
of requiring y t , to remain positive throughout the evolution, being the squares of Y t , Y^. 
This is of relevance if one wants to run the quantities between two low-energy scales 
choosing some initial conditions at one of these scales without referring explicitly to the 
initial values at the unification scale Mjj. [Such a possibility is being implemented in the 
fortran code SUSPECT, see section 3.2]. 



tan (3 


Y b (t = 0) 


Y t (t = 0) 


Y b (t) 
"exact" 


Y b (t) 
Runge-Kutta 


Ytit) 
"exact" 


Ytit) 
Runge-Kutta 


2 


0.0387453 


1.13007 


0.0145059 


0.0145050 


0.775788 


0.775974 


10 


0.174138 


1.01581 


0.0630978 


0.0631052 


0.54263 


0.542743 


50 


0.866544 


1.01097 


0.435682 


0.439526 


0.585453 


0.590258 



Table 1: Numerical comparison of the exact one-loop solution, truncated to the first 
iteration, with the Runge-Kutta RG evolution, obtained with only a one step evolution 
over 10 orders of magnitude. 
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3. The Physical Parameters 



3.1 Particle masses and couplings 
3.1.1 Mass matrices and couplings 

In this section, we discuss the general features of the chargino/neutralino, sfermion and 
Higgs boson sectors to set the conventions and the notations which will be used further on. 

a) The chargino/neutralino sector 

The general chargino mass matrix depends on the parameters M 2 ,/i and tan/3 |22j, ^5 



M 



c 



M 2 V2M W sin (3 

y/2Mw cos (3 [i 



(3.1) 



is diagonalized by two real matrices U and V, 
WMcV- 1 -> U = £L and V = 



0+ if detMc > 
a 3 0+ if det.M c < 



(3.2) 



where the 03 matrix and the 0± matrices are given by [with the appropriate signs de- 
pending on the values of M 2 , /i, and tan/?] 



^3 



with 



+ 1 
-1 



tan 26 '_ 
tan 20+ 



cos 6± sin 9± 
— sin 6± cos 6± 



(3.3) 



2 V / 2M W /(M 2 cos (3 + n sin (3) 

Mf -fi 2 -2M$ v cos ]} 
2\f2M w {M 2 sin (3 + \i cos (3) 



Ml -/i 2 + 2M^ cos (3 



(3.4) 



This leads to the two chargino masses, 
1 

V2 



Ml + /i 2 + 2M^ 



(3.5) 



=F {{Ml - /i 2 ) 2 + 4M^ cos 2 2(3 + 4M^(M 2 + /i 2 + 2M 2 /isin2/3)}- 

In the limit |/i| ^> M 2 , Mz, the masses of the two charginos reduce to [e M = sign(/i)] 

Mi 



~ Mo — 



X-2 



M 2 W 



(M 2 + /isin2/3) 



|Ml + — r e » (M 2 sm2(3 + fi) 



(3.6) 
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For I //. I — > oo, the lightest chargino corresponds to a pure wino state with m + ~ M 2 , 
while the heavier chargino corresponds to a pure higgsino state with m v + = \u\. 

In the case of the neutralinos, the four-dimensional neutralino mass matrix |22] de- 
pends on the same two mass parameters /1 and M 2 , if the GUT relation M\ = | tan 2 8w 
M 2 ~ §M 2 is used. In the (—iB, —iW 3 , H®, H 2 ) basis, it has the form [sly 
sin 2 # 



\v\ 



M 



N 



M x 
M 2 
-M z sw cos (3 M z cw cos /3 

M^s w sin (3 —M z c w sin /3 



-M z s w cos (3 M z s w sin /? 
M z cvk cos /3 — M z cvk sin (3 
-/i 




(3.7) 



It can be diagonalized analytically by a single real matrix Z\ the [positive] masses of 
the neutralino states m v o have complicated expressions which will not be given here. In 
the limit of large values, the masses of the neutralino states however simplify to 



m v o 



A2 



m v o 

A3 



Mi 



/i 2 

M 2 -^ 
/x 2 

1M 2 



(M x + //sin2/3) s 2 w 
(M 2 + /x sin 2/3)4- 



2 /i 2 



1M le (1 
2 t ^v 1 



2 /x 2 



e A*(l — sin 2(3) (ja + M 2 s 2 w + M lC 2 w 
sin 2(3) U - M 2 s 2 w - M lC 2 w 



Again, for — * 00, two neutralinos are pure gaugino states with masses m x o ~ Mi 



m x o = M 2 , while the two others are pure higgsino states, with masses m x o ~ m x o ~ 
b) The sfermion sector 



Assuming a universal scalar mass mo and gaugino mass at the GUT scale, one 
obtains relatively simple expressions for the left- and right-handed sfermion masses when 
performing the RGE evolution to the weak scale at one-loop order, if the the Yukawa 
couplings in the RGE's are neglected [for third generation squarks this is a poor approx- 
imation since these couplings can be large; in this case numerical analyses are needed as 
will be discussed later]. One has: 



m 2 f 

JL,R 



m 



I + J2F i (f)m 2 ± (Ij - e f s 2 w )M 2 z cos2(3 



(3.9) 



i=i 



If and e f are the weak isospin and the electric charge of the sfermion and Fi are the RGE 



coefficients for the three gauge couplings at the scale Q ~ M Z} given by 

12 \ 



Ci(f) 



1 h aU hJ QZ 

1 - I 1 - — hLog— j 



4:71 



'Mi 



u , 



(3.10) 
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The coefficients b i} assuming that all the MSSM particle spectrum contributes to the 
evolution from Q to the GUT scale M G , are given by: b\ = 33/5, b 2 = l,b 3 = —3. The 
coefficients c(f) = (ci, c 2 , c 3 )(j) depend on the hypercharge and color of the sfermions 



c(L) 



2 

V o 



,c(l 



c{Q) 



2 

V I / 



C{U R ) 




, c(d 



R 



( — 

15 



\ 3 



With the input gauge coupling constants at the scale of the Z boson mass oi\{Mz) ~ 
0M,a 2 (M z ) ~ 0.033 and a 3 (M z ) ~ 0.118, one obtains M v ~ 1.9 x 10 16 GeV for the 
GUT scale and o>u = 0.041 for the coupling constant ajj. Using these values, one obtains 
for the left- and right-handed sfermion masses 





= m l 


+ 6.28m 2 1/2 


+ 0.35Mf 


cos 2(3 




= m 2 


+ 6.28m 2 1/2 


- 0.42M| 


cos 2(3 


771? 


= m 2 


+ 5.87m 2 1/2 


+ 0.16M| 


cos 2(3 




= ml 


+ 5.82m 2 1/2 


- 0.08M| 


cos 2(3 


m k 


= m 2 


+ 0.52m 2 /2 


+ 0.50Mf 


cos 2(3 




= m l 


+ 0.52m 2 /2 


- 0.27M| 


cos 2(3 


m L 


= m l 


+ 0.l5m 2 1/2 


- 0.23M| 


cos 2(3 



(3.11) 

In the case of the third generation sparticles, left- and right-handed sfermions will mix 
P7H ; for a given sfermion f — t,b and f, the mass matrices which determine the mixing 
are given by 



M 



f 



m , 



7l 
m f( A f 



+ m 2 f mj(Af — firf) 



"f 



m ; + mi 

JR 1 



(3.12) 



where the sfermion masses ^j LR are given above, mj are the masses of the partner 
fermions and r& = r T = l/r t = tan/5. These matrices are diagonalized by orthogonal 
matrices; the mixing angles 8f and the squark eigenstate masses are given by 



sin 29 f 
1 



2m f {A f - /j,r f ) 



cos 29 



m ■ 



m ■ 



f 



m- t — m% 

JL JR 



m ; 



m. 



m ,- 



m 



fi h h h 

)l + m ) R T ~ m ) R ^ + Am % A f - ^fY 



(3.13) 
(3.14) 



Due to the large value of mt, the mixing is particularly strong in the stop sector. This 
generates a large splitting between the masses of the two stop eigenstates, possibly leading 
to a lightest top squark much lighter than the other squarks and even the top quark. 



c) The Higgs sector 
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We come now to a description of our parameterization of the MSSM Higgs sector 
23||. Besides the four masses, M h , M H , Ma and M H ±, the Higgs sector is described at 



the tree level by two additional parameters, tan/? and a mixing angle a in the CP- 
even Higgs sector. Due to SUSY constraints as discussed before, only two of them are 
independent and the two inputs are in general taken to be tan/3 and Ma- Radiative 
corrections, the leading part of which grow as the fourth power of the top mass and 
logarithmically with the common squark mass [351 HO], change significantly the relations 
between masses and couplings and shift the mass of the lightest h boson upwards. These 
radiative corrections are very important and should therefore be included in any analysis. 
Here we will, however, only discuss the leading part of this correction which in the simplest 



case can be parameterized in terms of the quantity |48 



3G f 



m; 



V2% 2 sin 2 /3 



Log 1 + 



(3.15) 



The CP-even Higgs boson masses are then given in terms of the pseudoscalar mass Ma 
and tan/5, and the charged Higgs boson mass in terms of Ma, are given by 



M 



h,H 



M\ + Mf + e T \j (Ml + M| + e) 2 - 4M|M§ cos 2 2/3 - 4e(M| sin 2 (3 + Mf cos 2 (5) 



M 2 H± 



Ma + M^ 



(3.16) 



The mixing angle a reads in terms of Ma and tan/? 

M\ + M| 



tan 2a 



tan 2/3 



7T 



Ml-Mf + e/ cos 2/5 



— < a < 
2 ~ ~ 



(3.17) 



Once tan/3 and Ma are chosen and the leading radiative correction is included in a, all 
the couplings of the Higgs bosons to fermions, gauge bosons and Higgs bosons are fixed. 
However, in the trilinear MSSM Higgs boson couplings, there are also large radiative 
corrections which are not entirely mapped into the angle a, but the leading part can also 



be expressed in terms of the leading correction e [|5Q| - 





9&UU 


9<l>dd 


9§vv 


9®VA 


h 


cos a/ sin (3 


— sin a/ cos (3 


sin(/? — a) 


cos(/3 — a) 


H 


sin a / sin f3 


cos a / cos (3 


cos(/3 — a) 


— sin(/3 — a) 


A 


l/tan/? 


tan/3 









Table 2: Higgs boson couplings in the MSSM to fermions and gauge bosons relative to 
the SM Higgs couplings, and coupling to Higgs and gauge bosons. 

The couplings of the charged Higgs boson to down (up) type fermions are (inversely) 
proportional to tan/?, as in the case of the pseudoscalar Higgs boson. For the CP-even 
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Higgs bosons, the couplings to down (up) type fermions are enhanced (suppressed) com- 
pared to the SM Higgs couplings [tan/3 > 1] ; the couplings to gauge bosons are suppressed 
by sin / cos(/3 — a) factors. A has no tree level couplings to gauge bosons. Note also that 
while the couplings of the h and H bosons to ZA and W + H~ pairs are proportional to cos 
and sin(/3 — a) respectively, the W + H~A coupling is not suppressed by these factors. The 
couplings of the MSSM neutral Higgs bosons to fermions and gauge bosons [normalized 
to the SM Higgs coupling g Hsmff = {^2G F fl 2 m s and g HmVV = 2( V / 2G F ) 1/2 M 2 ] and 
to gauge and Higgs bosons [normalized to (a/2Gf) 1 / 2 /M^(p$ + Pa)^ with p$ and pa the 
Higgs bosons 4-momenta] are given in Table 2. 

Let us turn now to the h boson couplings to stop squarks which will be a very important 
ingredient for the discussion of the next section 4. Normalized to 2M§(v2G f) 1 ^ 2 , they 
are given in the decoupling limit Ma 3> Mz, by 



9hi 2 t 2 



cos 2/3 



cos 2/3 



cos 6+ — -s w cos 26 



— sin Qi — -s w cos 26, 



m t 
Ml 

2 

m t 

Ml 



1 „ m+Af 

1 • OQ mtAt 

sin 26 f - ^ 

2 1 Ml 



(3.18) 



and involve components which are proportional to A t = A t — \ij tan (3. For large values of 
the parameter A t which incidentally make the t mixing angle maximal, | sin2^-| ~ 1, the 
last components can strongly enhance the g h n couplings and make them larger than the 
top quark coupling of the h boson, g htt oc m t /M Z - The couplings of the heavy H boson 
to stops involve also components which can be large; in the case of the lightest stops, the 
coupling reads in the decoupling limit: 



9uhh = sin2 P 



1 2 

-cos 2 Of - -s^ cos 26\ 



Ml tan/3 



+ \ sin 26, 



mt , A t 



Mf 



tan/3 



+ H) (3.19) 



For large tan/3 values, the mf and the A t components are suppressed; only the component 
proportional to /i is untouched. The pseudoscalar A couples only to iit 2 pairs because of 
CP-invariance, the coupling is given by: 

1 m t 



2 M 2 Z 



(Ai/tan/3 — fi) 



(3.20) 



In the maximal mixing case, | sin 26 f \ ~ 1, this is also the main component of the H boson 
coupling to tit2 pairs except that the sign of /i is reversed. 

Finally, let us exhibit the couplings of the h boson to sneutrinos and the lightest 
neutralinos that we will need in section 4.2. For sneutrinos, the couplings is just given in 
eq. (|3.18| ) and setting the fermion mass and the mixing angle to zero: 



9hc 



— cos 2/3 
2 H 



For the normalized couplings of the h boson to the LSP, one has 

Shxlxl = ( Z w ~ tan QwZii) (sin /3Z U - cos/3Z 13 ) 
with Z is the matrix diagonalizing the neutralino mass matrix discussed above. 



(3.21) 



(3.22) 
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3.1.2 Inverting the chargino/neutralino spectrum 

Once a few super-partners will be discovered at LEP/LHC, the next immediate task 
would be to reconstruct from the measured parameters, as precisely as possible, the 
structure of the soft-SUSY breaking Lagrangian. Although the relationship between 
physical parameters [mass eigenvalues, mixing angles and physical couplings] and e.g. the 
phenomenological MSSM Lagrangian is well established [^T], |22], ^3|, it would be useful to 
invert such a relationship, namely to derive Lagrangian parameters directly from physical 
parameters. This is however especially non-trivial in the neutralino sector, which involves 
the 4x4 mass matrix shown above to "de-diagonalize" . 

Here, we illustrate a relatively simple scheme for such an analytic inversion || for 
most of the Lagrangian parameters of the phenomenological MSSM, in terms of a minimal 
set of appropriately chosen physical input parameters. In the pure gaugino sector, the 
algorithm gives for a given tan /3 the values of the /x, M\ and M 2 parameters in terms of 
three arbitrary input masses, chosen indifferently among four, namely either two chargino 
and one neutralino masses or two neutralino and one chargino masses. 

Note that, in a more standard approach [i.e. from the Lagrangian to the physical 
parameters], one may obtain a similar information, e.g. by some systematic scanning of 
the whole parameter space. However, the advantage of having relatively simple [and thus 
fast] analytical expressions should be obvious, since in practice a complete scanning of the 
unconstrained MSSM parameters would be rather tedious and probably not particularly 
illuminating [for a recent detailed analysis, although in the more constrained mSUGRA 



see Ref. [ET1 for instance]. In addition, there are some subtleties in such a 



reconstruction, like the occurrence of possible discrete ambiguities which are most clearly 
seen via explicit inversion, as we shall illustrate. 

Let us now sketch our general procedure to reconstruct the gaugino sector parameters^; 
for more details we refer to Ref. ||. First, one should fix a specific choice of input /output 
parameters and a simple starting point is to assume that tan/3 is an input i.e. that it 
has been extracted from elsewhere prior to gaugino reconstruction [see e.g. section 5.2]. 
Then, we consider two basic algorithms or "scenarii", that we call Si and S 2 - 

Scenario Si', here, the input masses are assumed to be the two chargino masses m v ± 
and m x ± and one (but any) neutralino mass ±m x o. Then, the algorithm Si gives the 
parameters: /1, M 2 , Mi, plus the values of the three other neutralino masses m x o, (J 7^ i). 
More precisely, when assuming that tan (3 and the two chargino masses are given, the basic 
equations giving /j and M 2 are simply obtained by inverting the chargino mass expressions 
eqs. (pll), obtaining 



fx 2 = l -{m 2 x ± + mj± - 2M 2 , ± [(m£± + - 2M 2 V ) 2 - 4(M^ sin 2(3 ± m x ±m x ±) 2 ] 1/2 ) 

M 2 = [m 2 ± + m 2 ± - 2M 2 / - ii 2 } 1 ' 2 (3.23) 



6 Note that we restrict ourselves to real-valued parameters, but do not necessarily assume universality 
of gaugino masses; without loss of generality, one thus can choose M2 to be always positive, while the 
signs of Mi and /1 remain arbitrary, as a result of the phase re-parameterization freedom of MSSM fl52| . 
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with the sign of \i determined from 



M 2 /z = il% sin 2(3 ± m y ±m x ± (3.24) 



Note that in eq. ( |3.23| ), the ± in front of the square root simply reflects the spurious 



fi 2 ^h> M| ambiguity, while the ± inside the square root or in eq. ( p. 24 ) corresponds to 



a true ambiguity, i.e. when the expression under the root is positive (or zero) for both 
sign choice there are two possible solutions for (/i, M 2 ). The occurrence of this twofold 
ambiguity crucially depends, obviously, on the mass values w\±, m x ± an d tan/5, as will 
be illustrated. It is relatively easy to determine in which parameter domain one has either 
no solution, a unique or twofold solution. 

Concerning the neutralino mass inversion, we note first that since we restrict our- 
selves to the case where Mi, M 2 and \i are all real-valued, the neutralino mass matrix is 
symmetric and can be diagonalized through a similarity transformation, i.e. 

PMnP- 1 = A< iagonal (3.25) 

Now, rather than an analytically cumbersome inversion of the mass eigenvalue expressions, 
a simpler procedure is to start from the four basic invariants 

TrM , ^(TrM) 2 - -TrM 2 , -(TrM) 3 - -TrM TrM 2 + -TrM 3 , DetM (3.26) 
2 2 6 2 3 

under similarity transformations. These invariants contain the complete information on 
the relationship between the mass eigenvalues and the initial parameters in the neutralino 
mass matrix, but do not favor any particular set of parameters. Thus, the system may be 
solved in many different ways depending on the choice of input/output one is interested in. 
In fact, the necessary and sufficient conditions for the existence of solutions to eq. ( |3.26| ) 
take the form 



P 2 2 3 + (fi 2 + M z — M X M 2 + (M x + M 2 )S 23 - S 2 3 )B 



2.3 



+fiM 2 z (c 2 w M 1 + s 2 w M 2 ) sin 2(3 - fi 2 M x M 2 = (3.27) 



(Mi + M 2 - 5 23 )P| 3 + Gu 2 (Mx + M 2 ) + M 2 {c 2 w M 1 + s 2 w M 2 - fism2(3))P 23 
+/j(M 2 z (c 2 w M 1 + s 2 w M 2 ) sin 2(3 - fiM x M 2 )S 2?> = (3.28) 

where 5*23 = rh x o + m x o and P23 = ^x^Xjj [ w ith similar equations for all possible combi- 
nations of two neutralino masses, (m x o, m x o)]. These equations constitute our basics to 
invert the neutralino sector. For instance, in scenario Si we can extract Mi and the three 
physical masses m x o, m x o, m x o as functions of the mass m x o, from any one of eqs. ( |3.27| ) 
or ( p.28|) . Note that the mass m x o plays the role of any neutralino mass to be given as 
input, i.e. there will be a relabeling of neutralino states depending on the values of the 
other parameters. 

Scenario 5*2: here, we assume that /1 [alternatively M 2 ] is an input parameter together 
with two neutralino masses, say m x o and m x o. Then, the quadratic system eqs. ( |3.27| ) 



34 



and ( p. 28 ) gives M 2 and M\. The key point is that it is relatively simple to merge these 
two basic algorithms, S\ and S 2 , to also obtain u, M 2 , Mi consistently from m ±, mo, 
and m x o] for instance choosing an arbitrary initial guess value for m x ± [alternatively 
M 2 ], one simply has to use Si followed by S 2 , eventually iterating until a consistent, i.e. 
convergent, set of values is obtained. In most practical cases, convergence is very fast 
after 2 or 3 iterations. This peculiar decomposition, with this choice of input/output 
masses is deliberately chosen as the one giving the most algebraically tractable inversion 
in the gaugino sector. It does not imply however, a strong particular prejudice on the 
chronology of discovery of the gauginos. The most likely situation where one would 
presumably first discover two neutralinos and only one chargino, is precisely tractable 
from the combined Si + S 2 algorithm as explained above. The price to pay however, is 
that scenario Si + 5*2 [with only one chargino mass input m x ± , and without further model 
assumption] potentially gives more ambiguities than Si alone. The upshot is that up to 
four (at most) distinct solutions for (/x, Mi, M 2 ) can occur for some m x ±, m x %-> m x% m P u t 
choices once all constraints are taken into account [including in particular our necessary 
sign convention M 2 > 0]. 

Let us illustrate with some representative plots the results of the inversion in the 
gaugino sector according to the algorithms Si and S 2 . As it turns out, a number of 
general and interesting properties of the inversion can directly be seen irrespective of the 
precise values of the other parameters that have to be fixed, like tan (3 typically. 

(i) Two charginos plus one neutralino input 

We first discuss the basic algorithm Si in Fig. 2, where to exhibit as much as possible 
the dependence on the physical inputs, we fixed only one chargino mass, say m % ±, and 
varied the other mass m ±. Fig. 2 exhibits characteristics that are quite generic; there 
are three distinct zones as regards the existence, uniqueness, or possible ambiguities on 
the parameters /x, M 2 , Mi. 

a) The grey shaded region corresponds to m x ± + m x ± < 2M^ where there are no 

solutions for real \i that has to be rejected according to our basic assumptions^; if one 
takes a smaller or larger m ± value, this region around m ± will be simply displaced 
accordingly, b) In the left and right border zones are the regions of twofold ambiguities 
on /x, M 2 as indicated, c) Finally the two bands in between zones a) and b) correspond to 
the region where eqs. ( |3.23| ) give a unique solution for /x and M 2 ; note that those bands 
are narrower when tan f3 is increasing [tan /3 — 2 in Fig. 2] , irrespective of the m x ± values, 
becoming e.g. only a few GeV wide for tan/3 > 35. 

Note also that /x and M 2 are rather insensitive to tan /3, apart from the discontinuous 
change occurring for one of the solution at the border between zones b) and c). One can 
also see from Fig. 2 the relatively simple shape of zx and M 2 as function of m ±, with 
an almost constant or linear dependence apart in some narrow regions. This is easily 



7 Of course, more generally, one could be interested in complex /i solutions. However, the present 
algorithm is not entirely consistent if \x and M2 are assumed complex so that in the present context a 
complex /1 solution of eq. ( 3.23| ) has to be rejected. 
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Figure 2: The parameters u, M 2 and Mi as functions of m + for fixed m+ = 400 GeV 



m + = 50 GeV and tan/3 = 2. 
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understood since from eqs. ( |3.23| ), one obtains fi(M 2 ) ~ m x ±(m x ±) for m x ± <C m x ± or 
m ± 3> m ± . In Fig. 2 we also plot Mi for the corresponding values of u and M 2 and for 
fixed m x o [the almost constant behavior of M\ in this plot, apart from small m x ± ~ 100 
GeV, is not completely obvious but is explained in details in Ref. ||]. 

Thus, the information from the plots in Fig. 2 is that, apart from some small regions, 
for a very wide range of \m ± — m ±1 the dependence of u, M 2 [and even Mi to some 
extent] upon the latter mass difference is strongly correlated. It is straightforward to 
obtain some resulting values of the parameters /i, M 2 and Mi at the GUT scale, when 
a RG evolution of these parameters is applied after the inversion algorithm S\. The 
behavior of each parameter as a function of input masses remains essentially the same 
apart from a systematic shift due to the RG evolution. A comparison with SUSY-GUT 
model assumptions is then possible at this level. 

(ii) One chargino plus two neutralinos input 

Next we illustrate the probably more phenomenologically relevant combined scenario 
Si plus S 2 , namely where rn x ±, m x o and m x o are given as input. As expected, Fig. 3 
reflects the more involved inversion when combining algorithm Si [with unknown Tn x ±\ 
and 5*2 due to the possible occurrence of a larger number of distinct solutions for (/i, Mi, 
M 2 ) in this case. However, apart from relatively messy-looking but narrow zones where 
twofold solutions occur for this particular input choice, for a wide range of the m x 2 values 
the solution is unique at least for these input values. The shaded regions again corresponds 
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to a zone where one (at least) of the output parameters (//, Mi, M 2 ) becomes complex- 
valued (that is, for any possible solution). In Fig. 3 we only show on purpose a range 
of values such that all masses are relatively small, while for larger |m x g | the dependence 
of fi, M2, Mi upon the latter becomes simpler and almost linear, in accordance with 
the behavior in the previous figure for scenario Si alone. Also, the dependence of the 
scenario Si + S 2 upon tan/3 variations is relatively mild. In contrast, varying m x ± and/or 
m x o input values for plots similar to Fig. 3 has more drastic effects since in particular, 
the number of distinct solutions crucially depend on those inputs. More precisely, when 
varying those two input masses, some of the branches in Fig. 3 may disappear or on the 
opposite, extra branches may appear although the behavior of a given unaffected branch 
as a function of m x 2, remains essentially the same. 
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Figure 3: The parameters fi, M 2 and Mi as a function of m x o for fixed m x o = —100 GeV 
and m ± = 80 GeV for tan (3 = 2. 

Finally, let us make a remark on the other MSSM parameters inversion. In parallel 
to the reconstruction of the gaugino sector soft-breaking parameters from the physical 
masses, it is natural to attempt such a reconstruction for the remaining part of the soft- 
breaking Lagrangian. In contrast to the gaugino sector however, the de-diagonalization 
of the sfermion sector and the Higgs sector does not present much analytical difficulties, 
provided of course that one knows a sufficient number of physical masses and/or couplings. 
Here, we briefly mention that such an inversion is indeed possible and we refer for more 
details and illustrations to Ref. ||. We should only note that, in contrast with the 
gaugino sector inversion where the only difficulty was due to algebraically non-trivial 
de-diagonalization, for a realistic inversion of the Higgs parameters care should be taken 
with the correct account of one-loop corrections to the scalar potential and Higgs masses. 
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3.2 The program SUSPECTQ 
3.2.1 Introduction 

It is a well-known fact that the proliferation of Supersymmetry breaking terms in the 
unconstrained MSSM Lagrangian makes the Lagrangian-to-physical parameters [i.e. par- 
ticles masses and couplings, etc.] relationship a rather tedious task to derive in an ex- 
haustive manner. Although several systematic routines doing this work with different 
levels of refinement are available, it turned out to be highly desirable to develop our own 
tool within the GDR workshop, with the specific aim among other things, to fix a GDR- 
common set of parameter definitions and conventions once for all, and to have as much 
as possible flexibility on the input/output parameter choice. 

Such a program should necessarily include the possibility of renormalization group 
evolution of the relevant parameters, due to the all importance of the high/low energy 
scales interplay in most modern SUSY scenarii, with GUT models as a typical but not only 
example. Besides, it unavoidably poses numerous questions on the relevance of various 
possible "options" [or more precisely model assumptions] to be available, as well as the 
choice of an adequate level of approximation at different stages, such as typically, in the 
effective potential related to radiative corrections to Higgs boson masses. 

This brought us to develop the Fortran code SUSPECT which calculates the masses 
and couplings of the SUSY particles and the Higgs bosons of the MSSM. The code contains 
one source file [SUSPECT. f] and one input file [SUSPECT. in]; any choice and option is 
driven from this input file where one can change data almost at will. All results, including 
comments when useful, are to be found in the output file [SUSPECT. out] which is created 
at any run of the program. It is hoped that the code may be readily usable even with not 
much prior knowledge on the MSSM. 

The code is as yet at a rather preliminary stage [version 1.1]. As mentioned previously, 
its main purpose at present is to propose some conventions, definitions, and possible 
flexibility choices which are largely open to discussions/suggestions. In the present version, 
only two extreme "models" are readily available: the "phenomenological MSSM" [but 
with the possibility of RG evolution to arbitrary scales] and the mSUGRA model; many 
"intermediate" models can be very easily included [suggestions are particularly welcome 
here]. For instance the pMSSM with unification of the gaugino masses at the GUT scale, 
or with a common scalar mass are straightforward to implement by just setting constraints 
"by hand" in the input file. The output for the MSSM mass spectrum has been compared 
with several other existing similar codes [see the "tools for SUSY" working group report 



33[ for more details on those comparisons], and is in reasonably good agreement. Up 
to a few percent differences in the output is to be sometimes expected, however, due to 
different approximations used in different codes. 



8 During the last GDR-SUSY meeting in Montpellier, some people complained about the former name 
of the fortran code: MSSMSPEC seemed to be difficult to pronounce by some of our (presumably non 
Slavic) colleagues, probably due to a local cluster of consonants. We propose a change of name to 
SUSYSPECT, or to make short SUSPECT [since every code, is, a priori!]. 
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The rest of this subsection is organized as follows: we first describe the most impor- 
tant subroutines in the case of the "phenomenological" MSSM in particular those which 
compute the chargino/neutralino, sfermion and Higgs boson masses and couplings. We 
will then discuss the case of the constrained MSSM with the various choices of approxima- 
tions and refinements, paying a special attention on the RG evolution and the radiative 
electroweak symmetry breaking. Next, a description of the input and output files to be 
driven by the user is given. In doing so we shall try to be as close as possible to the 
definitions and conventions given in the previous subsections. We finally collect a list of 
various available [or not yet] options and/or model assumptions, with a clear mention of 
eventual limitations of the present version of the code. 

3.2.2 The "phenomenological" MSSM 

The purely phenomenological MSSM, that is, with the set of input parameters as defined 
in section 2.1, and without any RG evolution, is implemented by enforcing in the SUS- 
PECT. in input file the choice ICHOICE(1)=0. Another alternative of interest is to have 
the full set of parameters as defined by the phenomenological model, but with the further 
possibility to evolve those parameters to -or from- an arbitrary scale prior to calculating 
the physical spectrum. This latter case is readily implemented and enforced by the choice 
ICHOICE(l)=l. 

There are three main subroutines evaluating all possible SUSY physical masses, respec- 
tively chargino/neutralinos, Higgs boson and sfermion masses, for given input Lagrangian 
parameters: 

CALL GAUGINO (MU , Ml , M2 , M3 , BETA , ALPHA , GMC , GMN , XMN) 
CALL SUSYCP (TGBETA) 

CALL SFERMION (MSG) , MTR , MBR , MSL , MTAUR , AL , AU , AD , MU , 
GMST , GMSB , GMSL , GMSU , GMSD , GMSE , GMSN) 

For the couplings, only a subset is implemented at present: the couplings of the MSSM 
Higgs bosons to standard and SUSY particles [the angles a and (3 = atan(tan/3) of the 
MSSM Higgs sector are called by a common]. Note that the couplings are not readily 
default output in the SUSPECT1.1 version, but are straightforwardly available, from var- 
ious common blocks, as specified below. The remaining couplings [sfermion/ino couplings 
to gauge bosons and fermion-sfermion-ino] will be included soon. 

In the subroutine GAUGINO, the input parameters are //, the gaugino masses Mi, M 2 , M 3 
[which can be made related], the mixing angles f3 and a in the Higgs sector [the former 
is an input while the latter is calculated in the subroutine SUSYCP as explained below]. 
Signs and other relevant conventions are consistent with the ones discussed in section 3.1. 
The output is: 

GMC(l), GMC(2): ordered values of the chargino masses 
GMN(1)-GMN(4): ordered absolute values of the neutralino masses 
XMN(1)-XMN(4): neutralino masses including the eventual negative signs. 
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The neutralino mass matrix is diagonalized analytically. The negative signs in XMN(l)- 
XMN(4) are relevant for instance for the couplings to Higgs bosons; if needed, these 
couplings are available via a common defined in the subroutine: 

AC1(2,2),AC2(2,2),AC3(2,2): chargino-/i, H, A couplings 
AN1(4,4),AN2(4,4),AN3(4,4): neutralino-/t, H, A couplings 
ACNL(2,4),ACNR(2,4): chargino-neutralino-if ± couplings. 



The Higgs boson spectrum and couplings are calculated by the subroutine SUSYCP. 
It uses as two important input parameters: tan (3 and the mass of the pseudoscalar Higgs 
boson Ma which may be either given as direct input or calculated from the two soft- 
SUSY breaking scalar masses Mh v and M# d and \i as in eq. (|2.28|) . In addition, SUSYCP 
uses as input the third generation soft scalar masses and trilinear couplings as well as fi 
which enter in the radiative corrections. In the case where the pole Higgs boson masses 
are calculated it needs the parameter M 2 , as well as the masses of stops and sbottoms 
and their couplings to Higgs bosons [which are calculated similarly as in the subroutine 
SFERMION]. SUSYCP calls the subroutine SUBH [[54]], which calculates the renormal- 
ization group improved values of the Higgs boson masses and the couplings A in the scalar 
potential^. The SUSYCP output gives the SM and A, h, H, H ± Higgs boson masses and 
the running A mass [the SM Higgs boson mass is in fact the input M&\ 

AMSM,AMA,AML,AMH,AMHC,AMAR 

Depending on the flag ICHOICE(IO) one has either the running (=0) or pole (=1) Higgs 
boson masses masses [in the former case AMAR=AMA]. The output also provides via 
common blocks the self-couplings Ai_7 in the scalar potential, and the Higgs couplings 
to fermions, gauge bosons and the self-couplings [the notation being obvious, with T,B,L 
standing for t, b, r and L,H,A for h, H, A; GLPM=sin(/5 — a) and GHPM=cos(/3 — a)} as 
well as the running angles /3 and a: 

COMMON/HMASS/AMSM,AMA,AML,AMH,AMCH,AMAR 
COMMON/HSELF/LAl,LA2,LA3,LA4,LA5,LA6,LA7 

COMMON/COUP/GAT, GAB, GLT,GLB,GHT,GHB,GZAH,GZAL,GHHH,GLLL, 
GHLL,GLHH,GHAA,GLAA,GLVV,GHVV,GLPM,GHPM,B,A 



In the subroutine SFERMION, the input parameters are the soft-SUSY breaking mass 
parameters for left-handed and right-handed squarks and sleptons, the trilinear couplings 
At, At, and A T for the third generation and the parameter /i. The output are 2-vectors: 

GMST,GMSB,GMSL,GMSU,GMSD,GMSE,GMSN 

for the masses of the two t, b, f eigenstates, their first/second generation partners as well 

9 Note that compared to the original version [M, in the one used in the program some bugs [e.g. 
in the calculation of the b contributions to the running Higgs boson masses] have been corrected and 
some improvements [such as using the standard parameters given in the input file and calculating the 
running quark masses and a s using the subroutines RUNM and ALPHAS, as well as calculating the 
Higgs couplings to squarks as in SFERMION, see below] have been made. Furthermore, the function F0 
which calculates the loop functions is given explicitly. 



40 



as the mass of the P's. The D-terms are calculated by the subroutine and full mixing 
in the third generation is implemented. Again, the couplings to Higgs bosons are not 
explicitly written but are available if desired from a common defined in SFERMION: 

GLEE(2,2),GLTT(2,2),GLBB(2,2): h couplings to f,t,b 
GHEE(2,2),GHTT(2,2),GHBB(2,2): H couplings to f,t,b 
GAEE,GATT,GABB: A couplings to fif 2 , tjz, hh 
GCEN(2,2),GCTB(2,2): H± couplings to fu, tb. 

The conventions for the signs of Af,fi as well as the mixing terms are as introduced 
previously. [The sfermion masses are running masses, since they are obtained from the 
Lagrangian parameters of the scalar sector; note indeed that for stops and sbottoms the 
running quark partner masses are implemented in the mass matrices]. 

Of course, the three subroutines need in addition some standard input parameters. 
Apart from "electroweak" coupling constant inputs [such as the Fermi constant or the 
QED coupling defined at zero-momentum] and some particles masses [such as the W,Z 
gauge boson and fermion masses] which are precisely known and which are included in 
the source code, there are other less precisely known parameters which are very important 
[especially in the constrained MSSM, since e.g. the GUT and SUSY thresholds are quite 
sensitive to the values of some of them]. These constants [defined at the low scale, where 
they are experimentally measured, currently Mz] are: 

ALFINV: the inverse QED constant l/a(M z ) in the MS scheme (=127.9) 
SW2: sin 2 6 W (M Z ) in the MS scheme (=0.2315 for m t = 175 GeV); 
ALPHAS: the value of a s {M z ) at the M z scale (=0.119 ); 
MT,MB,MC: the t, b, c pole masses (=175, 4.7, 1.42 GeV) 

The reference values for those parameters as indicated in parenthesis above correspond to 
central values quoted by the most recent data |55| ; they may accordingly be varied within 
the known error bounds. 

3.2.3 Constrained MSSM 

The case of mSUGRA is dealt with by the flag ICHOICE(1) = 10 in SUSPECT.in. In the 
mSUGRA choice, all the SUSY parameters are determined in terms of four arbitrary and 
one discrete parameters which are given at a "high" scale which, if unification of the gauge 
couplings is imposed, corresponds to the GUT scale M\j. These parameters, as discussed 
in section 2.1, are: 

TGBETA(MZ) , mO(GUT), mHALF(GUT) , AO(GUT), sign(MU) 

Other parameters also enter the game as "starting guess" values, and must be specified 
in some cases: 

MU, B, SUSYM, EHIGH, ELOW 

Typically, the actual values of fi and B are determined consistently from EWSB, but 
the latter almost unavoidably needs some iteration procedure, and thus some starting 
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"seed" fi, B values are necessary. If the unification of the gauge couplings scenario is not 
chosen two additional inputs must be given: SUSYM corresponding to an initial guess 
of the SUSY threshold [a single universal one to simplify], and EHIGH corresponding to 
the high energy end of the RG evolution [if gauge unification is enforced, both SUSYM 
and EHIGH can be safely set to arbitrary values, as more adequate values will be then 
calculated]. Finally, one needs to specify ELOW which corresponds to the lower-energy 
end of the RG evolution equations; this is always necessaryp ] and important as it fixes 
the final scale at which all physical masses will be computed. Moreover, ELOW should 
be below or at most equal to the EWSB scale [see the discussion in section 2.2]. Thus 
ELOW should be roughly of the order of the electroweak scale, with some flexibility, 
say Mz < ELOW < 1000 GeV. In the mSUGRA case all other parameters are then 
determined by GUT universality and RG evolution. 

In addition, independently of the model choice driven by ICHOICE(l), there are 
several possibilities for the RG evolution equations and for electroweak symmetry breaking 
[which are however not always all relevant, depending on the choice of models]: 

• ICHOICE(2): elaboration level of the RG evolution equations. 

= 11 (21): one (two)-loop evolution with simple [single SUSY-threshold scale] thresh- 
old effects. 

=12 (22): one (two)-loop evolution, with more realistic thresholds [i.e. step func- 
tions for each particle specy thresholds]. 

• ICHOICE(3): Gauge coupling and Yukawa coupling unification. 
= 0: no unification 

= 1: only gauge coupling unification [then approximate GUT and simple SUSY- 
threshold scales are automatically determined]. 
= 2: gauge and bottom-tau Yukawa coupling unification 

• ICHOICE(4): Accuracy of RG evolution. 
= 0: fast but less accurate RG evolution. 
= 1: accurate evolution but rather slow. 

• ICHOICE(5): for electroweak symmetry breaking. 

=0: no radiative EWSB [electroweak breaking is of course there and determines 
consistent values of B and /i at the EWSB scale; the difference between "radiative" 
and "no radiative" EWSB was explained in section 2.2]. 

=1: radiative EWSB is implemented, i.e. B and // are fixed at the EW scale and by 
other parameters to be consistent with RG evolution from a given (generally EHIGH 
~ Mu) scale. [NB: ICHOICE(5) = l is automatically enforced for the mSUGRA case 
ICHOICE(1)=10] 

10 For all mass and coupling calculations performed by the routines GAUGINO, SUSYCP and 
SFERMION described above, a scale is implicitly involved since the input Lagrangian parameters are 
implicitly defined at some scale, the latter being eventually the one resulting from an RG evolution; see 
the discussion below. 
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• ICHOICE(6)-ICHOICE(9) are reserved for later use. 

• ICHOICE(IO): =0 (1) computes running (pole) Higgs boson masses^]. 

Note that we anticipated several possible steps between the pMSSM and mSUGRA 
and therefore left the possibilities ICHOICE(6-9) reserved for later use. Besides the 
physical spectrum calculation for input Lagrangian parameters as described above, the 
two other main tasks performed by the code are the RG evolution of the parameters, and 
the consistent implementation of EWSB, which we discuss is some more details now. 

(i) Renormalization Group evolution: 

ODEINT is a subroutine returning a set of masses and coupling parameters at a speci- 
fied scale Qout> when given at an initial scale Q in . It is based on a Runge-Kutta numerical 
algorithm solving differential equations by Numerical Recipes f56fl . RG evolution is triv- 
ially reversible, so that the choice of Q in and Q out is in principle quite flexible; it is also 
possible in particular to use this evolution back and forth if needed. 

CALL ODEINT (y , n , xl , x2 , eps , hi , 1 . d-5 , nok , nbad , derivl , rkqc) 

The input are the values y(n) of all evolving couplings and masses relevant to the phe- 
nomenological MSSM, given at the scale Q in = e Xl , the latter being also given as input; 
see section 2.3 for the precise assignment of y(n) components in terms of MSSM couplings 
and masses. [Quite obviously, when specific constraints on the y(n) are implemented at 
some scale, like in the mSUGRA model typically, it simply implies specific initial condi- 
tions for those. But the evolution of each mass or coupling parameters is driven by its 
own specific beta function and we therefore should keep the latter as model-independent 
as possible]. The output are the values of y(n) at the scale Q out = e X2 to be specified. 

According to the previous discussion ODEINT also needs the beta functions: /3(y)(= 
d(y)/dha.Q), for all the relevant y(n), which are provided by the subroutines: 

subroutine DERIVl(x,y,dydx) 
subroutine DERIV2(x,y,dydx) 

for the one and two-loop RG beta functions respectively, including threshold effects. 

(ii) Electroweak Symmetry Breaking: 

There is a set of subroutines [in fact essentially one] which evaluates the one-loop 
contribution to the effective potential, whose detailed shape study settles the occurrence 
[or not] of the SU(2) x U(l) spontaneous symmetry breaking, irrespective of particular 
model building assumptions. What is however, a model assumption chosen as option 

11 Physical masses are, strictly speaking, defined as the pole masses: M polc = m(M polc )+ radiative 
corrections, where m(Q) is the corresponding running mass, whose evolution with energy scale Q is 
entirely dictated by the RG; the radiative correction parts are all non-RG finite loop corrections. However, 
apart from the exception of the Higgs boson and gluino masses, for most other MSSM particles the non- 
RG corrections are assumed to be negligible with respect to the present level of approximation and are 
not included in SUSPECT. 
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flag, is whether the EWSB is radiative or not; namely whether the relevant parameters of 
the Higgs potential at the EWSB scale are consistent or not with their respective values 
chosen at a higher scale, given their specific RG evolution properties [see again section 
2.2 for a more complete discussion on EWSB]. 

call VL00P2 (Q2 , mt2 , mst22 , mst 12 , mb2 , msb22 , msbl2 , 
mtau2 , mstau22 , mstaul2 , -rmu , dVdvd2 , dVdvu2) 

The inputs are: 

Q2: the scale Q 2 (> ELOW) at which EWSB is supposed to happen; 
mt2, mst22, mstl2: the t,ti,t 2 masses; 
mb2, msb22, msbl2: the b,bi,b 2 masses; 
mtau2, mstau22, mstaul2: the r,fi,f 2 
—RMU: the parameter fi 

The output are dVdvd2, dVdvu2, which are [up to some appropriate overall constants] 
the derivatives d{y e s)/d{y J), d(V c f{) / d(v 2 ) respectively, the basic quantities entering Higgs 
mass corrections and EWSB consistency conditions. As already mentioned, our present 
algorithm calculates values of the parameters B and \i consistent with EWSB by iterating 
the two expressions of fi 2 and B, eqs. ( |2.28| ) of section 2.2. The iteration is non trivial 



because in eqs. (|2 .28 



M 2 Hd = M 2 Hd + d(V cS )/d(v 2 d ) and M 2 Hu = M 2 Hu + d(V cS )/d(v 2 u ) (3.29) 

involve other MSSM particle masses and couplings, which themselves depend on /z [in most 
practical cases, however, the algorithm converges very fast, after 2 or 3 iterations]. Besides 
eqs. Q2.28D , there are other constraints to fulfill in order to assert that one reached a local 
minimum of the potential, as discussed in section 2.2. This is simply checked numerically 
in SUSPECT. In contrast, there is no consistency checks about the possible occurrence 
of charge and color breaking [CCB] minima in the present version 1.1. Note, however, 
that a more refined different algorithm, partly based on analytical studies of the complete 
potential at its minimum M, will be implemented later. 

3.2.4 Example of input/output files 

The following example of SUSPECT. in input file, corresponds to the mSUGRA so-called 
"SNOWMASS" point 1: 

1) choice of different options/models, etc flags: 

ICHOICE(l) = 10 

ICH0ICE(2) = 21 

ICH0ICE(3) = 12 

ICH0ICE(4) = 

ICH0ICE(5) = 1 

ICH0ICE(6) = 

ICH0ICE(7) = 
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ICH0ICE(8) = 
ICH0ICEC9) = 
ICH0ICE(10)= 1 

2) Initialize SM parameters: 
ALFINV =127. 9d0 

SW2 =.2315d0 

ALPHAS =.12d0 

MT =175. dO 

MB = 4.7d0 

MC =1.42d0 

3) Initialize SUSY parameters 
SUSYM = 200. dO 
EHIGH(GUT=2.dl6 

ELOW = 433. dO 

(here SUGRA case with universality) : 

mO = 400. dO 

ml/2 = 200. dO 

AO = O.dO 

TGBETA = 2.d0 

SGNMUO = -l.dO 

4) pMSSM case (at HIGH scale!) 



MHU2 




.5D4 


MHD2 




5.d5 


Ml 




200. DO 


M2 




200. DO 


M3 




200. DO 


MSL 




1.D2 


MTAUR 




1.D2 


MSQ 




l.d2 


MTR 




l.d2 


MBR 




1.D2 


MEL 




1.D2 


MER 




1.D2 


MUQ 




1.D2 


MUR 




1.D2 


MDR 




1.D2 


AL 




200. DO 


AU 




200. DO 


AD 




200. DO 



The corresponding SUSPECT. out output file is shown below [the comments are self- 
explanatory] : 
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HIGH (GUT) , LOW (Final) and SUSY-threshold SCALES: 

0.274650E+17 433.000 175.000 
EVOLVE FIRST TIME from MZ to HIGH (GUT) scale 

gauge cpl~2 1,2,3 at HIGH (GUT) scale 0.529596 0.526736 0.517345 
RUNNING mtau,mb,mtop at HIGH scale 1.58621 1.18050 76.0442 
EVOLVE DOWN to Low-energy 



RESULTING values at Low-energy scale: 
tan (beta) at Low scale: 1.93532 

mtau,mb,mtop at Low scale: 1.76582 2.69798 151.422 

Atau,Ab,Atop at Low scale: -135.228 -633.479 -368.151 

m(phi_u)~2, m(phi_d)~2 -137416. 179919. 
m_tauR, m_L, m_bR, m_tR, m_Q: 

407.225 424.685 607.580 399.425 530.837 

rel. % err. in B, mu: 0.280874E-03 0.338679E-03 
(indicate departure from consistent radiative EW breaking) 
loop-level EW breaking + stability tests: 
ml~2*m2"2-B"2*MU~2, ml~2+m2~2 +/-2*B*MU: 
-0.767807E+09 994235. 100948. 
FINAL SOFT-BREAKING PARAMETERS: 
mu(Q_EWB), B(Q_EWB): -504.916 -442.295 
gauginos 1,2,3 masses: 83.3875 165.443 524.285 



FINAL PHYSICAL PARAMETER RESULTS: 
tgbeta M_A M_h M_H M_H+ 

2.000 735.902 79.555 742.020 743.895 

chi_+~l chi_+~2 chi_0~l chi_0~2 chi_0~3 chi_0~4 
172.336 515.130 85.597 172.466 506.938 516.171 

stop_l stop_2 sup_l sup_2 
429.182 557.382 621.448 608.878 

sbottom_l sbottom_2 sdown_l sdown_2 
532.761 607.879 624.448 608.114 

stau_l stau_2 selec_l selec_2 sneut 
408.361 426.429 426.354 408.756 421.948 
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3.2.5 Discussions and outlook 

In summary, the main core of the code SUSPECT gives the physical masses and couplings 
of the SUSY particles and the MSSM Higgs bosons as functions of input Lagrangian pa- 
rameters, taking into account renormalization group evolution if needed, and consistent 
electroweak symmetry breaking. In addition, there are a number of options, which might 
be useful for many practical purposes, that we briefly describe now. 

INVERTOR: 

This is a subroutine [based on the discussion in section 3.1.2], which will be included soon 
in the next version, as an option in SUSPECT and whose purpose is to determine the 
inverted spectrum relationship [i.e. recovering the Lagrangian parameter values directly 
from physical masses and/or couplings]. The algorithm in its present form essentially 
deals with the non-trivial inversion in the gaugino parameter sector, where the input can 
be either two charginos and one neutralino, or two neutralinos and one chargino physical 
masses. The output are the Lagrangian gaugino parameters, /i, Mi, M 2 [the M 3 to M g i uirio 
relation obviously does not need inversion]. 

RGEXACT: 

This is a subroutine also to be included in the next version as an option in SUSPECT 
which implements an exact RG evolution solver [limited to one-loop approximation]. We 
refer to section 2.3 for a more detailed discussion of the procedure. Let us simply men- 
tion here, that the exact solutions of the relevant SUSY RG equations which have been 
recently derived for arbitrary values of tan j3, should not only be useful to improve the 
general RG evolution algorithm, but more importantly, should provide a better control 
on some non-trivial issues of the RG evolution, such as the possible occurrence of Landau 
poles in the Yukawa couplings typically. 

Finally we emphasize that the code is largely in a developing stage, and many other 
extensions of the present algorithms and options are foreseen, as well as interfaces with 
other numerical codes for SUSYQ. Some of the present limitations and approximations 
of the version SUSPECT1.1 are: 

(i) Model choice: 

Only two extreme "models", phenomenological MSSM on one side and mSUGRA with 
all universal terms on the other side, are treated at present. Many "intermediate" models 
can be enforced by hand by choosing appropriately the input parameters and will be very 
easily included later. 



12 For instance, the three subroutines GAUGINO, SFERMION and SUSYCP are also used in the 
program HDECAY which calculates the branching ratios of the MSSM Higgs bosons and that SUSYCP 
is used in HPROD [see next section] for the Higgs production cross sections in e + e~ collisions: an interface 
is therefore easily possible and will be done soon. 
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(ii) Scale choice: 

The input scale at which to define parameters is not at present completely flexible: in 
the phenomenological MSSM, all SUSY parameters should be entered at the same scale 
[the latter scale is however arbitrary, provided it is chosen between 1 TeV and say 10 18 
GeV] . An important exception is tan f3, which should be given at the Z-boson scale and 
of course the SM parameters such as a s , etc.. [see the example of input file]. 

(Hi) Approximations in specific calculations: 

There are approximations in the RG evolution, for instance an option on one-loop versus 
two-loop evolution, where in the latter only the gauge and Yukawa couplings truly in- 
volve two-loop contributions. Also, there are approximations in the contributions to the 
one-loop corrections to the scalar potential [relevant for the radiative EWSB consistency 
conditions]: that is, only the third generation fermion and sfermion contributions are 
included. This implies, strictly speaking, some inconsistency with the RG evolution since 
in the latter much more contributions are taken into account. The inconsistency however, 
should not be very significant numerically. A more complete and consistent treatment of 
the one-loop scalar effective potential will be available in the next version of the program. 

(iv) Threshold effects: 

These are not yet fully realistic at the present moment: only a single, average SUSY 
threshold scale can be chosen. Inclusion of more realistic thresholds [i.e. for each particle 
species] will be soon available. 

(v) Electroweak symmetry breaking: 

As explained previously, the EWSB conditions are consistently implemented numer- 
ically by iteration on the parameters ji and B, and the occurrence of a local minimum 
is checked numerically. There is however, no check about the possible occurrence of in- 
consistent charge and color breaking minima in the present version. These should be 
implemented in the next version, together with a complete and more analytical treatment 
of the EWSB conditions. 

(vi) Known bugs: 

Finally to be fair one cannot omit mentioning that one serious bug that has been iden- 
tified in this 1.1 version, is not yet completely understood. The observation is that, for 
some choices of the low scale input parameters, and in particular for a relatively low value 
of tan/5, < 1.3-1.4, there is an uncontrollable growth of the Yukawa/gauge couplings, 
when running the RGE in the upwards direction. This seemingly "Landau pole" effects 
is not yet under control and we hope to solve this problem for the next version [the exact 
RG solver RG EXACT will certainly help us to resolve the issue]. 

The next version SUSPECT1.2 should therefore be hopefully more complete. 
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4. Higgs Boson Production and Decay 



One of the main motivations of supersymmetric theories, as discussed in section 2, is the 
fact that they provide an elegant way to break the electroweak symmetry and to stabilize 
the huge hierarchy between the GUT and Fermi scales |21, 22 1. The probing of the Higgs 
sector of the MSSM [^] is thus of utmost importance. The search for the CP-even Higgs 
bosons h and H, the pseudoscalar Higgs boson A and the charged Higgs particles H ± 



of the MSSM is therefore one of the main entries in the LEP2 
and future e + e~ colliders agendas. 



571, Tevatron PJ, LHC 



39. 



In the theoretically well motivated models, such as the mSUGRA scenario discussed 
previously, the MSSM Higgs sector is in the so called decoupling regime |62|] for most of 
the SUSY parameter space allowed by present data constraints ||55|| : the heavy CP-even, 
the CP-odd and the charged Higgs bosons are rather heavy and almost degenerate in 
mass, while the lightest neutral CP-even Higgs particle reaches its maximal allowed mass 
value M h < 80-130 GeV @ 



depending on the SUSY parameters. In this scenario, 
the h boson has almost the same properties as the SM Higgs boson and would be the sole 
Higgs particle accessible at the next generation of colliders. 

In this section, we discuss the production of the MSSM Higgs bosons at the next 
generation of colliders, mainly at the LHC and a future e + e~ linear collider as well as 
their decays modes, in particular the SUSY channels. The new material presented in this 
section is based on Refs. R, H, 0, [TO, O. 



4.1 MSSM Higgs boson production at the LHC 
4.1.1 Physical set— up 



At the LHC, the most promising channel [59, ^Kj for detecting the lightest h boson is the 
rare decay into two photons, h — > 77, with the Higgs particle dominantly produced via 
the top quark loop mediated gluon-gluon fusion mechanism |63|, EM 



99 



h 



(4.1) 



In the decoupling regime, the two LHC collaborations expect to detect the narrow 77 
peak in the entire Higgs mass range, 80 < Mh < 130 GeV, with an integrated luminosity 
J C ~ 300 fb _1 corresponding to three years of LHC running fl59"| ; see Fig. 4. 

Two other channels can be used to detect the h particle in this mass range: the 



production in association with a W boson |65|] or in association with top quark pairs p6 



PP 



hW and pp — ► tth 



(4.2) 



with the h boson decaying into 2 photons and the t quarks into b quarks and W bosons [for 
the latter process, the Higgs boson detection with h — > bb final states looks also promising; 
see Ref . |67j for instance] . Although the cross sections are smaller compared to the gg — > h 
case, the background cross sections are also small if one requires a lepton from the decaying 
W bosons as an additional tag, leading to a significant signal. Furthermore, the cross 
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section a(pp — > ith) is directly proportional to the top-Higgs Yukawa coupling, the largest 
electroweak coupling in the SM; this process would therefore allow the measurement of this 
parameter, and the experimental test of a fundamental prediction of the Higgs mechanism: 
the Higgs couplings to fermions and gauge bosons are proportional to the particle masses. 

The additional vector boson fusion mechanisms |68| . 



pp -> W*W*/Z*Z* -> qqh (4.3) 

is less interesting to detect the h boson than in the SM, but can be useful to test the h 
particle properties [for instance to measure the hWW and hZZ couplings]. 

The heavy CP-even H and CP-odd A bosons can be searched for at the LHC through 
their decays modes into t + t~ pairs with the Higgs bosons produced in the gg fusion 
mechanism or in association with bb pairs: 

gg — > H/A and gg, qq^bb + H/A (4.4) 

This needs large values tan/3 > 5 for the Higgs boson masses M H ^ ^ 300 GeV, to enhance 
the cross sections of the above processes and the t + t~ decay branching ratios, which in 
the case where only standard modes are allowed reach the asymptotic value of ~ 10% 
||69|| . The decays of the H and A bosons into muon pairs, H/A — > n + fi~, give a rather 
clean signal and can be used despite of the very small branching ratios, ~ 4.10 -4 in the 
asymptotic limit. For lower values tan/3 < 3 [most of which will be covered by the upgrade 
of LEP to y/s = 200 GeV] and not too large Mh,a values, the decays H/A — > tt can be 
also used. The decays H — > hh — > 6677 and A — > Zh — > l + l~bb can also be used in a tiny 
area of the MSSM parameter space; see Fig. 4. 

At the LHC, the only way to detect the charged Higgs boson is when M#± < m t — mj 
and the H ± particles can be then produced in top quark decays, 

pp -> tt with t -> H + b (4.5) 

The usual signature is to look at a breakdown of lepton universality by selecting the decay 
H + — > t + u t , which together with the decay H + — > cs for small tan/5 values, is considered 
as the main decay mode [see however the discussion in section 4.2]; Fig. 4. 

In the next subsections we will discuss the effects of supersymmetric particles in the 
production of the lightest h boson at the LHC: first the contributions of light stops in 
the production of the h boson in the gg — > h mechanism which can significantly alter 
the expected production rate and then the associated production of the h boson with 
i-squark pairs for which the cross section can be rather large, exceeding the one for the 
SM-like process pp tth 0. 

One of the most important ingredients of these discussions is that stops can alter 
significantly the phenomenology of the MSSM Higgs bosonsQ. The reason is two-fold: 
(i) as discussed in section 3.1, the current eigenstates, ti and in, mix to give the mass 

13 This might also be the case of the sbottoms for large values of tan/3 and the parameters \i and A; 
however this will not be discussed here and we will assume that the mixing is zero in this sector. 
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eigenstates t\ and i%\ the mixing angle % is proportional to m t A t , with A = A t — \xj tan /3, 
and can be very large, leading to a scalar top squark t\ much lighter than the t-quark and 
all other scalar quarks; (ii) the couplings of the top squarks to the neutral Higgs bosons h 
given in section 3.1, involve components which are proportional to A t and for large values 
of this parameter, the coupling can be strongly enhanced. 

The strong h couplings^ to stops would result in a stop contribution to the hgg 
and hjj vertices that is comparable or even larger than the top contribution, altering 
significantly the rate for the process gg —>■ h —>■ 77. It also leads to a possibly substantial 
cross section for the production of the h boson in association with stops, pp — > tth. 




m A (GeV) 



Figure 4: MSSM parameter space with the contours of the Higgs bosons visible at the 
LHC with a luminosity / Cdt = 300 fb" 1 and combining the experimental data of ATLAS 
and CMS; from Ref. ]7|. 



14 The couplings of the neutral heavy Higgs bosons H, A to t-squarks and of the charged Higgs bosons 
to £6-squarks involve also large components as discussed in section 3.1. For small tan/3 values, the 
decays of these heavy Higgs bosons to stop squarks if kinematically allowed can be therefore enhanced 
by these strong couplings, as will be discussed in the next section. 
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4.1.2 Higgs production in the gluon fusion mechanism 



In the SM, the Higgs-gluon-gluon vertex is mediated by heavy [mainly top and to a lesser 
extent bottom] quark loops, while the rare decay into two photons is mediated by W- 
boson and heavy fermion loops, with the jy-boson contribution being largely dominating. 
In the MSSM however, additional contributions are provided by SUSY particles: q loops 
in the case of the hgg vertex, and H ^ , / and \ loops in the case of the h — > 77 decay. 
In the latter case JnJ, the contributions of bosons, sleptons and the scalar partners 
of the light quarks are in general small given the experimental bounds on the masses of 
these particles [[55]. Only the contributions of relatively light t squarks [and to a lesser 
extent b for large tan (3 values and xf f° r masses close to 100 GeV, which could contribute 
at the 10% level] can alter significantly the loop induced hgg and ^77 vertices. 

The expressions of the partial width for the decay h — > gg, can be found in Ref. [|64}] ; 
the cross section cr(gg — > h) is directly proportional to the decay width Y(h — > gg). The 
latter cross section is affected by large QCD radiative corrections |6^| ; however they are 
practically the same for quark and squark loops, and if only deviations compared to the 
standard case are considered, they drop out in the ratios. The partial width for the 
decay h — > 77 can be found e.g. in Ref. |7l| ; the QCD corrections are small and can be 
neglected. The 77 and gg decay widths of the h boson are evaluated numerically with 



the help of an adapted version of the program HDECAY [72]. 

Figs. 5 and 6 show, as a function of A t for tan/3 = 2.5, the deviations from their SM 
values of the partial decay widths of the h boson into two photons and two gluons as well 
as their product which gives the cross section times branching ratio o(gg — > h — > 77). 
The quantities R are defined as the partial widths including the SUSY loop contributions 
[all charged SUSY particles for h — * 77 and squark loops for h — > gg] normalized to 
the partial decay widths without the SUSY contributions, which in the decoupling limit 
correspond to the SM contributions: R = T MSSM /T SM . 




52 




In Fig. 5, the stop mass is set to = 200 GeV. For small values of A t there is no 
mixing in the stop sector and the dominant component of the hit couplings, eq. ( 3.18|) , 
is the one proportional to m^/M| [here, both t\ and ti contribute since their masses and 
couplings to h are almost the same]. The sign of this component, compared to the hti 
coupling, is such that the t and t contributions interfere constructively in the hgg and 
hjj amplitudes. This leads to an enhancement of the h — > gg decay width up to 60% in 
the MSSM. However, the h — > 77 decay width is dominated by the W amplitude which 
interferes destructively with the t and t amplitudes, and the t contributions reduce the 
h — > 77 decay width by an amount up to —20%. The product R(gg — > 77) in the MSSM 
is then enhanced by a factor ~ 1.2 in this case. 

With increasing A t , the two components of g^ti [which have opposite sign since 
sin 26$ oc m t A t in eq. ( p. 18 )] interfere destructively and partly cancel each other, resulting 
in a rather small stop contribution. For larger values of A t , the second component of 
9hhti becomes the most important one, and the t\ loop contribution [£2 is too heavy 
to contribute] interferes destructively with the one of the t-quark. This leads to an 
enhancement of R(/i — > 77) and a reduction of K(gg — > h). However, the reduction 
of the latter is much stronger than the enhancement of the former [recall that the W 
contribution in the h — > 77 decay is much larger than the t contribution] and the product 
K(gg — > 77) decreases with increasing A t . For A t values of about 1.5 TeV, the signal for 
gg — ► h —>■ 77 in the MSSM is smaller by a factor of ~ 5 compared to the SM caseQ. 

Fig. 6 shows the deviation K(gg — > 77) with the same parameters as in Fig. 5 but 
with different t\ masses, mf = 165, 400 and 600 GeV. For larger masses, the top squark 
contribution oc will be smaller than in the previous case. In the no-mixing case, 

the enhancement (reduction) of the hgg(h'yy) amplitude is only of the order of 10% for 

15 Note that despite of the large (t, b) mass splitting generated by large At values, the contributions of 
the isodoublet to high-precision observables stay below the experimentally acceptable level p3| . 
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~ 400 GeV, and leads to an almost constant cross section times branching ratio 
for the gg — > h — > 77 process compared to the SM case. Again the stop contribution 
vanishes for some intermediate value of A t , and then increases again in absolute value for 
larger A t . However, for m^ ~ 400 GeV, the effect is less striking compared to the case 
of m~ tx = 200 GeV, since here a(gg — > h) x BR(/i — > 77) drops by less than a factor of 2, 
even for extreme values of A t ~ 2.5 TeV. As expected, the effect of the top squark loops 
will become less important if the t\ mass is increased further to 600 GeV for instance. In 
contrast, if the stop mass is reduced to mi ~ 165 GeV, the drop in K(gg — > 77) will be 
even more important: for A t ~ 1.5 TeV, the gg — > 77 cross section times branching ratio 
including stop loops is an order of magnitude smaller than in the SM. For A t ~ 1.3 TeV, 
the t amplitude almost cancels completely the t/b quark amplitudes; the non-zero value 
of R(gg — > 77) is then due to the imaginary part of the 6-quark contribution. 

One should recall that M h varies with A t , and no constraint on M h has been set in 
Figs. 5-6. Requiring Mh ?t 90 GeV, the lower range A t < 350 GeV and the upper ranges 
A t > 1.5(2.3) TeV for m^ = 200(400) GeV for instance, are ruled out. This means that 
the scenario where K(gg — > 77) > 1, which occurs only for A t < 300 GeV for = 200 
GeV is ruled out for M h > 90 GeV. Therefore, the rate for the gg — > 77 process in the 
MSSM will always be smaller than in the SM case, making more delicate the search for 
the h boson at the LHC with this process. 

For large values of tan/5, tan/3 3> 1, the mixing in the sbottom sector can also be very 
large, leading to possibly rather small, and a large g h ii coupling which can also 
generate large bi loop contributions to the hgg and /177 vertices. Indeed, for tan/3 ~ 50 
and m~ bi = 200 GeV, the deviations of the K(h — > gg) and thus K(gg — > h — > 77) 
observables from unity are substantial for large values of For instance, for ~ 1 
TeV the gg — > 77 cross section in the MSSM can be suppressed compared to the SM case 
by a factor of 5. When m~ hi is increased (reduced) the effect becomes less (more) striking. 

In the case where the decoupling limit is not yet reached, the hWW and htt couplings 
are smaller than in the SM, and both the gg — > h cross section and h — > 77 widths are 
suppressed compared to the SM case, even in the absence of the squark loops. Including 
light t contributions will further decrease the amplitudes for large A t . For large tan/3 
values, the hgg amplitude can be enhanced by the 6-loop contribution, but the h — > 77 
branching ratio is strongly suppressed due to the absence of the PF-loop and the increase 
of the total decay width oc ml tan 2 (3. 

Finally, in the case of the heavy H boson, squark loop contributions to the cross sec- 
tion gg — > H can be even larger since because of the larger value of Mh, more room will 
be left for the t [and b] squarks before they decouple form the Hgg amplitude. In addition, 
for M H values above the squark pair threshold, the decays H — > tit\ or H — > b\bi will 
be kinematically allowed and could have large branching ratios, therefore suppressing the 
other decay modes including the H — > t + t - channel. For the pseudoscalar A boson, how- 
ever, squark loops will not have drastic effects on a(gg — > A): because of CP-invar iance, 
the A boson couples only to tit 2 or b\b 2 pairs while the gluon coupling to different squarks 
is absent; the Agg amplitude cannot be built at lowest order by scalar quark loops. 
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4.1.3 Higgs production in association with light stops 



If one of the stop squarks is light and its coupling to the h boson is enhanced, an addi- 
tional process might provide a new important source for Higgs particles: the associated 
production with t states, 

PP ~~ * 99 + QQ ~~ * tth (4.6) 

At lowest order, i.e. at 0(Gpa^), the process is initiated by 12 Feynman diagrams: 10 
diagrams for the gg mechanism [including those with the quartic gluon-squark interaction 
and the three-gluon vertex] once the various possibilities for emitting the Higgs boson 
from the squark lines and the crossing of the two gluons are added and 2 diagrams for 
the qq annihilation process; some generic diagrams are shown in Fig. 7. Due to the larger 
gluon luminosity at high energies, the contribution of the gg-fusion diagrams is much 
larger than the contribution of the qq annihilation diagrams at the LHC. 



g 




Figure 7: Generic Feynman diagrams for the production of the h boson in association 
with top squarks via gg fusion and qq annihilation. 

In Fig. 8, the pp — > tit\h cross section [in pb] is displayed as a function of for 
tan/3 = 2, in the case of no-mixing [A t = 200, /i = 400 GeV], moderate mixing [A t = 500 
and \x = 100 GeV] and large mixing [A t = 1.5 TeV and fi = 100 GeV]. We have used 
mi L = mi R = irtq as is approximately the case in GUT scenarios and for illustration, 
tan/5 = 2 and 30. Note for comparison, that the cross section for the standard-like 



pp — > tth process is of the order of 0.6 pb for M h ~ 100 GeV |[74[ ; m t = 175 GeV, and 
the CTEQ4 parameterizations of the structure functions [[H| are chosen. 

If there is no mixing in the stop sector, t\ and have the same mass and approximately 
the same couplings to the h boson since the mf/M| components are dominant. The cross 
section, which should be then multiplied by a factor of two to take into account both 
squarks, is comparable to the a(pp — ► tih) in the low mass range mf < 200 GeVQ. For 

16 If the t masses are related to the masses of the light quark partners, rriq, the range for which the 
cross section is rather large is therefore ruled out by the experimental constraints on m§ p5| . 
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intermediate values of A t the two components of the ht\t\ coupling interfere destructively 
and partly cancel each other, resulting in a rather small cross section, unless m t - ~ (9(100) 
GeV. In the large mixing case A t ~ 1.5 TeV a(pp — > t\tih) can be very large. It is above 
the rate for the standard process pp — > ith for values of m t - smaller than 220 GeV. If t\ is 
lighter than the top quark, the i\iih cross section significantly exceeds the one for ith final 
states. For instance, for = 140 GeV corresponding to M h ~ 76 GeV, <r(pp — > iiiih) 
is an order of magnitude larger than a(pp —>■ tih). 



A t = 1.5 TcV, /( = 100 GeV 

A t = 0.6 TeV, II = 100 GeV ■ ■ ■ ■ 
A t = 0.2 TcV, II = 400 GeV 




Figure 8: The production cross section a(pp — > t\tih) [in pb] as a function of the t\ mass 
and three sets of A t and \i values and tan/3 is fixed to tan/3 = 2. 



In Fig. 9, we fix the lightest top squark mass to = 165 GeV ~ m^ s and display 
the pp — > gg + qq — > titi/i cross section as a function of v4 t . For comparison, the * and 
• give the standard-like pp — > ith cross section for = 100 GeV and tan/3 = 2 and 
30, respectively. For tan/3 = 30 the cross section is somewhat smaller than for tan/3 = 2, 



a mere consequence of the increase of the h boson mass with tan/3 As can be seen 



again, the production cross section is substantial for the no-mixing case, rather small for 
intermediate mixing [becoming negligible for A t values between 200 and 400 GeV], and 
then becomes very large exceeding the reference cross section for values of A t above ~ 1 
TeV. For instance, for the inputs of Fig. 8, a{pp — > iii\h) exceeds a(pp — > tih) in the SM 
for the same Higgs boson mass when A t > 1(1.05) TeV for tan/3 = 2(30). 

For the signal, in most of the parameter space, the stop decay is t± — > b\ + , if < 
m t + m x o where Xi is the LSP, or i\ — > txl, in the opposite case [see section 5.4]. In the 
interesting region where a(pp — > titih) is large, i.e. for light ti, the decay i\ — > bx + is 
dominant, unless — m x + is very small, in which case the loop induced decay, i\ — > cxI, 
can become competitive. Assuming that sleptons are heavier than the chargino, xt w iU 
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Figure 9: The production cross section a(pp — * t\tih) [in pb] as a function as a function 
of A t for fixed m fl = 165 GeV and for tan/5 = 2, 30 (b). 



mainly decay into bW + + missing energy leading to t\ — > 6H^ + final states. This is the 
same topology as the decay, t — > &H^ + , except that in the case of the t there is a large 
amount of missing energy. If sleptons are also relatively light, charginos decays will also 
lead to l^Xi fi na l states. The only difference between the final states generated by the 
tth and tih processes, will be due to the softer energy spectrum of the charged leptons 
coming from the chargino decay in the former case, because of the energy carried by the 
invisible LSP. 

The Higgs boson can be tagged through its h — > 77 decay mode. In the decoupling 
limit, and for light top squarks and large A t values, the branching ratio for this mode can 
be substantially enhanced compared to the SM Higgs boson as discussed in the previous 
subsection. Therefore, 77+ charged lepton events can be more copious than in the SM, 
and the contributions of the pp — > tth process to these events can render the detection of 
the h boson much easier than with the process pp — > tih alone. 

Although a detailed Monte-Carlo analysis will be required to assess the importance 
of this signal and to optimize the cuts needed not to dilute the contribution of the tth 
final states, it is clear that in a substantial area of the MSSM parameter space, the 
contribution of the top squark to the r y r jl ± signal can significantly enhance the potential 
of the LHC to discover the lightest MSSM Higgs boson in this channel. This would be 
a new and very interesting means to search for top squarks at the LHC, which due to 
the large QCD background from it production, are otherwise difficult to detect in other 
channels. Last but not least, and as welcome bonus, this process would allow to measure 
the htt coupling, the potentially largest electroweak coupling in the MSSM, opening thus 
a window to probe directly the trilinear part of the soft-SUSY breaking scalar potential. 
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4.2 Higgs boson decays into SUSY particles 



4.2.1 SUSY decays of the neutral Higgs bosons 

As discussed previously, in mSUGRA scenarii the MSSM Higgs bosons H, A, are 
rather heavy and degenerate in mass, while the lightest h boson is in the decoupling 
regime. We have seen in the previous section 4.1 that the search of the h boson at the 
LHC relies heavily on the rare h — > 77 decay mode, while the almost only possibility 
of detecting the heavy H, A bosons for masses beyond 300 GeV and large enough tan/3 
values, tan/5 > 5, would be to look for the decays into tau pairs, H/A — > t + t~. With 
high enough luminosity, these signals will be visible at the LHC [see Fig. 4]. 

However, in these analyses, it is always assumed that the heavy H/A bosons decay 
only into standard particles, and that the SUSY decay modes are shut. But for such 
large values of Mh,a, at least the decays into the lightest neutralinos and charginos, and 
possibly into to light t and b, can be kinematically allowed. These modes could have large 
decays widths, and thus could suppress the H/A — > t + t~ branching ratios drastically. 
For the h boson, because of its small mass, only a little room is left for decays into SUSY 
particles by present experimental data |55| . However, the possibility of h decays into 



neutralinos is not yet completely ruled out, especially if one relaxes the gaugino mass 
unification; decays into sneutrinos are also still possible. When these invisible decays 
occur, they can be dominant, hence reducing the probability of the h — > 77 decay to 
occur. These SUSY decays should therefore not be overlooked as they might jeopardize 
the detection of the Higgs particles at the LHC. 

a) Invisible decays of the h boson 

Despite the lower bound of 91 GeV on the mass of the lightest chargino xt an d the 
constraints from X0X0 searches at LEP2 the decay of the lightest h boson into a pair 



of lightest neutralinos is still kinematically possible. Even in the constrained MSSM with 
a common gaugino mass at the GUT scale, leading to the well-know relation between 
the wino and bino masses Mi = |tg 2 #vyM2 ~ 1^2, the lower bound on the LSP mass is 



only m x i > 40 GeV ||55|| . Since the upper bound on the lightest h boson in the MSSM is 
M h ~ 130 there is still room for the invisible decay h — > X0X0 to occur. 

Although high values of tan/3 are required to be closer to the upper bound of Mh, the 
hbb coupling is SM-like if the h boson is in the decoupling regime: in this case the hxlxl 
coupling can be much larger then the hbb Yukawa coupling and the decay of the h boson 
into the lightest neutralinos can be dominant, resulting in a much smaller BR(/i — > 77) 
than in the SM. Far from the decoupling limit, the coupling ghbb ~ tan/5 is strongly 
enhanced for tan/5 > 3, while the h boson couplings to W bosons and top quarks [which 
provide the main contributions to the /177 loop vertex] are suppressed. This again will 
result in a strong suppression of BR(/i — > 77). 

The partial width for the decay h — > XoXo is given by [f3^ = 1 — 4m^ /M| and the 
coupling gh xx is given in eq. ( p.22|) of section 3.1] 
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The decay is important only for moderate values of M 2 and /i [with a preference for 
/i > 0] since the ft, boson prefers to couple to neutralinos which are a mixture of gauginos 
and higgsinos. In this range, the decay h — > XoXo is dominant if Mh is above the 2m x i 
threshold; close to the threshold, the decay width is strongly suppressed by the 0^ factor. 

As an illustration of this possibility we show in Fig. 10 the fraction BR(/i — > 77) 
as a function of \i for two values of tan/5 = 2, 30. We choose Ma = 1 TeV [to be in 
the decoupling regime] and the "maximal mixing" scenario A t = y/6m,q with the common 
squark mass parameter fixed to 1 TeV [to maximize the h boson mass]; this gives M h ~ 126 
GeV for tan/3 = 30 and M h ~ 106 GeV for tan/3 = 2 [the variation with /x is almost 
negligible]. Focusing first on the tan/3 = 30 and M 2 = 140 GeV case, for \/j,\ > 200 GeV 
the channel h — > X1X1 is kinematically closed and the 77 branching ratio is SM-like, BR 
~ 2.3 x 10~ 3 . In the range 110 < < 200 GeV, the LSP is lighter that M h /2 while 
the chargino is still heavier than 91 GeV, the decay h — > X1X1 is thus allowed to occur 
and suppresses the 77 branching ratio. The suppression is stronger with decreasing 
since on the one hand the phase-space becomes more favorable, and on the other hand 
the LSP tends to be an equal mixture of higgsino and gaugino. The maximum drop of 
BR(/i — > 77) is a factor of three and two for /1 > and /i < respectively. For values 
\fi\ < 110 GeV, m ± exceeds its experimentally allowed lower bound. 
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Figure 10: Branching ratios in units of 10 
tan/3 = 2(30) and M 2 = 100(140) GeV. 



for the decays h — > 77 as a function of /i for 



In the case tan/3 = 2 and M 2 = 100 GeV, the only experimentally allowed region is 
fj,\ > 110 GeV with /i < 0, since elsewhere the chargino is heavier than 91 GeV. In this 
fi\ range, the decay h — > X1X1 is kinematically allowed, but the branching ratio is very 
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small, less than 0.5%. This is due to the fact that in this area xl is a pure bino state and 
its couplings to the h boson are strongly suppressed. What makes the h — * 77 branching 
ratio drop by almost a factor two compared to the previous case is first, the smaller mass 
of the h boson [the decay width grows with the third power of the Higgs mass] but also 
because of the contribution of the chargino loops to the h — > 77 decay which interfere 
destructively for \i < with the dominant contribution due W loops [the reduction is 
nevertheless very mild, at most 15% in this case]. 

If the constraint on the unification of the gaugino masses at the GUT scale is relaxed, 
there is practically no lower bound on the LSP mass. Indeed, for relatively large values of 
the Higgs-higgsino mass parameter /i, the lightest chargino xt an d the next-to-lightest 
neutralino x 2 are wino-like with a mass ~ M 2 while the lightest neutralino is bino-like 
with a mass ~ M x . Since Mi is a free parameter, it can be as small as possible leading 
to a possibly very light LSP. The decay h — > xlxl wm then have more room to occur. 
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Figure 11: Branching ratios for the decays h — > 77 in units of 10 3 as a function of M1/M2 
for tan/5 = 2 and two sets of M 2 , fi values. 

The branching ratio for the decay h — > X1X1 can be rather large thus suppressing the 
77 branching ratio. This is exemplified in Fig. 11, where BR(/i — * 77) is shown as a 
function of the ratio M 1 /M 2 for tan/3 = 2 and for two sets of M 2 and /i values; M 2 = 140 
GeV, fi = 220 GeV leading to m x + ~ 96 GeV, and M 2 = fi = 250 GeV leading to 

m + ~ 175 GeV: the remaining inputs are as in the previous figure. In the first scenario 

1 

[solid lines], when the LSP is very small BR(/i — > 77) drops to the level of 5.10 -4 , a strong 
reduction compared to the expected rate ~ 2.10~ 3 . With increasing M\jM 2 and hence 
with increasing LSP mass, it stays almost constant until the 2m x o threshold is reached 
for Mi ~ M 2 /2 and the branching ratio recovers its standard value. 

In the second scenario [dotted lines], BR(h — > 77) starts at the same level as previously, 
but increases more rapidly and reaches approximately the standard value for Mi ~ M 2 /4 
which corresponds to the kinematical limit for the decay h — > XiXv When the LSP decay 
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is shut, the difference between the 77 branching ratios in the two scenarios is due to a 
constructively interfering chargino loop contribution [the sign of the xt contribution goes 
with the sign of fi] in the case where m x + ~ 96 GeV and which enhances the 77 decay 
width by 20% or so. This picture is expected not to be altered significantly for larger 
values of tan/3 if the h boson is in the decoupling regime as discussed previously [in fact 
for large tan/3 values and for some moderate values of the parameters M2 and /1, even the 
decays into the lightest and the next to lightest neutralinos is possible]. 

Another kinematically still possible SUSY mode for the lightest h boson is the decay 
into sneutrinos. Indeed, the experimental lower bound on the v masses is still rather 
low, my > 45 GeV |55|, leaving some room for the decay h — > vv to occur. However, 
because of SU(2)l invariance, the sneutrino and the left-handed charged slepton masses 
are related and one should avoid being into conflict with the experimental bounds on the 



II mass which are stronger, 



mr 



> 70 GeV. However, even in this case one can obtain a 



rather light sneutrino since a splitting between the v and II masses can be generated by 
the D-terms. Indeed, recalling eqs. (|3.11| ) and denoting the common scalar mass by m, 
one has: 



ml ~ m 2 + 0.50M|cos2/3 



TTlr 



~ m 2 - 0.27M|cos2/? 



(4.8) 



For small values of m, the slepton masses are governed by the D-terms, and for large 
values of tan/?, cos 2(3 — > — 1 and the D-terms become maximal. Since they tend to 
increase m~ and decrease m„, relatively low masses for sneutrinos can be kept while still 
having rather heavy left-handedQ sleptons [note however, that the v should not be lighter 
than the lightest neutralino which is expected to be the LSP]. 

In the decoupling limit, the h boson coupling to sneutrinos is also proportional to 
cos 2/3 [eq. ( |3.21| ) of section 3.1], and for large tan/3 values it becomes also maximal. And 
since it is a "gauge" coupling, it is much larger than the hbb Yukawa coupling, and the 
decay h vv is always largely dominating once it is kinematically allowed. The partial 
decay width for the decay, summing over the three sneutrinos, is given by 
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(4.9) 



Modulo the velocity factor f3„, the partial width is larger than the otherwise dominant 
bb decay width by a huge factor: M A z /{2m 2 b M 2 ) ~ 230 for M h = 130 GeV. Thus, if the 
h — ► vv decay mode is allowed, all the branching ratios for the other decay channels 
including the h — > 77 mode, will be suppressed by two orders of magnitude! Since the 
sneutrinos will decay invisibly in this mass range [m^, < vn x ± and the only possible mode 
is z> — ^ vxi], the h boson would be then also very difficult to detect at the LHGQ. 

17 The D-terms for right-handed charged sleptons are approximately the same as for the left-handed 
ones and tend also to decrease the mass. However, in GUT scenarii such as mSUGRA, the In tends to be 



lighter than the sneutrinos for reasonable values of the gaugino mass my 2 ; see eqs. ( 3.11 ). In this case 
the decay h — > vv is forbidden because of the experimental bound mr > 70 GeV. 

18 At e + e~ colliders missing mass techniques allow for an easy detection in the process e + e _ — > hZ. 
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b) H/A decays into SUSY particles 



If the CP-even and the CP-odd Higgs bosons H and A are heavy, M Hj a ^ 300 GeV, 
at least the decays into the lightest neutralinos and possibly charginos should be kine- 
matically allowed. If the couplings to bb and t + t~ pairs [which together with ti states 
account for most of the total decay width in the absence of SUSY modes] are not strongly 
enhanced, and hence for not too large values of tan /3, these decays might be dominant and 
suppress drastically the branching ratios for the H/A — > t + t~ signals. This is exemplified 
in Fig. 12 where BR(H/A — > t + t _ ) are plotted as a function of the H/A masses for three 
values tan/5 = 5, 10 and 30. The choice M 2 = —pi = 200 GeV has been made leading to 
m v o ~ 90 GeV and m + ~ 160 GeV [with a small variation with tan/31. 



o.oi 



1 1 1 


l l 


i 

H 

A 


■ .:. tan/3 = 


30 




' ■ . . \. tan/3 = 


10 




tan/3 = 


5 




BR(H/A -> t+t-) 






M 2 = -fi = 200 GeV 
i i i 


i i 


i 



300 400 500 600 700 800 900 1000 



Figure 12: Branching ratios for the decays H/A — > t + t as a function of M H:A for 
tan/3 = 5, 10 and 30 and for the values M 2 = -/j, = 200 GeV. 



The branching ratios for H and A decays are almost the same except for small values 
of tan/3 and relatively small Higgs masses: in this case, the decoupling limit is not yet 
reached and additional [and different] decay modes occur for the H and A bosons; see 
section 4.1. For tan/3 = 5 the H/A couplings to down-type fermions are not very strongly 
enhanced and the decays into charginos and neutralinos have large branching ratios: they 
decrease BK(H/A — > t + t~) from the standard ~ 10% value for small Higgs masses 
[where only a few SUSY channels are open and some are suppressed by phase space] 
to less than 0.4% for very heavy Higgs boson masses M h ,a ~ 1 TeV [here most of the 
neutralino/chargino channels are open and they are not suppressed by phase space], thus a 
reduction by more than a factor of 20 compared to the branching ratio without the SUSY 
decays. For tan/3 = 10, the couplings to 6-quarks and r-leptons are more enhanced and 
BR(H/A — > t + t~) are larger by slightly more than a factor two compared to the previous 
case. For even larger values of tan/3, tan/3 = 30, the decays into charginos and neutralinos 
are not dominating anymore, and the branching ratios for the H/A decays into tau pairs 
are suppressed only slightly, less than a factor of two. The branching ratios for the decays 
into ji + [i~ can be obtained from Fig. 12 by rescaling the numbers by m 2 /m 2 .. 
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In the preceding discussion, the decays of the H and A bosons into sfermions were 
assumed to be shut. However, as discussed previously, at least one of the stops can be 
rather light and its couplings to the Higgs bosons enhanced; the decays of the heavier 
Higgs bosons H, A to stop squarks might be therefore kinematically accessible and could 
dominate over the standard decays, and even over the decays into charginos and neutrali- 
nos^|. Indeed, in the decays of the H boson into stops, the partial widths up to mixing 
angle factors are proportional to Gpmf /(M H tan 2 /?) or/and Gprn^fi — A t /tanP) 2 /M H ; 
for small tan/3 values and not too large Mh and for intermediate tan/? values and for large 
H and A t , the widths for the decays H — > it can be very large and can compete with, 
and even dominate over, the other [standard and SUSY] decay channels. The branching 
ratios for the H decays into r pairs would be then further suppressed. 
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Figure 13: Branching fractions for H — > t\t\ [left] and H — > t + t~ [right] as a function of 
M H for tan/? = 2.5, 5, 10, 30 and M 2 = fi = m f ~/2 = 250 GeV and A t = 1.5 TeV. 

This is illustrated in Fig. 13, where BR(if — > txti) is shown as a function of Mh for 
tan/? = 2.5, 5, 10, 30 and ~ 200 GeV [for tan/3 = 2.5 this is achieved by setting mj L = 
m f ~ R = 500 GeV and A t = 1.5 TeV]; M 2 = fx = 250 GeV. As expected BR(H -> tJi) 
decreases with increasing tan/3 values and increasing Mh- However, it is still at the level 
of ~ 50% for tan/3 = 5 and Mh = 1 TeV. For tan/3 = 30, the channel H — ► t\t 2 opens 
up for Mh ~ 900 GeV, and the curve shows the sum BR(if — > itfi + iii 2 )- But for this 
large tan/3 value, the branching ratio barely exceeds the level of 20% in contrast to lower 
tan/3 values where it can reach almost unity for small Mh- For larger Mh, the decays into 
charginos and neutralinos become more important and will dominate; so BR(if — > t + t~) 
is reduced anyway. For the A boson the only important decay into sfermions is A — >• tit 2 
[and maybe b±b 2 for tan/? ^> 1]. Thus both stops must be light for the decay to be 
allowed by kinematics. This happens only in a small area of the parameter space, unless 
all squarks are relatively light. For instance, in the scenario of Fig. 13, m 4 ~ 2 ~ 700 GeV 
and the decays A,H — > t\t 2 occur only for masses close to 1 TeV. 

19 The decay widths of the H bosons into the light fermion partners are proportional to 
GfM^t sin 2 2(3 /Mh for Mh ^> rnjr. They are thus suppressed by the heavy H mass and cannot com- 
pete with the decays into fermions [t, b, t and possibly \ states] for which the widths grow as Mh- The 
pseudoscalar A boson cannot decay into the partners of light fermions, if the fermion mass is neglected. 
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4.2.2 Decays of the charged bosons 



Charged Higgs bosons are searched for at the Tevatron and will be searched for at the 
LHC in the mass range M#± < m t — mb, i.e. when produced by a decaying t-quark. Some 
experimental limits from Tevatron [a tan/3-M#± exclusion contour] are already available 
HfBH - However, the detection techniques rely on the assumption that no H + decay mode 
other than cs and t + u t is kinematically significant, and in particular that no decay into 
SUSY particles is possible. However, decays into the lightest charginos and neutralinos as 
well as decays into sleptons [|1(J are still allowed by present experimental data; for instance, 
the decay channel H + — > XiXi is certainly possible for M H ± > 130 GeV, despite the lower 
bound of m x + > 91 GeV |5|. 

In Fig. 14a, the relative branching ratio for the decay H + — ► xfXi is shown as a 
function of M#+, for tan/3 = 2, M 2 = 110 GeV, \fi\ = 500 GeV, and all remaining soft 
masses large enough to forbid other supersymmetric decay modes. For these values of 
parameters, one has m x + = 96.5 GeV and m x o ~ 50 GeV. Also shown is the case where 
the assumption of gaugino mass universality is relaxed: keeping fixed all other parameters 
to the previous values, M\ is set to 20 GeV, and the lightest neutralino mass becomes 
m x o ~ 16 GeV. In this case, the decay channel XiXi opens already at ~ 115 GeV. For 
small values of tan/3, whenever kinematically allowed, the mode XiXi is the dominant 
one. For large values of tan/3, however, only the channel H + — > r + u T survives. 
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Figure 14: Branching fractions for the decays H + — > XiXi (left) and H + — > l + v (right) 
as a function of M^±, for tan/3 = 2, and some M2, Mi, [i and A values. 

The previous discussion applies only to supersymmetric models with heavy sfermion 
masses. The existing lower bounds on slepton masses, however, are still rather modest: 
~ 65 GeV for f's and > 45 GeV for z>'s |]55|| . The decay H + — > f + v T is therefore 
kinematically possible for M H ± > 110 GeV. Fig. 14b shows the relative branching ratio for 
tan/3 = 2, M2 = 120 GeV, Mi ~ M2/2 and two choices of parameters in the slepton mass 
matrices: (a) mj L = mj R = mi = 75 GeV, /i = 500 GeV, A T = and (b) mj = 90 GeV, 
\i = —500 GeV, A T = 1 TeV. This leads to a slepton spectrum: (a) m s = 56 GeV, 
rrig ~ rrijx = 83 GeV and the two f masses are 10 GeV below and above this value, and 
(b) Tfiy = 75 GeV, m g ~ = 97 GeV and = 63 GeV, and m f2 = 121 GeV. 
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The prominence of the H + — > f + z/ T mode above threshold is explained by the fact that 
the H + u* L l L coupling oc g/\^2Mw sin 2/3 is very large compared to the H + u* l t r coupling 
g / v2(m T / My/) . Due to the sin 2/3 dependence, this term quickly dies off for increasing 
tan /3. In this case, however, there exists other directions of parameter space where this 
decay mode has still a branching ratio close to unity When A T ~ /itan/3, in fact, the 
left-right mixing in the slepton mass matrix tends to vanish, but the coupling of the 
Lagrangian term H + v* L f R : —g/V2(m T /M w )(fi+A T tanf3) acquires a 1/ cos 2 (3 dependence, 
which increases with increasing tan/3. For instance, tan/3 = 10, A T = 2 TeV, /i ~ 200 GeV 
and mi ~ 80 GeV, give a branching ratio above 90% for the decay channel H + — > f + u T , 
when kinematically accessible. Such a decay mode produces a final r + plus missing energy 
as the r + z/ r mode, but with a much softer r + than that searched for at hadron colliders. 

Finally, let us make a few remarks on other possible decay of in non-SUSY two- 
Higgs doublet models (2HDM). Although the decay b — > 57 excludes masses M H ± up 
to ~ 165 GeV, irrespectively of tan/3 JF7], direct searches still allow for the possibility 



Mh± < Myy- The states are searched for at LEP2 and the Tevatron, relying again 
on the cs and rv signals [see Ref. for details and references]. 



(i) In 2HDMs, there is no lower bound on the pseudoscalar mass Ma coming from 
LEP: since the mixing angle a is in this case a free parameter, it can be chosen in such a 
way that the coupling ZhA vanishes and the process Z — > hA does not occur. This makes 



that the LEP2 bound M4 > 75 GeV |55j does not hold, leaving open the possibility of 
a very light A bosonQ. [In this case the cross section for the process e + e~ — > Z* — > hZ, 
which is complementary, is not suppressed compared to the SM Higgs boson, and the 
bound M h > 88.6 GeV |55] applies here]. Now since gn+w-A is a "gauge" coupling, it 
is clear that the decay H + — > AW + can be rather important in a 2HDM. Indeed, this 
remains true even for an off-shell P^-boson, in spite of the virtuality and the additional 
weak coupling suppressions. Even for masses Mn± ~ 55 GeV, i.e. roughly the value 
excluded at LEP2, the branching ratio can be still at the level of 50% for small enough 
tan/3 and Ma values. For heavier if ± bosons that can be searched for at the Tevatron, 
the decay tf* -> AW& can be by far the dominant one. 

(ii) For charged Higgs boson masses above ~ 140 GeV, even if the decays into AW* 
bosons are suppressed [for instance the A boson is too heavy], the H + — > tv t and cs 
channels are still not the dominant ones for small tan/3 values. Indeed, for tan/3 ~ 1 the 
three-body decay mode H + — > bt* — > bbW + is already the dominant decay mode: despite 
the virtuality of the top quark and the fact that the process is of higher-order, the H ± 
coupling to tb quarks is much larger than the rv and cs couplings, leading to a large 
compensation. For instance, for M H ± = 140 GeV and tan/? = 1, the branching ratio is 
already at the level of 50%. Note that this decay mode [as well as the three-body decay 
mode H^ 1 — > hW* — > hff] can also be important in the MSSM for large enough H 
masses [but still below m t — nib] and small tan/3 values. 

20 Two other production mechanisms are possible at LEP1, since they require very large statistics: the 
fcrmion loop-mediated decay Z — > and the radiation off the bb and t + t~ lines; however the rates are 
small especially for moderate values of tan/3. One remains therefore with the rather modest bound from 
the decay T — > A7 which has been searched for by the Crystal Ball Collaboration Q , Ma > 5 GeV. 
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4.3. MSSM Higgs production in e + e~ collisions 
4.3.1 Production mechanisms 

The main production mechanisms of the neutral MSSM Higgs bosons at e + e~ colliders 
are the Higgs-strahlung process and pair production J7j| |50j , as well as the WW and ZZ 



fusion processes |68 



(a) Higgs-strahlung e + e _ — > (Z)^Z + h/H 

(b) pair production e + e _ — > (Z)^A + h/H 

(c) fusion processes e + e" — > vv (WW} — > vv + h/H 

e + e- -> e + e~(ZZ) -> e + e~ + h/H (4.10) 

The CP-odd Higgs boson A cannot be produced in the Higgs-strahlung and fusion pro- 
cesses to leading order since it does not couple to vector boson pairs. 



The charged Higgs particle [5(| can be pair produced through virtual photon and 
Z boson exchange in e + e~ collisions, and also in top-quark decays for masses below 
m t — m b ~ 170 GeV: 

(d) charged Higgs e + e~ — > (7, Z*) — > H + H~ 

e + e~ -> tt with t H + b (4.11) 

The production cross sections for the neutral Higgs bosons are suppressed by mixing 
angle factors compared to the SM Higgs production, 

a(e + e- -> Zh) , a(VV -> /i) , a(e + e" -> Aff) ~ sin 2 (/3 - a) 
a(e+e- -> ZH) , (r(yy -> if) , a(e + e- -> A/i) ~ cos 2 (/3 - a) (4.12) 

while the cross section for the H ± particles does not depend on any parameter other than 
Mh± when pair produced in e + e~ collisions; a tan/3 dependence is however present for 
H ± production in top decays, due to the branching ratio BR(i — ► H + b). 

Modulo the mixing factors, the cross section for the fusion process, e + e~ — > i> e v e § 
with $ = h, H is enhanced at high energies since it scales like log s/M|; the cross 
section for e + e~ — > e + e~$ is ~ 16cos 2 ^vy i- e - one order of magnitude smaller. The 
cross sections for the Higgs-strahlung and pair production processes, e + e~ — > $A 
and H + H~ scale like 1/s and hence are smaller at high energies. Close to the decoupling 
limit, the factor cos(/5 — a) — > M\ sin2/3/(2M^) vanishes, and the only relevant processes 
for the production of the heavy states [when kinematically allowed] are the pair production 
of A and H and the charged Higgs pair production: e + e _ — ► HA and e + e~ — > H + H~ . 
The cross section for the fusion process, e + e~ — > uuH, which in principle is enhanced at 
high energies, is only relevant in the mass range of a few hundred GeV and for small tan/? 
values where the factor cos 2 (/3 — a) is not prohibitively small. For the lightest h boson, 
the only remaining production processes are the bremsstrahlung and vector boson fusion 
mechanisms: e + e~ — > hZ and e + e~ — > vv/e + e~ + h, exactly like for the SM Higgs bosons. 
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The cross sections for processes the (a)-(c) and (d) are shown in Figs. 15 and 16 as 
functions of the Higgs boson masses for tan/3 = 2.5 at a cm. energy yfs = 500 GeV. 
They are obtained with the help of the program HPROD to be discussed later. 

There are several additional processes for the production of the neutral Higgs particles, 
and in particular of the lightest h boson, that we will take as an example. Although these 
processes are of higher order, the production rates can be substantial for high-luminosities, 
and they might open a window to the determination of several fundamental properties of 
the Higgs particles. 

(i) The associated production with tt pairs which allows to measure directly the 



tt-Higgs Yukawa coupling, the strongest coupling in the standard sector: 

e + e" -> ( 7 ,Z*) -> ti<!> (4.13) 

(ii) The double Higgs boson production, either in the bremsstrahlung or in the fusion 
processes, which allows to determine the trilinear Higgs boson couplings, and reconstruct 
the all important scalar potential j82], ^3 



e + e- -> Zhh and WW/ZZ -> hh (4.14) 



(Hi) The production of the Higgs bosons in association with a photon [84], 85J , which 
allows to measure the /177 and hZj loop induced vertices and probe the effects of heavy 
particles which couple to the Higgs boson 

e+e" -> hj (4.15) 

In addition, if stop squarks are light, the production of the h boson with stop pairs, 
e + e~ — > tth ||86|| , which allows to measure the tth couplings, the largest electroweak 
couplings in the MSSM, opening a window to the probing of trilinear couplings. 

Finally, e + e~ linear colliders can be made to run in the 77 mode with the high energy 



photons coming from Compton back-scattering |87j of laser photons. This lead to the 
production of the Higgs bosons as s-channel resonances: 77 — > h, H, A, allowing to mea- 
sure the Higgs-photon couplings [similarly to the process e + e~ — > 7+Higgs]. Furthermore 
charged Higgs bosons can be pair-produced in this mode. 

4.3.2 The program HPROD 

There are many generators which deal with the production of the MSSM Higgs bosons at 
e + e~ colliders Here, we briefly describe the program HPROD [[TTJ which calculates 



the cross sections of the SM and MSSM Higgs particles at e + e machines. It includes: 

(i) In the SM: the strahlung and the WW/ZZ fusion processes as well as the higher 
order processes e + e~ — > ttH°, H°H°Z and VV — > H°H° in the longitudinal vector boson 
approximation. The associated production with a photon will be included soon. 

(ii) In the MSSM it includes all the processes (a)-(d) above: the strahlung and the 
pair production as well as the WW/ZZ fusion processes for neutral Higgs production, 
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and the e + e and top decay processes for the charged Higgs boson. The higher order 
processes will be included soon. 

(Hi) In the MSSM, the complete radiative corrections in the effective potential ap- 
proach with full mixing in the t,b sectors; it uses the RG improved values of the masses 
and couplings, and the relevant leading next-to-leading-order corrections are also imple- 
mented. Both the pole and the running Higgs boson masses are calculated. 

(iv) Off-shell Z and Higgs bosons in the bremsstrahlung and the Higgs pair production 
processes [the latter only in the region of parameter space where they are relevant: close 
to the production threshold and for sizeable Higgs boson decay widths]. The Higgs boson 
total widths [without the decays into SUSY particles] are included in the routines. 

(v) The initial state radiative corrections for the bremsstrahlung, pair production 
processes and the WW/ZZ fusion processes. The effect of beamstralung for the TESLA 
machine by making an interface with the program CIRCEE |88| will be implemented soon. 

The basic input parameters, fermion and gauge boson masses and their total widths, 
coupling constants and in the MSSM, soft SUSY-breaking parameters can be chosen from 
an input file. In this file several flags allow to switch on/off or change some options [e.g. 
choose a particular Higgs boson, include/exclude initial state radiation, use running or 
pole Higgs masses]. The parameters in the input file are: 

y/s the cm. energy of the collider; IHIGGS: an integer which chooses the Higgs 
boson [0,1,2,3 for #°, /i, if, i/^ and 4 for all MSSM Higgs bosons]; tan/? and M A : the 
two basic parameters in the MSSM Higgs sector; m b , m t the bottom and top masses; 
a s{Mz), Gp, a(0): the strong, Fermi and QED coupling constants; M z , M W ,T Z - the 
W, Z boson masses and the Z total width; M 2 , A*, m,g , m^- , m_g : the SUSY break- 
ing mass parameters; Ab, A t : SUSY breaking trilinear coupling for stops and sbottoms; 
ISHELZ/ISHELH: integers which for 1(0) include (exclude) off-shell Z/Higgs bosons in the 
bremsstrahlung/pair production processes; IPOLE: an integer which for 0(1) calculates 
running (pole) MSSM Higgs masses; ISR: an integer which for 1(0) excludes (includes) 
initial state radiation; IHIGH: an integer which for 0(1) includes (excludes) the higher 
order processes for the SM Higgs boson. An example of input file for Higgs boson pro- 
duction in the MSSM is shown below. 



SQRS 




500. DO 


ALPH 




137.036D0 


MDR 


= 1000 


DO 


IHIGGS 




4 


GF 




1 . 1664D-5 


AU 


= 2400 


DO 


TGBET 




2.5D0 


GAMZ 




2 . 489D0 


AD 


= 2400 


DO 


MABEG 




50. DO 


MZ 




91. 187D0 


IPOLE 


= 




MAEND 




250. DO 


MW 




80 . 33D0 


ISHELZ 


= 




NMA 




5 


MU 




500. DO 


ISHELH 


= 




MB 




4.62D0 


M2 




500. DO 


ISR 


= 1 




MT 




175. DO 


MSG) 




1000. DO 


IHIGH 


= 1 




ALS(MZ) 




0.118 


MUR 




1000. DO 









Tab. 3: Example of input file for the MSSM Higgs production at yfs = 500 GeV. Figs. 15 
and 16 are obtained with the inputs given in this file. 
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Figure 15: Production cross sections for the neutral MSSM Higgs bosons in e + e collisions 
as a function of the h and H masses at a cm. energy y/s = 500 GeV and for tan/? = 2.5. 
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Figure 16: Production cross sections for the charged Higgs boson in e + e collisions as 
a function of the H ± mass for a cm. energy y/s = 500 GeV and for the decay t — > 
H + b, tan/3 = 2.5. 
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The program is written in FORTRAN and has been tested on several machines. All 
the necessary subroutines [e.g. for integration] are included. The program is lengthy but 
rather fast [especially if some options, e.g. ISHELH, are switched off]. The results for the 
many production cross sections in femtobarns are written into several output files [with 
headers indicating the processes and giving the inputs]: 

zps.out : M H o , a(H°Z) , a(H°uu) , a(H°e + e-) 

zpt.out : M H o , a(WW -> H°H°) , a(ZZ -> H°H°) , a(H°H°Z) , a(tiH°) 

zpl.out : M A , M h , a(hZ) , a(huu) , a{he + e~) , cr(hA) 

zph.out : M A , M H , a {HZ) , a(Hvv) , a(He + e~) , a (HA) 



zpc.out : 


M A , M H± 


, <r(e + e — 


H + H~) , a(tt - 


-> H + H-) 




MA 


Mh 


hZ 


hnunu 


he+e- 


h A 


50.0000 


48.0416 


22.09 


52.00 


5.483 


32.98 


100.000 


79.7747 


37.30 


72.05 


7.658 


17.48 


150.000 


97.2542 


52.94 


91.32 


9.741 


4.564 


200.000 


103.804 


57.82 


95.56 


10.21 


1.019 


250.000 


106.497 


58.98 


95.78 


10.23 


. 2473 




Tab. 4a: The output file zpl.out for the input file in 


Tab. 3. 


50.0000 


131.191 


37.36 


46.92 


5.032 


13.66 


100.000 


143.719 


22.15 


25.60 


2.750 


20.38 


150.000 


173.380 


6.506 


6.176 


0.6653 


20.11 


200.000 


215.038 


1.625 


1. 164 


0.1257 


8.726 


250.000 


261.096 


0.4601 


. 2384 


0.2575E- 


01 0.000 




Tab. 4b: The 


output file zph.out for the input file in 


Tab. 3. 


MA 


MH+ 


ee->H+H- 


tt->H+H- 






50.0000 


92.1414 


98.79 


53.73 






100.000 


126.452 


78.94 


24.77 






150.000 


168.790 


49.36 


0.5294 






200.000 


214.453 


16.70 


0.000 






250.000 


261.706 


0.000 


0.000 







Tab. 4c: The output file zpc.out for the input file in Tab. 3. 

In a future version of the code, we will implement the remaining higher order processes 
[discussed previously] of the SM and for the h boson in the MSSM, and include the 77 
option of the collider. We also plan to make an interface with the programs SUSPECT 
[section 3.2] and HDECAY [which calculates the total decay widths in a more precise 
way]. Beamstralung will also be incorporated by making and interface with the program 
CIRCEE. The Higgs boson production at hadron machines will be included in the next- 
to-leading version of the code. 



70 



4.3.2 Higgs boson production in association with light stops 



As previously discussed [section 3.1 and above], if the off-diagonal entries in the t mass 
matrix is large, the eigenstate t\ can be rather light and at the same time its couplings 
to the Higgs bosons strongly enhanced. The htit\ coupling would be then the potentially 
largest coupling in the electroweak sector of the MSSM, and its measurement would 
allow to pin down some of the soft-SUSY breaking parameters of the MSSM Lagrangian, 
in particular the trilinear coupling A t . This coupling can be measured in the process 
PP — > Htih at the LHC [see section 4.1.3]. However, this coupling can be best measured 



in the associated production process e + e~ — > tit±h fll2fl . This is first due to the fact 
that the cleaner environment of high-energy electron-positron colliders allows for a more 
efficient search of this complicated final state, ans also because of the very high-luminosity 
expected at such machines [for instance for the DESY-machine TESLA: / C ~ 500 fbr 1 in 
a year] , which compensates for the fact that the process is of higher order in perturbation 
theory and has a small cross section in priciple. 

At future linear e + e~ colliders, the final state t\tih may be generated in three ways: 



a) Two— body production and decay: 

If the center of mass energy of the e + e~ collider is high enough, one first produces 
a mixed pairs of top squarks, e + e~ — > tit 2) through the exchange of a virtual Z-boson, 
and then makes the heaviest top squark decay into the lightest stop and the Higgs boson, 
t 2 — > tih, if the splitting between the two stops is larger than the h boson mass. In 
principle, if phase-space allowed, the cross section for the two-body production process 
times the branching ratio for the two-body decay should be large enough for the final 
state to be copiously produced. 

However, the Zt\t 2 coupling is proportional to sin 29 t while the ht\t 2 coupling is pro- 
portional to cos 29 1, since in a large part of the MSSM parameter space [as discussed 
in section 4.1.3] sin2# t is either small [no-mixing case] or close to one [maximal mixing 
case], the cross section times branching bratio, which is then proportional to sin4^ g , is 
always rather small. [In addition, the decay width t 2 — > ht\ is in general much smaller 
than the t 2 decay widths into chargino and neutralinos] . Nevertheless, there are limited 
regions of the MSSM parameter space where the rate for this process is visible for the 
high luminosities J C ~ 500 fb _1 expected at the linear colliders. 

b) Production in the continuum in e + e collisions: 

There are three types of Feynman diagrams leading to t\txh finale states in e + e~ 
collisions: Higgs boson emission for the t\ lines which proceeds through s-channel photon 
and Z-boson exchange, Higgs boson emission from the t 2 lines which proceeds through 
Z-boson exchange, and Higgs boson emission from a virtual Z boson which then splits 
into t\ti pairs. When the g^uu coupling is large, the two last types of Feynman diagrams 
give negligible contributions since the virtuality of t 2 is large and the Zt\t\ coupling not 
enhanced. The Dalitz plot density of the process [in terms of the scaled energies of the 
two stops] can be then written in a very simple form. 
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As an illustration, the total cross section for the process e + e~ — > ht\tx is shown in 
Fig. 17 as a function of the lightest stop mass for a cm. energy ^/s = 500 GeV, for 
the values At = 1.5 TeV and tan/3 = 3,30. One can see that below masses of the order 
of mj ~ 160 GeV, the cross section is larger than 1 fb, leading to ~ 500 events with 
the expected luminosity quoted above. The final state topologies have been discussed in 
section 4.1.3; however, here the main decay mode of the h boson, h — > bb, can be used 
[this calls for good micro-vertex detectors]. With the relatively large sample of events, 
this final state should be experimentally possible to be detected even after efficiencies 
have been included. Higher cm. energies would allow to probe larger stop masses. The 
cross section for this process has also been calculated in Ref. p9| using the package SUSY- 
GRACE; it has been shown that combined with the e + e~ — > titx production cross section, 
the SUSY parameters of the stop sector can be determined. 
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Figure 17: The cross section for the process e + e~ — > ht\ti as a function of the lightest stop 
mass for a cm. energy y/s = 500 GeV and for the values A t — 1.5 TeV and tan/3 = 3, 30. 



c) Production in the continuum in 77 collisions: 

In the 77 option of future high-energy e + e~ linear colliders, with the high energy 
photons coming from Compton back-scattering of laser beams, the final state t\tih can 
be generated by emitting the h boson from the stop lines in the process 77 — > tit±. 



In Ref. ||12|| , the cross section for this reaction has been calculated, and the analytical 
expression has been found to be slightly more complicated than the one obtained in the 
e + e~ mode. 

However, because the cm. energy of the 77 collider is expected to be only ~ 80% of 
the one of the original e + e~ machine, the process is less phase-space favored. In addition, 
the collected luminosity is expected to be somewhat smaller than the one of the e + e~ 
mode, leading to a smaller number of events. Nevertheless, the cross section for the 
t\t\h final state is of the same order as in the e + e~ mode for cm. energies not too close 
to the kinematical threshold, and the process might be useful to obtain complementary 
information since it does not involve the Z-boson and exchanges. 
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5. SUSY Particle Production and Decays 



5.1. Virtual SUSY effects 

Virtual effects of supersymmetry, for example through gauge boson self-energy or gauge 
boson-fermion-fermion vertex corrections, are known to be invisibly small at the Z reso- 
nance unless the masses of the SUSY particles are very small, a situation which is by now 
excluded or unlikely f55j. It has been noticed |K| that this property does not apply for 



energies beyond the Z peak due to the raise of box diagrams. At LEP1 energies, the box 
contributions are strongly suppressed as they do not benefit of the enhancement due to 
the Z resonance. Beyond this energy, one can expect them to become relatively impor- 
tant. This is what happens in the standard electroweak case with the box diagrams due 
to WW formation and neutrino exchanges, where their effect reaches the percent level 
in the cross section for e + e~ — > at y/s ~ 200 GeV. A second source of local box 

enhancement is the threshold effect appearing around s = (mj + mj) 2 where and rrij 
are the masses of the particles formed in the intermediate state. In the aforementioned 
WW box diagrams the contribution peaks at 1.2% around 161 GeV. 

With this standard situation in mind, one can look at the SUSY box contributions due 
to charginos or neutralinos and sleptons for the production of light fermion pairs [[HJ. The 



case of Xi ~X°-slepton boxes is particularly favored for several reasons: the usual gaugino- 
fermion-sfermion couplings have full electroweak strength, and the box contribution is not 
mixed with other significant virtual corrections [self-energy, vertex corrections] involving 
neutralino pairs since the gaugino components are decoupled from gauge bosons and their 
higgsino components are decoupled from light fermions, so it should be easier to single out 
this contribution. In order to disentangle such contributions from other effects affecting 



the Z boson tail, it is convenient to use the so-called "Z-peak subtraction method" [91 
which consists in using as inputs the measured Z-parameters [mass, partial widths and 
asymmetries] and to describe the additional contributions to the observables beyond the 
Z-peak in terms of four functions constructed with a subtraction at s = M|. In this 
way, one eliminates all known or unknown effects contributing at the Z-peak and to the 
Z boson tail, which do not have a strong s-dependence. 

A first numerical exploration, assuming that neutralinos as well as the left- and right- 



handed sleptons are degenerate, has been made [[13] . This choice allows to greatly simplify 
the computation as in this case one can use sum rules for reducing the lengthy expressions 
in the box amplitude. For simplicity neutralinos were also taken as pure wino and bino 
states. Results are shown in Fig. 18 for the e + e" — > cross section and forward- 

backward asymmetry, for a common neutralino mass m x o = 100 GeV and a slepton mass 
mg = 60 GeV. As expected, a peak associated to the threshold effect appears at s = ^rn?^, 
which reaches the level of 1.5% in the cross section and 0.7% in the FB asymmetry. The 
half-width is of about 15 GeV. 

The level reached by these contributions is very encouraging since such a precision for 



the cross section is certainly possible at LEP2 [92]; in the case of the FB asymmetry the 



observability is more difficult. Increasing the e mass to 110 GeV, the peak cross section 



73 



however decreases roughly by a factor of two. In the case of chargino boxes, assuming again 
m x + = 100 GeV and mj = 60 GeV, one obtains a similar effect on the cross section and on 
the forward-backward asymmetry [in this case, the sign is opposite]. However, in contrast 
to the neutralino case, there now exist other important chargino effects in gauge boson 
self-energies and vertex corrections; cancellations among these various contributions occur 
and depend on the masses and couplings of the SUSY particles. 

These results are promising since from a general point of view they suggest that virtual 
SUSY particle effects may be observable at LEP2. A complete study for both neutralino 
and chargino contributions at LEP2, including all possible effects within the MSSM and 
varying the parameters in their allowed domain, is under way. 
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Figure 18: Neutralino box effects on muon pair production at LEP2 energies. Relative 
effect in permille on the cross section (solid) and on the forward-backward asymmetry 
(dashed) as a function of the cm. energy. 



5.2 Correlated production and decays in e + e — ► XiXi 

Once charginos will have been discovered, the experimental analysis of their properties, 
production and decay mechanisms, should reveal the basic structure of the underlying 
low-energy supersymmetric theory. This needs precise experimental measurement, which 
are more likely to be performed in the clean environment of e + e~ colliders. In the case of 
the chargino sector, three parameters are needed to describe it completely: tan (3, /i and 
M 2 . One needs therefore several experimental observables to pin down these parameters. 

In e + e~ collisions, charginos are produced either in diagonal or in mixed pairs . 
Since the second chargino is generally expected to be significantly heavier than the first 
state, at LEP2 and potentially even in the first phase of e + e~ linear colliders, the lightest 
chargino xt may be, for some time, the only state that can be studied in detail. It is 
therefore wiser to focus on the diagonal pair production of the lightest chargino e + e~ — > 
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XiXi ■ The production process is generated by s-channel 7 and Z exchanges, and t- 
channel v exchange. It will depend on the two angles 0l and <pR [more precisely on 
cos20i and cos2<Pb\ which define the chargino-chargino-Z and the electron-sneutrino- 
chargino vertices, the chargino mass m + and the sneutrino mass m^. Thus at least four 

Xi 

measurements are needed to determine these parameters. 

However, the situation is complicated by the decays of the charginos. Indeed, the 



charginos will decay into the lightest neutralino and light fermion pairs |93[ and since 
the two neutralinos are stable and escape undetected, it is not possible to make a com- 
plete reconstruction of the events; in particular, one cannot measure the xt production 
angle. Furthermore, the chargino decays occur through W boson and sfermion [and also 
a marginal contribution from charged Higgs boson] exchanges, and the knowledge of the 
sfermion masses and couplings is in principle also required. 

Recently it has been shown JTJ]] that even in this situation the underlying SUSY 
parameters can be extracted from the mass m v ±, the total production cross section, and 

Xi 

the measurement of the polarization with which the charginos are produced. The x ± 
polarization vectors and x~X spin-spin correlation tensors can be determined from the 
decay distributions of the charginos. Beam polarization is helpful but not necessarily 
required. No detailed information on the decay dynamics, nor on the structure of the 
neutralino, is needed to carry out the spin analysis. There are already several analyses 
dealing with polarized chargino production [see Ref. [Q or the contribution of Katsanevas 
et al. |5SJ for instance]; Ref. [Q however attempts to explore analytically the event 
characteristics as will be summarized below. 

Since the xt lifetime is very small, only the correlated production and decay can be 
observed experimentally 

e+e" -> Xi(P,n)xi(p,n) 

^ X°x + (/1/2) 



- x? + (hh) 

where n and n are the spin 4-vectors and p, p the momenta of the charginos. In covariant 
language the final state distributions are found by combining the polarized cross section 



da = (da)- [l - - + Q^n^n v 



(5.1) 



with the polarized distributions for the decays into a neutralinos and light fermion pairs 



dr = (dr) 



dr = (dr) 



Inserting the completeness relations Y^ n u,n u = —Quv + PnPu/ m2 -- 

Xi 

event topology can then be calculated from the formula 

1 



(5.2) 

r]^ u etc, the overall 



da final = (da)(dr)(dr)- 



(5.3) 



In these expressions V{V) and V(V) are the polarization vectors of the produced (de- 
caying) states, while is the chargino spin-spin correlation matrix. 
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In the case of CP-invariance [the CP non-invariant case has been discussed in Ref. 
the overall topology is determined by seven independent functions: the unpolarized cross 
section cr unpo i, the two components of the polarization vector of one of the charginos, and 
four correlation functions. Except for o" unpo i, these observables will depend on the final 
state, and hence on the decays of the charginos into neutralinos and fermion pairs. How- 
ever, it was shown [14| that it is possible to construct two observables, denoted by V 2 /Q 
and V 2 /y, which do not depend on the chargino decays and hence, reflect unambiguously 
the properties of the chargino system, being not affected by the neutralinos. The energy 
dependence of the ratios V 2 /Q and V 2 /y is shown in Fig. 19 for tan/3 = 2 and for the 
(M 2 , p) parameters chosen as: gaugino region (81, —215) GeV, higgsino region (215, —810) 



GeV and mixed region (92, —93) GeV, which lead to a chargino mass m^± ^ 
The two ratios are sensitive to the couplings at sufficiently large cm. energies. 
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Figure 19: The energy dependence of the ratios V 2 j Q and V 2 /y for tan/5 = 2. 

Therefore, in addition to the chargino mass m^± which can be measured very precisely 
near the threshold where the production cross section a(e + e~ — > XiXi) rises sharply with 
the velocity, and the v mass which can be determined with the energy variation of the cross 
section [in case where the determination of mp is poor, longitudinal beam polarization 
might be helpful to eliminate the t-channel v exchange], we have two observables which 
allow to determine the two mixing angles cos 2(f) l and cos 20^. A representative example 
of determining these two parameters is shown in Fig. 20 at a cm. energy y/s = 500 
GeV for the "measured" values: m x ± = 95 GeV, <r(e + e~ — > XiXi) = 0-38 pb, V 2 /Q = 
—0.25,V 2 /y = —5.0. The three contour lines meet at a single point (cos2</> £ , cos 20^) = 
(-0.60, -0.49) for m p = 250 GeV. 

Once cos 26 L , cos 2(f> R and m ± are known, the SUSY parameters tan/3, M 2 , u can 
be then determined with at most a two-fold discrete ambiguity. Introducing the two 
triangular quantities p = cot(<pR — 4>l) and q = cot((pR + 4>l) and solving the set 

2,2 o 4 + 4 2 2 4S LSR C L + Cr 

P +1 = 2 7 v7 , V - Q = 7 ^ , PQ = (5.4) 

{Cl - C R ) 2 [C L - Cr) 1 C L - Cr 



76 



1.0 



0.5 



t o.o 

o 
o 



-0.5 



-1.0 







\ 




- V 2 

- \ P7Q=-0.25 




\ 




\ 




\ 

\ 


5=0.38 pb 


\ 








\ 


^ P7Y=-5.0 


\ 




\ 

\ 





-1.0 -0.5 0.0 0.5 
cos2(|) L 



1.0 



Figure 20: Contours for the "measured values" of the total cross section, V 2 / Q and V 2 /y 
for m x ± = 95 GeV [m 5 = 250 GeV]. 

with cl,r = cos20£ j /j and sl,r = sin 20^, the solutions (p, q) in point (1) and point (2) 
of Fig. 21 are found for sin20£ sin20# > or < 0, and a second set is found by reversing 
the signs of the solutions pairwise. From this, one obtains for tan/5: 



cos 2(f) R > cos20£ 
where x 2 



tan f3 



p< 



-g 2 ±2 v / X 2 (p 2 + g 2 + 2-x 2 



( V T+F-v / T+7) 2 -2x 2 



tan/3 >1 (5.5) 



If the denominator is positive, there are either up to two solutions 

for tan/3 > +1 in point (1) and none in point (2), or at most one in point (1) and at 
most one in point (2). The possible solutions can be counted in an analogous way if the 
denominator is negative; the roles of point (1) and point (2) are just interchanged. The 
same counting is also valid in the second case cos20k < cos2</>£ : tan/3 — > 1/ tan/3. Thus, 
only a two-fold ambiguity is inferred from all solutions in point (1) and point (2). 

The gaugino and higgsino mass parameters are given in terms of p, q by: 




Mi 



w 



V2 



(p ± q) sin (3 — (p =F <?) cos (3 



(5.6) 



The parameters M 2 , \i are uniquely fixed if tan/3 is chosen properly in point (1) and/or 
point (2). Since tan/3 is invariant under pairwise reflection of the signs in (p,q), the 
definition M2 > can be exploited to remove this additional ambiguity. Returning to the 
"experimental values" of mass, cross section and spin correlations introduced above, the 
following parameters are extracted for point (2) [tan/3; M 2 , //] = [1.0; 53 GeV, —52 GeV] 
or [1.4; 240 GeV, 137 GeV]. Point (1) gives negative values for tan/3 so that the solution 
derived from the "experimental values" is unique in this case. 

Thus, the fundamental SUSY parameters tan /3; M 2 , /x can be derived from the observ- 

ables m ± and cos20#, cos20^ up to at most a two-fold ambiguity, by considering the 

Xi 
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production of the lightest chargino pair and using the partial information on the chargino 
polarizations. Moreover, from the energy distribution of the final particles in the decay 
of the chargino, the mass of the lightest neutralino can also be measured; this allows to 
determine the parameter M\ so that also the neutralino mass matrix can be reconstructed 
in a model-independent way 
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Figure 21: Determination of (p, q) from p 2 + q 2 , pq and p 2 — q 2 for pq > 0. 



5.3 Chargino/neutralino production at hadron colliders 
5.3.1 Theoretical cross sections 

Because they are strongly interacting, squarks and gluinos are the most copiously pro- 
duced SUSY particles at hadron colliders, and the processes pp — > gg,qq,qg have been 



extensively discussed in the literature; for a review see Ref. [^|. However, the charginos 
and neutralinos might be the lightest SUSY particles, and would be the first to be kinemat- 
ically accessible. In addition, even if the existence of supersymmetry would be established 
by the detection of squarks and gluinos at the LHC, the direct production of charginos 
and neutralinos would allow various tests of the model, and is therefore very interesting 
to investigate. 

Except for squark and gluino production, the next potentially largest cross section for 
SUSY particles at the LHC is the one for the production of the lightest chargino and the 
next-to-lightest neutralino: 

pp -> qq -> x°2 Xt ( 5 - 7 ) 

Although it is an electroweak process, the production rates are large enough to allow for 
a copious number of events at the LHC, even at a low luminosity C = 10 fb _1 . The 
process is of the Drell-Yan type, and proceeds via s-channel W boson exchange and 
i-channel squark exchange. In this section we revisit the tree-level calculation for the 
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partonic cross section of this process, and this for two reasons: (z) this is a good warming- 
up exercise using the semi-automatic program FeynMSSM that we are developing for 
SUSY computations, which we will describe briefly later on, and which we plan to use, in 
the future, to calculate the one-loop QCD and electroweak radiative corrections to this 
process; (ii) compare our independent calculation with two existing and disagreeing results 
in the literature |)7, 98] to settle the issue. By doing so we also took the opportunity 



to compare our notations [mainly based on Ref. H5[] and subsequent papers] with the 
notations of the previous references [which, of course, are also distinct from each other... 
For more details the reader is referred to Ref. JT^] . 



The differential cross-section for the subprocess has been given in both Refs. |57| and 



38]. However the only place where the result of Ref. [0 can be found in an analytical 



form is in the program ISAJET. When this expression is extracted from the code and 



compared with the one given in Ref. [9*8"], one finds a discrepancy. The formula from 
ISAJET reads 

da/dt = 1/(1Qtts 2 )[A s (U + T)/V s + A"(U -T)/V s 
+A' S (m x +m x o)s/V s + A u U/V 2 U + A t T/V' t 2 
+A st (s - Ml)T/(V s V' t ) 

+A' st m x + m x os(s - M^)/(p 8 V t ) + A su U(s - M^)/{V S V U ) 

+A' SU m x + m x os(s - M^)/(V S V U ) + A tu s m x + m x o/(V u V' t )} (5.8) 



where 



U = (m 2 x + - u)(m x o - u) , T = (m^+ - t)(m*a - t) (5.9) 



m x + and m x o are respectively the chargino and neutralino masses and u, t and s the 
Mandelstam variables. T> s and T> u are given by 

V s = {s-M^Y + M^T 2 W , V u = {u-m\ L ) , (5.10) 
K = , 1> t = (t-m 2 ii ). (5.11) 

where mg L , m Si are the left-handed down and up squark masses respectively. The factors 
A^'' ) encapsulate the information on the chargino-neutralino couplings to W bosons and 
sfermions. Eq. ( |5.8j ) agrees with the corresponding one in Ref. [^8| apart from the term A" 
which is missing in the latter. We repeated the calculation in a completely independent 
way and found full agreement with eq. ( |5.8| ) thus re-ensuring the validity of the expression 
used by ISAJET. 

This calculation has been performed with the program FeynMSSM, which is based 
on the packages FeynArts [which generates the full set of diagrams automatically] and 
FeynCalc [which calculates the amplitudes, and reduces them to the standard forms] 
developed |99| for a semi-automatic calculation of amplitudes in the Standard Model. 



FeynMSSM, incorporates the MSSM Feynman rules [following the notations of Ref. 
and thus generates the amplitudes for any SUSY process in principle. The program is 
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running under Mathematica2.0, and has been checked for processes with up to three 
particles in the final state, as well as loop diagrams [up to two loops] which is in fact the 
main purpose of the program. When FeynMSSM will be fully operational, the hope is 
that it would be very efficient and useful when dealing with higher-order processes, either 
multi-loop or multi-final state, in the MSSM. 



5.3.2 Searches at the LHC 

Let us now turn to the discussion of this process at the LHC. To suppress the huge QCD 
backgrounds, one needs to look at signatures with isolated leptons in the final state, and 
the many leptons there are the smaller will be the backgrounds. The one-lepton channel 
suffers from the large number of singly produced W bosons, while the two-lepton channel 
suffers from single Z boson and WW pair production. One should then look at the 
signal with three leptons plus missing energy. The leptons of the signals will come from 
the decays of the chargino into a W boson and the LSP, and the decay of the next-to- 
lightest neutralino into a Z boson and the LSP. The branching ratios are sizeable in a 
large area of the MSSM parameter space, despite of the small branching ratios of the 
gauge boson [especially the Z boson] decays into charged leptons. The main background 
reactions are due to W Z and ZZ production, as well as the semi-leptonic decays of heavy 
quarks and the Zbb background. There are also backgrounds from the production of other 
SUSY particles, and in particular squark and gluino production [which then decay, among 
other final states, into the studied neutralino and chargino via cascades] and associated 
chargino/neutralino production with a squark or a gluino; there are also backgrounds 
from the pair production of charginos and neutralinos. 

The signal can be distinguished from the background events, mainly due to the lower 
hadronic activity and not from the amount of missing energy; applying a set of standard 
[ATLAS] cuts ||59|| , one can obtain a 5a signal at the low luminosity LHC for gluino masses 
smaller than ~ 500 GeV; see Ref. [|100j| for details. Furthermore, the mass difference 
between the produced neutralino and the LSP, m x o — m x o, can be measured precisely 
from the invariant mass distribution the lepton pair coming from Z decays. 

In this section, we propose to use the charge asymmetry 

A = (N + — AL) / (N+ + AL) (5.12) 

where A + (A_) is the number of positively (negatively) charged leptons, as a new observ- 
able in LHC analyses in processes where charged particle final states are produced []TH]. 
What is exploited is the fact that, contrary to LEP and the Tevatron, the LHC is an 
asymmetric machine as far as the electric charge is concerned. Indeed, at the LHC we 
expect to produce more W + than W~ bosons, since the proton contains as valence quarks 
two u and one d quarks and the process ud — > W + is more favored than ud — > W~. 



It has been noticed in Ref. [T6|] that the charge asymmetry increases linearly with 
increasing mass of the formed final state. It can be therefore used for indirect measurement 
of particle masses, based only on the structure functions. This has been checked with 
PYTHIA 6.114 for the events W — > ev e or fiu^ using the structure functions CTEQ4 



SO 



75], and varying artificially the W boson mass. This exercise has been repeated for the 



process pp — > W^Z — > 2 jets + ev e and 2 jets + fxu^, which exhibits the same behavior. To 
check that this effect was not due to the particular choice of the structure functions, we 
used a different set of parton densities and the result was the same, except of course for 
the different slope of the charge asymmetry as a function of the mass. 

These results are illustrated in Fig. 22, where the charge asymmetry in percent is 
plotted against the mass of the final W or WZ-system at the LHC with y/s = 14 TeV, 
using the CTEQ4M parton densities. The two dark points at M ~ 80 and 171 GeV 
correspond to the asymmetries of the W and WZ systems, which are respectively 19% 
and 21%. The brighter circles and crosses correspond respectively to the asymmetries of 
the W and WZ systems when the effective mass is varied [in the latter case, the ratio 
Mw/Mz is kept constant]. One notices that these asymmetries are quite substantial, 
~ 15 to 35%, and can be therefore experimentally measurable. 
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Figure 22: The charge asymmetry [in %] as a function of the mass of the final W or 
WZ-system at the LHC with \fs = 14 TeV, using the CTEQ4M parton densities. 
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In Ref. [16], it has been therefore proposed to apply this method to the process 
PP — * X2X1 3/ ± + missing energy. Let us first recall that the production amplitude is 
dominated by the s-channel W boson exchange for heavy enough squarks. In this dia- 
gram, the udW + vertex depends only on the proton structure function, and is completely 
independent of the final SUSY state. However, this vertex decides of the W boson charge, 
and thus on the global charge of the chargino and the trilepton event. One can measure 
the charge asymmetry of the set of events which passes the cuts needed to extract the 
signal, and determine the sum of the masses m^o +m +. In Fig. 23, this mass was varied 
by scanning the MSSM parameter space, and the obtained asymmetry is shown. As can 
be seen, the charge asymmetry is always proportional to the mass of the final system. 

To confirm the usefulness of this method to measure the \ masses, a detailed Monte- 
Carlo simulation including the signals and backgrounds which pass the selection cuts is 
required. Although not very precise [there is a systematic bias from the used structure 
functions, and there is still a small dependence on the SUSY parameters due to the 
small interference between the s-channel W exchange and the t-channel squark exchange 
amplitudes!, the results are encouraging and the method deserves further investigations. 
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Figure 23: The charge asymmetry [in %] a function of m x o + m x + at the LHC with 

yfs = 14 TeV, using the CTEQ4M parton densities. 
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5.4 Two— and three— body sfermion decays 



The dominant decay modes of scalar fermions are the two-body decays into their partner 
fermions plus neutralinos or charginos: 

fi^fX°j , fi^f'xf (5-13) 

In general at least the decays into the lightest neutralino [which is the lightest SUSY 
particle] and chargino is kinematically available. In the case of squarks, if their masses 
are larger than the sum of the partner quark and the gluino masses, the decay 

q — > qg if rriq > m q + nig (5-14) 

largely dominates since it is a strongly interacting process. These decays are well-known; 
see for instance Ref. for a discussion. 



In the case of the third generation sfermions, and in particular for stops, the mixing 
can be so important that a large splitting between superpartners is generated; see section 
2.4. This opens the possibility of decays of sfermions into lighter sfermions and gauge 
bosons or Higgs bosons: 

fi - f 3 + W ± , % + 

- fj + Z, fj + h,A,H (5.15) 

These decays, including the mixing between sfermions have been discussed recently; see 
Ref. ||101|| for instance. 



Finally, the mass of the lightest stop could be smaller than the sum of the top and 
LSP masses, and smaller than the lightest chargino mass. In this case the only allowed 
decay mode will be into a charm quark and the lightest neutralino through loops |p.02j| 



ti^cxl (5.16) 



There are many situations, however, where these two-body decays of sfermions are 
suppressed or kinematically inaccessible. This happens for instance for the lightest stop 
and sbottom squarks in the case of large sfermion mixing, and for the scalar partners of 
the light quarks in models where the gaugino mass unification constraint is relaxed, and 
in gauge mediated SUSY breaking models. In this case, various three-body modes might 
become relevant. A list of three-body decays of sfermions which might be important 
includes: 



(z) Decays of the stops into a 6-quark, a neutralino and a If or H + boson: 

h,2 -> t*,b*, X + * -> b X °i W + or b X °i H + (5.17) 

The decay t — > bx\W + is especially important in the case of the lightest stop, when it 
is lighter than m& + m + and m t + m o. In this case, the lightest stop t\ will decay into 
a charm quark plus the lightest neutralino, eq. (|5.16| ). Because it is loop mediated, this 



S3 



decay has a partial decay width that is suppressed by four powers of the electroweak 
coupling and by the CKM angle V c b, compared to the usual two-body decay widths. In 
this case, the three-body decay into bx°W with a virtual exchange of a top quark, a 
sbottom and a chargino might become competitive. In fact, the decay is also suppressed 
by two powers of the electroweak coupling as well as by the virtuality of the exchanged 
particles; however, there are some kinematical regions where this mode dominates over 
the t\ — > Cxi decay channel. This decay mode has been also discussed in Ref . [ |103| | , where 
however only the matrix elements in terms of the four momenta of the particles involved 
have been given. 

The decay t — > bx\H + is similar to the previous one with the W boson replaced by 
the charged Higgs boson H + . In the MSSM the charged Higgs boson is always heavier 
than the W and this in principle disfavors this process kinematically compared to decay 
into Wbxi- However, the couplings of the H + boson to the third generation quarks and 
squarks can be strongly enhanced, leading to a possible compensation. For relatively 
small masses of the H ± boson, the decay branching ratio of this channel can be sizeable. 

(ii) Decays of the heavier stop into the lighter stop or sbottom and a fermion pair: 

t 2 -> Z*, h*, A* - tjf , h -> W*, H + * - bff (5.18) 

This decay occurs when the states t\ or b are rather light and there is a splitting between 
the two stops and/or between i 2 and the sbottoms, while the mass differences m^ 2 — 
or m ( - 2 — mi are still smaller than M /. .1//,. M \ or M W ,M H ±, respectively. When i 2 is 
lighter than m b + m x ± and m t + m x o, the two-body decays [except for the loop decay 

t — > cXi] are forbidden, and these three-body modes become relevant. For high values 
of tan f3 the corresponding decays of the heavier sbottoms [i.e. the stops replaced by the 
sbottoms in the decays above] might also be of some relevance. 

(iii) Decays of squarks into quarks and a stop or a sbottom via gluino exchange: 

q — ► g* —>■ q 1 1 or qbb (5.19) 

The decay q — > qbb of a scalar partner of a light quark occurs when gluinos are heavier 
than the squark [but not too much because the large virtuality of the exchanged gluino 
will strongly suppress the decay] and when there is a strong mixing in the sbottom sector 
which makes that the lightest sbottom is lighter than all other squarks. One also needs 
that the lightest chargino and neutralinos are rather heavy and/or of almost higgsino type 
to suppress their couplings to the squarks. However, in models with the unification of the 
gaugino masses at the GUT scale, the gluino and chargino/neutralino masses are related 
and for not too heavy gluinos [not to suppress the decay width by the g* virtuality] the 
decays into the other chargino and/or neutralinos are still allowed kinematically and would 
dominate since these particles would be then of the gaugino type. However, relaxing the 
assumption of gaugino mass universality, one can find a large area of the parameter space 
where the mode eq. (|5.19|) can be sizeable. 



For small values of tan f3, it is the mixing in the t sector which can be strong and 
ti possibly lighter than the other squarks; a similar situation as previously occurs then 



84 



for the decay q — > gtt despite of the kinematical complication of having a top quark in 
the decay product which requires larger decaying squark masses than in the previous case. 

(iv) Decays of squarks into quarks and a slepton-lepton pair via chargino and/or 
neutralino exchange: 

q^X°* ,X ±if ^qll (5.20) 

This decay occurs when sleptons are lighter than squarks and the squark decays into 
neutralinos and charginos are kinematically shut or suppressed by the small squark- 
quark-higgsino couplings. This happens for instance in gauge mediated supersymmetry 
breaking models where, first the gravitino which couples only weakly to squarks is the LSP 
and in addition there is a hierarchy which makes the squarks heavier than sleptons. The 
next-to-lightest particles could then be the sleptons [and in particular the stau's because 
of the mixing for large values of tan/3]. If squarks have masses smaller than the chargino 
or neutralino masses so that the two-body decays q — > x°Qi X l' are kinematically closed, 
the modes eqs. ( 5.20 ) become the dominant ones. 



However, even in the "phenomenological MSSM" discussed before, these decays are 
possibly important if the lightest neutralinos and charginos are pure higgsinos [thus sup- 
pressing the squark-quark-ino couplings and the two-body decays into quarks+inos] and 
the other charginos and neutralinos are not much heavier than the sleptons [not to be 
hurt by the virtually of these states]; the relaxation of the gaugino mass unification as- 
sumption will help in enlarging the MSSM parameter range where this situation occurs. 
The three-body decay modes of first /second generation sleptons into third generation 



sleptons+leptons have also been discussed in Ref. [|104|| in the context of gauge mediated 
breaking models. 



All these decays are discussed in detail in Ref. |T7]]. The Dalitz plot densities for the 
partial widths AY / (dxida^), where x\ and x<i are the scaled energies of two final particles, 
are given including the full dependence on the final fermion/sfermion masses and on the 
couplings [sfermion mixing is included for instance]. The two remaining integrations are 
then performed numerically to obtain the partial decay widths and the branching ratios. 
In the cases where there is only one massive particle in the final state [such as in case 
{iv), (Hi) for sbottoms and in case (ii) for both stops and sbottoms] the integration can 
be done analytically and is given. 

A detailed discussion of the relative magnitude of the decay modes (i)-(iv ) is given in 



Ref. [IT]]- in addition, the fortran code calculating the partial widths and branching ratios 
is available. It is based on the subroutines SFERMION, GAUGINO and SUSYCP [which 
calculate the masses and couplings of sfermions, charginos/neutralinos and Higgs bosons 
respectively] included in the program SUSPECT discussed in section 3.2. Therefore, an 
interface with SUSPECT is easily possible and will be done in a near future. 
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5.5 Stop and sbottom searches at LEP200 



As discussed previously, squarks of the third generation, namely the stops and the sbot- 
toms, have a special place in the SUSY spectra due to the particular Yukawa couplings 
of their partners, the top and bottom quarks; stops and sbottoms could even have masses 
accessible at LEP200. In this section we discuss the searches of these squarks at LEP200 



H|, focusing on the scenarii where R-parity is conserved and in both the unconstrained 
MSSM and the mSUGRA models. The notations and conventions are those introduced 
in Section 3. 



5.5.1 Squark production at LEP200 



At LEP200, squarks of the same mass are pair-produced through 7/Z s-channel exchange. 
One notes that the squark-squark-Z boson coupling can vanish for the special value of 



the mixing angle: 9 q = A q cos(y e^s^/Jf), which corresponds to 0.98 rad for the stop and 
1.17 rad for the sbottom. At the Born level, the production cross section is given by ||105|| : 



a 



Born 



ft 3 



4 + 



\vl + a 2 e )vl 



StM* - M z) 



M\f + r|Mf 



(5.21) 



with v e = 2s 2 y — 1/2, a e = —1/2 and v q = 2(/| cos 2 9 q — Qg sin 2 Ow)] ft is the squark 
velocity. This cross section is minimal for: 

L e + R e 



cos 6L 



T q 
is 



1 + f 1 - ^ 



(5.22) 



with L P 



'14' 



1/2 and R e = s^. The mixing angle values are equal, up to a level of 5 
%, to those corresponding to a vanishing squark-squark-Z boson coupling. 

The QCD corrections are factorizable and one obtains for pure gluon exchange and 
emission in the final state [ 106 , 107 ]: 

4a 



a 



QCD 



a 



Born 



3 IT 



f(ft) 



2ft 



l + ft fn 2 



-3 



(5.23) 



[the expression of f(ft) is given in the Schwinger approximation which reproduces the 
complete formula at the 1.5 % level]. In the high energy limit ft — > 1, f(ft) — > 3 leading 
to a correction four times higher than the one corresponding to quark production in the 
same limit. 

The electromagnetic corrections can be taken into account by convoluting the cross 
section with ||106|| : 



[x 



ft em (1 



- 1 (l + ^/3 em )-i/3 em (l + x) 



7T 



1 + a em I - - — 

Z7T 



(5.24) 



where ft em = 2a em /n (Log(s/m 2 ) — 1) and a em the QED coupling constant defined at 
zero-momentum transfer. One obtains then the total cross section: 



a tot (e + e 



qq 



L ee (x)a QCD {xs)dx 



(5.25) 
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The various cross sections: <7 Born , <7^ CD , cr^ ED and the total cross section <r tot including 
all corrections are shown in Fig. 24 at a cm. energy y/s = 172 GeV in the case of the 
lightest stop and sbottom squarks with masses of 80 GeV and as function of the squark 
mixing angles. One can see that the QCD corrections, which are always positive, increase 
when rriq — > y/s/2, reaching 60% for a 80 GeV stop mass. They are constant and equal to 
15% far from the production threshold. QED corrections are negative near the kinematical 
threshold and positive when otherwise. Close to threshold, the QED corrections tend to 
compensate the QCD corrections and the result is a total cross section equal at the level 
of 5% for 771^=80 GeV to that of the Born approximation. 
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Figure 24: Stop and sbottom production cross sections at i/i=172 GeV as a function of 
the mixing angles and for a squark mass of 80 GeV. 



5.5.2 Squark decays 



In this section we discuss the decays of the stop and sbottom squarks in the framework 
of the uMSSM and mSUGRA. Let us first concentrate on the case of the lightest stop. 
According to the discussion in the previous section on sfermion decays, there are only a 
few decay modes relevant for stop masses that are kinematically accessible at LEP200, 
i.e. m h < 100 GeV: 

- The loop induced flavor changing two body decay mode into a charm quark and the 



lightest neutralino, t\ — 
decay width is given by 



ex?, always occurs since the lightest neutralino is the LSP. The 



r(ti 



a 



\e\ 2 mz 



hi 



2 

777j 
tl 



(5.26) 



where bn depend on the neutralino parameters /i, Mi, M2, tan/3. The e parameter takes 
into account the various possible loop contributions and is estimated [ |102j| to be e ~ 
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(1 — 4) x 10~ 4 . This small value has the consequence that the stop decay time is far longer 
than the strong- interaction time scale, tqcd ~ 10~ 23 s. If this channel dominates, the 
stop hadronises first before decaying. The corresponding decay into charm plus a gluino 
is ruled out for stop masses accessible at LEP200, due the experimental bound on the 
gluino mass from Tevatron [55]|. 

- The three-body decays t\ — > bl + v and t\ — > bl + v through the exchange of an off-shell 
chargino are allowed if mi x > m„ + m;, or > m\ + m&. The corresponding widths 
can be found in ||102 |. The first of these decays is more favored by kinematics since the 
experimental limit [55] on the sneutrino mass, > 37 GeV is weaker than the one for 
charged sleptons, mj > 65-85 GeV depending on the flavor. The other possible three- 
body final states are: t\ — ► bW + Xi and t\ — > bH + Xi through exchanges of sbottom, 
chargino or top quark. Nonetheless, in the MSSM, the H + mass is larger than the W 
boson mass, and the experimental bound on the LSP |55| tends to close this channel. 

- The decay into a bottom and the lightest chargino t\ — > bxt if charginos are lighter than 
stops. However, since the experimental limit on the lightest chargino mass is rather strong 



m x + > 90 GeV [pq| , this decay occurs only for stop masses close to 100 GeV. However, 
when it occurs, this decay mode largely dominates since it occurs at the tree-level and is 
a two-body decay. 

Thus the most likely decay modes of stops at LEP200 are the t\ — > cx\ and t\ — > 
blv. To illustrate the relative magnitude of these two channels we make the following 
assumptions: {%) The lightest neutralino Xi is the LSP, (ii) the lightest stop is lighter than 
the lightest chargino, and (Hi) the lightest stop is heavier than the sleptons. Furthermore, 
we shall use the following experimental bounds on the sparticle masses from direct and 
indirect searches: a) rriy > 37 GeV, b) mj > 85, 65 et 67 GeV for selectrons, smuons and 
staus respectively and c) m- g > 150 GeV [fffil . We vary the parameters M 2 and /i in the 
range [0,1000] GeV and [-1000, 1000] GeV, respectively . 

Fig. 25 shows the accessible region in the plane (fi, M%) for which we have m x + > 
80 GeV and m x o < m~ tl = 80 GeV. The stop mixing angle has been fixed to 9 t =0 and 
6 t = 7r/2. The stop partial widths for the two decay modes t\ — > cx\ and t\ — > bru 
have been calculated for two sneutrino masses: mj>=42 GeV and 70 GeV. Contours in 
the [M2,/i] plane for the branching ratios of the decay mode t\ — > cx\ are shown. One 
can see that this decay mode is largely dominating, with a branching ratio larger than 
99%, except for small sneutrino masses [which makes the decay mode t\ — > brv more 
phase-space favored] and small stop mixing angle [which favors one of the higgsino or 
wino components of the charginos]. Indeed, for 9 t = 0, the tbxt coupling depends both 
on the wino and higgsino components. Nonetheless, if the chargino is higgsino-like, the 
Xilv coupling is proportional to the tau Yukawa coupling which is in that case very small 
[remember that we are dealing with small tan/5 values]; this means that for 9t=0, the 
decay t\ — > blv has non negligible values only for a wino-like xt that is for \fi\ ^> M 2 and 
this is precisely the region where this decay dominates. On the other hand, for 9 t = tt/2, 
the coupling stop-bottom-chargino only depends on the higgsino component of xt and 
the decay t\ — > blv is then suppressed because the chargino-lepton-sneutrino coupling is 



in that case governed by the lepton Yukawa coupling and is negligible. 
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Figure 25: Contours for the branching ratio £i — > c\\ in the (/i, M2) plane for: tan/3=2, 
mfi =80 GeV, m x + > 80 GeV, m,>=42 or 70 GeV, 6 t =Q or tt/2. 

The analysis of the stop branching ratio is simpler in mSUGRA due to the reduced 
number of free parameters. Let us compare below the two stop decays t\ — > cx\ and 
t% — > in the framework of mSUGRA for the following range of parameters: «[/ = 
1/24.3, 2 < tan/3 < 8, -2.9 < A < 2.9 [in TeV], sign(/x) = ±, < m < 1000 and 
< mi/2 < 1000. The regions in the (7710,7711/2) plane for which the following conditions 
on the sparticles masses are fulfilled are shown on Fig. 26: m x o < < 80 GeV (25a), 
m x + > 80 GeV (25b) and m x? < < m h < 80 GeV (25c). 

The intersection of the domains in Figs. 26a and 26b gives the accessible space for 
which the decay t\ — > cxI is possible with a chargino heavier than the scalar top. One 
observes that the condition m x + > 80 GeV strongly reduces the possible space; this 
is due to the universality condition in the gaugino sector which prevents a large mass 
difference between Xi an d xf ■ The intersection between the domains in Figs. 26b and 
26c is negligible, which implies that the decay t\ — > blv is highly improbable if m + > 80 

A 1 

GeV. If the electroweak symmetry breaking conditions are imposed, the parameter space 
in the (m , m^j) plane in which this decay can occur is strongly reduced as it is observed 
in Fig. 26. In conclusion, the decay ti — > cx\ has a branching ratio near unity in the 
constrained MSSM; however the corresponding accessible space in the plane (mo,™^) is 
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small, especially if a proper EWSB is required. 
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Fi gure 26: Accessible domain for the decay t\ — > cx\ 
of parameters discussed above, and with the requirements 

for a), mjj <80 GeV and m x + >80 GeV for b) and m x o < m„ < with < 80 GeV 
for c). The points correspond to cases for which no correct EWSB is required and the 
bullets to accessible domains with a correct EWSB. 



Finally, let us briefly discuss sbottom decays. The only two accessible decay modes 
from sbottoms which can be produced at LEP200, are b\ — > &X12 [the decay b\ — > bg is 
forbidden by Tevatron results]. The relative magnitudes of these two decays depend on 
the sbottom mixing angle 6^ and on fi, M 2 and tan f3 parameters [we assume here the GUT 
relation between M\ and M2]. More precisely, if M2 S> Xi an d X2 are higgsino-like 
and the couplings b\bx° do not depend on 9b and are similar; the decay b± — > 6%? then 
dominates over the decay 61 — > • ^ I/ 1 1 ^ ^2 then x^ is photino-like and x\ is zino- 
like; the coupling b\bZ is minimal for 9^ = tt/2 and in this case b\ — > bx\ dominates over 
bi — > bx2- On the other hand, the coupling b\bZ is maximal for 6^ = and is larger than 
the 6167 coupling, leading to a domination of bi — > over ^1 — > ^Xv The fact that the 
sbottom hadronises or not, which is a question of experimental interest when searching 
for sbottom squarks which are mass-degenerate with neutralinos, also depends on 9f, and 
\i, M 2 ,tan/3 but is complicated by the lack of knowledge of the precise QCD time scale. 
Studies concerning this subject can be found in [IT3]. 
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6. Experimental bounds on SUSY Particle Masses 



6.1. Introduction 

A wide range of searches for supersymmetric particles are performed at present colliders 
and no deviation with respect to the Standard Model predictions was observed yet, unfor- 
tunately. Therefore new limits on the masses of these particles, assuming different models, 
are set by the various experiments. In this section we focus on the experimental results ob- 
tained both at LEP and the Tevatron, interpreted in the constrained MSSM or mSUGRA 
framework [but for LEP2, without the constraint of correct radiative electroweak symme- 
try breaking and without assuming any mass unification for Higgs bosons with the other 
scalar particles] assuming that the LSP is the lightest neutralino Xi an d that R-parity 
is conserved. In this model, as discussed in section 2, all SUSY particle masses, their 
couplings and their production cross sections and decay widths are predicted in terms of 
only five free independent parameters: m^, m , /i, tan/3 and A; the renormalization 
group equations are used to determine the parameters at low energies. Instead of mi/2, 
the LEP experiments usually derive their results as a function of the parameter M2, the 
SUSY breaking mass term associated with the SU(2)l gauge group (mu2 ~ M2/O.8I). 
At the Tevatron, the previous assumptions are made in the interpretation of the results, 
with the additional constraint of the radiative electroweak symmetry breaking. In this 
case, is no more a free parameter and is determined when m , tan/?, mu^ and the Z 
boson mass Mz are fixed; the sign of \x remains free; see section 2. 

The searches for supersymmetric particles at LEP2 concern sleptons, stops, sbottoms, 
charginos and neutralinos as detailed in ||108|| . These various SUSY particles decay to SM 



particles and two LSPs; therefore, SUSY signatures consist of some combination of jets 
or/and leptons and missing energy since the LSP escapes detection. The signal topology 
and the background conditions are in practice affected by the SUSY particle and the 
LSP mass difference (AM = msusY — m x°) which controls the visible energy. In the low 
AM(= 5 — 10 GeV) range, the expected topologies for the signals are characterized by a 
low multiplicity and a low visible energy and the background is dominated by two-photon 
interactions. For large AM(= 50 — 60 GeV) values, the signal signatures are very similar 
to those of PU-pair production. Since all the background sources are due to well calculable 
processes with reasonable production cross sections compared to the signal, most of the 
decay channels are studied at LEP 2. 

Another important key domain concerns the searches for the MSSM Higgs bosons. 
Presently, only the lighter neutral Higgs bosons h and A can be discovered at LEP2 [ |10S 



since the CP-even Higgs boson H and the charged Higgs bosons H ± are expected to be 
too heavy. However searches for charged Higgs bosons are also performed at LEP2 in the 
framework of a non-SUSY two-Higgs doublet model. 

After collecting ~ 150 pb _1 at LEP1 (1989 -> 1995), each experiment (ALEPH, 
DELPHI, L3, OPAL) has accumulated data at LEP2: ~ 5.5 pb _1 at = 133 GeV 
(1995), ~ 10 pb" 1 at y/S = 161 GeV (1996), ~ 10 pb" 1 at yfs = 172 GeV (1996) and 
~ 55 pb" 1 at y/s = 183 GeV (1997). This year, an amount of 150 pb" 1 at y/s = 189 
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GeV (October 1998) is already recorded by each experiment, but the results including all 
the statistics for 1998 are not yet available; a fraction of the overall integrated luminosity 
is only used corresponding to the statistics collected just before the summer conference 
time (-> July 1998 : 30 -40 pb" 1 ). 

At the Tevatron the main sources for SUSY are squarks and gluinos, abundantly 
produced due to the color factors and the strong coupling constant. Squarks or gluinos 
are produced in pair, and decay directly or via cascades to at least two LSP's. The classical 
searches rely on large missing transverse energy (E£) caused by the escaping LSPs. In 
addition, charginos and neutralinos are searched for via their leptonic decay channels by 
the two Tevatron experiments CDF and DO. Finally, searches for the MSSM charged 
Higgs boson are performed at the Tevatron, and bounds on the charged supersymmetric 
Higgs boson mass have been set by both experiments. Each experiment has collected an 
integrated luminosity of about 110 pb _1 at ^/s = 1.8 TeV. 

Since no evidence for production of supersymmetric particles has been found, exper- 
imental limits on their production cross sections and masses has been derived. In what 
follows, we will pay attention to all the decay channels used to extract the limits both 
at LEP2 and Tevatron. The most recent results have been used, including preliminary 
results reported at the last summer conference in Vancouver and the last LEPC meeting. 

6.2 The scalar particle sector 
6.2.1 The Higgs bosons 

Due to the expected large mass of the heavier CP even Higgs boson H, only searches for 
the neutral h and A bosons will be reported in what follows. These neutral Higgs bosons 
can be produced at LEP via two complementary processes: 

• The Standard Model-like Higgs-strahlung process, e + e~ —>■ hZ; the cross section 
for this process is equal to its SM analogue, reduced by a factor sin 2 (a — (3) which 
means that this process occurs mainly at low tan/3 values. 

• The associated pair production process, e + e~ — > hA, with a cross section propor- 
tional to cos 2 (/3 — a); this process is dominant at large tan/? values. 

In the first channel, searches are similar to those performed for the SM Higgs boson. 
For tan/3 >1 the main decay modes of the h and A bosons are into bb, and to a lesser extent 
t + t~ . Most of the experimental analyses required therefore at least two b quarks in the 
final states, as detailed in Tab. |5|. Clearly, 6-tagging plays a crucial role in all the Higgs 
searches at LEP; it allows to reach very high sensitivities by rejecting most of the WW 
background. Peculiar decays such as h — > AA or h — > XiXi have been also considered; the 
sensitivities achieved in the latter case are better compared to those obtained with the 
standard decay channels at LEP1 but worse at LEP2. Typical efficiencies and background 
expectations are listed in Tab. where observation and expectation from SM processes 
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1111-1116 , for two "benchmark" sets of the 



agree very well. 

Therefore, limits have been extracted 
MSSM parameters where M susy , representing the stop mass mean value in Figs. [26 
fixed to 1 TeV and m t = 175 GeV; Ma varies up to 2 TeV and tan/? from 0.5 to 50, and 
either minimal or no stop mixing [A = TeV, /i = —0.1 TeV] or maximal squark mixing 
[A = V6 TeV, n = —0.1 TeV] is assumed. 
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hA — > qqr~T + , r~T + qq 
hA^ AAA 
AAA^ bbbbbb 


Zh— > uuqq 
Zh-^ El + qq 
Zh^ t~t + qq 
Zh— > qq t~t + 
Zh.-> qqh(^ Xi Xi) 


136 


hA— > bbbb 
hA^ AAA 
AAA^ bbbbbb 




161,172,183 


hA^ bbbb 

hA^ bbr~r + , r~r+ bb 
AAA-^ bbbbbb 


Zh^ i/z/(h^all) 
Zh-> ri+(h-»all) 
Zh^ r-r+(h^all) 
Zh^ (h -> all) t-t + 
Zh— ■> qqbb 

Zh— > qq or l~l + h(-^ Xi Xi) 



Table 5: Search channels for the neutral Higgs bosons at LEP. 



The limits obtained by each experiments within these assumptions are listed in Tab. [7|. 
The combined limit, obtained as described in Ref. ||110|| leads to a gain of about 3-4 GeV 
with respect to the individual experiments. Four methods for the combination are used, 
they all agrees within a spread of ± 1.9 GeV; conservatively the lower limit is quoted in 
Tab. ^ (method C of Ref. [ 11C|| ), while the higher limits (with method D) are given in 
Figs. ^6[ These limits are valid for tan/3 greater than 0.8, irrespectively of the mixing value 
or of the method used for the combination. For the no mixing scenario, the range of tan/3 
between 0.8 and 2.1 is excluded independently of the method. It has to be noticed that 
the limits have been improved, for each LEP experiment, by almost 4 GeV/c 2 with the 
recent data collected at 189 GeV 1 1 1 7| | , reaching therefore the level of the combined limit 
at 183 GeV. It has been pointed out recently in Refs. [|116j| , ||ll4|1 , ||112|| , that selected sets of 
the MSSM parameters lead to less stringent exclusions, either because of an unexpectedly 
small bremsstrahlung cross section or because of a reduced decay branching ratio into bb. 
The analysis of Ref. [|112| shows however that the probability of occurrence of such sets is 
extremely small, typically at the 10~ 3 level. The theoretical status of these pathological 
parameter sets is currently under investigation. Meanwhile, the limits obtained in the 
benchmark cases can be regarded as sufficiently robust for practical purposes. 
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0.27 
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Table 6: Typical efficiencies and expected/observed number of events for the neutral 
Higgs bosons searches performed by the ALEPH experiment |111|| at ^=183 GeV. 





M hiA lower limit (GeV/c 2 ) 




Up tO y/s = 


183 GeV 




Range of 


lower 




Individual Limit 


LEP combined 




(ADLO) 




h (Obs.) 


70.7-74.4 


78.8 


(Exp.) 


67.4-70.3 


76.3 


A (Obs.) 


71.0-76.1 


79.1 


(Exp.) 


68.4-72.0 


76.3 



Table 7: Individual and LEP combined mass limits for the neutral h and A bosons ||110|| . 



At the Tevatron, the most widely studied mechanism for producing neutral Higgs 
bosons is the associated production of the CP-even h boson with a W or Z boson. 
In addition, one has to require the gauge boson to decay leptonically; this reduces the 
observable rate by a factor greater than 75 %. Therefore a collected luminosity of several 
fb _1 is needed to reach a sensitivity competitive with those achieved at LEP. Recently, 
in Ref. ||119|| , it has been made use of the fact that the Abb coupling is proportional 



to tan /3, i.e. the cross section production grows as tan 2 (3. Analyzing the events from 
CDF containing r _ r + + 2 jets in the final state [used to set limits on third generation 
leptoquarks], the authors derived new bounds in the plane (Ma, tan/3); in particular if the 
mass of the CP-odd Higgs boson is just beyond the present limit, the tan/3 region above 
80 is excluded. 

Let us now turn to charged Higgs bosons. The H ± bosons can be produced in e + e~ 
interactions via the process e + e~ — > 7*, Z* — > H + H~ and are expected to decay mainly 
into the heaviest lepton kinematically allowed and its associated neutrino, or into the 
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Figure 26: LEP combined exclusion domains (method D) in the plane (Ma, tan/3) (top) or 
(Mh, tan/?) (bottom) by the LEP direct searches for e + e~ — > hZ and hA at center-of-mass 
energies up to 183 GeV assuming either a maximal mixing or no mixing |110|| and [111 8 
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heaviest kinematically allowed quark pair whose decay widths are not Cabibbo suppressed; 
the relative branching ratios are model dependent. Therefore searches are performed in 
the three possible decay modes: e + e _ — > H + H~ — > t~v t t + v t , t + v t sc and esse. For 
each experiment |[L20|| - [|123j| and each decay channel, the number of selected events in the 
data is consistent, with the number of expected events from SM processes, as can be seen 
in Tab. || Compared to the neutral Higgs boson searches, the sensitivities achieved are 
lower; this is due to the high WW background contamination. The 95 % CL. individual 



M H ±=60 GeV/c 2 


r + u T r + u T 


T + U T SC 


esse 


Efficiency 


24 % 


42 % 


40 % 


Expected events 


9.2 


30.1 


99.4 


Observed events 


6 


28 


93 



Table 8: Typical efficiencies, number of expected and observed events obtained by the L3 
experiment in the charged Higgs boson searches at cm. energy of 183 GeV [|122j| . 





M H ± Lower limit (GeV/c 2 ) 




up to v/s = 183 GeV 




Range of 


Lower 




Individual Limit (ADLO) 


LEP combined 


H± (Obs.) 


56.6-59.0 


68.0 


(Exp.) 


56.0-62.0 


69.0 



Table 9: Individual and LEP combined mass limits for the charged Higgs bosons. 



limits on the charged Higgs boson mass is shown as a function of the branching ratio 
Br(iJ ± — > t^Vt) in Fig. |27|. Combining all LEP results, a lower limit on the charged 
Higgs boson mass of 68 GeV is established at the 95 % C.L., assuming just that the sum 
BR(/J + — > t + u t ) + Br(H + — > cs) is equal to one. It represents a gain of almost 10 GeV 
with respect to the individual limits, Tab. |9|. 

In the MSSM, the charged Higgs bosons are expected to be heavier than the W boson, 
and can be searched for in the decay products of the top quarks that are produced at the 
Tevatron. If the charged Higgs boson is light enough, the two possible decay channels for 
the top are either t — > W + b or t — > H + b; as an example for M H ± = 100 GeV, the decay 
t — > H + b becomes dominant for both large ( > 30) and low ( < 1) tan/3 values. Moreover 
for low tan/3 values H + decays dominantly into cs leading to final states containing six 
jets without isolated leptons. The measurement of the top branching ratio into W + b, 
tagged with leptonic events, allows to set an upper limit on the branching ratio of the 
decay t — > H + b, assuming an expected tt production cross section not larger than 5.6 
pb; CDF [ |124j| and DO [|125|| set a limit of about 25 % which can be turned into excluded 



regions in the (M#±, tan/3) plane as depicted in Fig. |28|. For instance, for M#± = 100 
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Figure 27: LEP combined exclusion domains in the plane (M#±, i?r(H ± — ► rz/ T )) by the 
LEP direct searches for e + e~ — > H + H~ at center-of-mass energies up to 183 GeV ||110|| . 
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Figure 28: DO exclusion domains ||125|| in the plane (tan/3, M^) from top decay. 
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GeV, all tan/5 values below 1.1 and above 40 are excluded by Tevatron searches [|118|| . 
These limits may be also turned into excluded region in the plane (M h , tan/?) as depicted 
in bottom Fig. |2§ JTT8J and []TT(J. 



6.2.2. Scalar leptons 

LEP experiments excluded scalar neutrinos with masses up to 43 GeV, from the compar- 
ison of the measured Z widths with the SM expectation [127|. 



In e + e~ collisions, the production of scalar muons /t, scalar taus f proceeds via 7 or 
Z exchange in the s-channel only, whereas scalar electrons e can also be produced by 
exchanging neutralinos in the i-channel. At the Tevatron, charged scalar leptons can 
be pair produced or produced in association with a scalar neutrino [for the left handed 
charged sleptons only] via the Drell-Yan mechanism. The detection of scalar leptons at 
the Tevatron is difficult due to the high background sources, the low production cross 
sections and the possibly important cascade decay branching ratios. Presently, charged 
scalar leptons can be searched for only at LEP2. 

The scalar leptons dominantly decay into their SM partners and the lightest neutralino 
Xi- If the scalar lepton is not the next-to-light SUSY particle (NLSP), cascade decays 
into X2 or i n ^° the lightest chargino xf [f° r only for a pure gaugino-type chargino] 
may be possible. The topologies arising from scalar lepton production are then usually, 
acoplanar leptons plus missing energy. As already explained in the introduction, the 
sensitivities for these searches depend strongly on the mass difference between the scalar 
lepton and the x\i AM. Typical efficiencies and expected events from SM processes for 



different AM ranges are given in Tab. |T0 
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Table 10: Scalar electron, scalar muon and scalar tau efficiencies (e) and the number of 
events expected from SM processes (iV ex p). Results are obtained at yfs = 183 GeV by the 
ALEPH experiment ||128|| as a function of AM for m~ l± = 75 GeV. 



No excess with respect to the number of events expected from SM processes was 
observed at LEP up to -\/s=183 GeV [ |128| |-[ |i"3"2]j . Since in general, the cross section for 
the pair production of Ir is smaller than for II, limits on masses are given by default taking 
into account IjJr production only. The expected background from W-pair production is 
subtracted for all LEP analyses. All limits are derived in the plane (m x o,m / - jj ). These 
limits depend only slightly on [i or tan/3 for jln and tr, since these parameters may just 



98 
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79-83 


85 
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75 
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75 


(Exp.) 


45-56 


67 



Table 11: Individual and LEP combined mass limits for eR, p,R, tr up to y/s = 183 GeV 
for AM greater than 20 GeV, and assuming tan/? = 2 and fi = —200 GeV for e R . 

change their couplings to x\ [relevant only for low masses of Xi i n t ne region covered 
anyway by the xt searches] but do not affect their cross section for pair-production. 
For £r, the production cross section in the Xi t-channel exchange depends on these 
parameters. The combined exclusion contours of all LEP analyses ||126|| up to 183 GeV 



are shown in Figs. |29|. The individual limits are derived for different tan/5 values 1.4 — > 
2, and for // = -200 GeV. 

Since we are far from the kinematical limit, more luminosity helps and a gain ranging 
from 2 to 10 GeV, depending on the scalar lepton flavor, is obtained in the combination 
as shown in Tab. O. We should mention also that the individual limits obtained recently 
at -\/i=189 GeV with an integrated luminosity of 30-40 pb _1 for each experiment reach 
the level of the LEP combined [up to 183 GeV] limit 

Assuming a common scalar lepton mass at the GUT scale, the relation between the 
masses of right- and left-handed sleptons can be used to combine results of the searches 
for acoplanar leptons (muons, electrons together) [|128|1 coming from the e^ei, process 
when the e~R and the Xi are mass degenerate (AM below 3 GeV). The result is shown in 
Fig. ^9), for tan/5 = 2 and \i = 100 GeV; this value of \i minimizes the product of cross 
section times branching fraction for the process e + e~ — > e~Re~L with cl — > exl for vanishing 
AM. In this case, a scalar lepton mass limit of 65 GeV/c 2 is set independently of AM; 
this limit holds for higher values of tan/3. 



6.2.3 Scalar quarks 

The squarks of the two first generations should be abundantly produced at the Tevatron 
which places already limits far above the kinematical reach of LEP2. Due to a large 
Yukawa coupling, the mass of one of the top squarks can be significantly smaller compared 
to those of the other scalar quarks; the large mixing between the left and right stops, 
proportional to m t , leads to a large splitting of the two mass eigenstates. The lighter 
stop t 1 could be within the discovery range of LEP2. In addition, for large tan/5 values 
( > 30), the mixing angle in the sbottom sector can also be large and the lighter sbottom 
b± could also be within the reach of LEP2. 
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Figure 29: Combined exclusion contours for scalar electrons, muons and taus ||126 



as a 



function of the LSP mass and the limit for sleptons obtained by the ALEPH experiment 



128|| when combining all selectron and smuon searches up to 183 GeV (bottom-right). 
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Stop and sbottom production at LEP proceed via Z/7 exchange in the s-channel; the 
production cross sections depend on the squark mass and the squark mixing angle cos 9 q 
and at co$Q q ~ 0.57(0.39) the stop (sbottom) decouples from the Z-boson and the cross 
section is minimal; see section 5.5. The dominant decay modes of the t± are expected to 
be ti — > CXi or t\ — > b&l + ] both of these decay modes have been searched for at LEP. The 
dominant decay mode of &i is expected to be bi — > bxl ■ 

Under the assumption of R-parity conservation, the Xi an d the v are invisible in the 
detector; thus stop or sbottom pair events are characterized by two acoplanar jets or two 
acoplanar jets plus two leptons, with missing energy. When the t\ — > cx\ decay mode is 
dominant, the corresponding decay width is small enough for the top squark to hadronize 
into a colorless "stop hadron" before decaying; this feature has been implemented by 
the LEP experiments in their Monte Carlos. No excess of events was observed by any 
of the four LEP experiments [ |133j |-[ |l~3"6|| , in the searches for stop and sbottom. Typical 
sensitivities obtained by the OPAL [|136|| experiment in the various channels, are listed in 



Tab. 12 
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Table 12: ti and bi efficiencies (e) and the number of events expected from the SM (iV ex p) 
obtained at = 183 GeV by OPAL ||136| as a function of AM for nig = 80 GeV. 



All LEP results have been combined [[126|| and are summarized in Tab. [13] for AM 
greater than 15 GeV/c 2 . Since a slight lack of events appears for sbottom searches, con- 
servatively, the combination in this channel is done without any background subtraction. 
The LEP combined excluded regions in the plane (m x o,mq) are shown in Fig. |3D|. 

At the Tevatron, top squarks can be directly pair produced or can be produced in 
the top quark decay mode t — > Xiii- The stop is searched for in the t\ — > CXi decay 
channel but also in the mode ti — > bxt when kinematically allowed. Recently, the CDF 
[L37|| experiment has reported a higher limit on the ti mass of about 122 GeV/c 2 for m x o 
lower than 50 GeV/c 2 , the sensitivities in the low AM ranges are weak, due to the huge 
irreducible QCD background. The complementary domains excluded both at LEP2 and 
the Tevatron are shown in Fig. 

At hadron colliders, one expects the strongly interacting SUSY particles, gluinos and 
squarks, to be produced with the largest cross sections. Therefore the search for qq, 
qg and gg final states has been performed at the Tevatron by both the CDF [138] and 
DO [139] collaborations. Depending on the mass hierarchy between squarks and gluinos, 
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Figure 30: LEP2 (top) [|1 26|| and CDF (bottom) |137| exclusion domains in the plane 
(m^m^) from stop decay t x — > c\i searches. 

102 











, lower limits 


(GeV/c 2 ) 










up to ^ = 183GeV 












Individual Limit (ADLU) 


LEP combined 










or AO only: * 


obs 


[exp] 


J- 

h - 


-* c Xi, 


COS# t 


= 1 


80-85 


86 ( 


35) 


ii - 


■+ CXi, 


COS# t 


= 0.57 


72-81 


83 ( 


30) 


ti - 


-> bul, 


COS# t 


= 1* 


84-85 


87 ( 


36) 


h - 


-> 6z>Z, 


COs6> t 


= 0.57* 


80-82 


85 ( 


34) 


h- 


-to?, 


COS^b 


= 1 


78-84 


86 ( 


33) 


h- 




COS^b 


= 0.39 


45-68 


75 (60) 



Table 13: Individual and LEP combined mass limits for stop and sbottom up to y/s = 183 
GeV for AM > 15 GeV, assuming maximal or minimal production cross sections. 



the subsequent decays will be q — > qg or g — > qq; then decays of squarks and gluinos 
proceed according to q — > qxl and g — > qqXi [the so called direct decays] while for 
sufficiently massive squarks and gluinos, the decays proceed through charginos Xi,2 or 
heavier neutralinos X234? which in turn decay to quarks, leptons or neutrinos and one 
LSP [the so called cascade decays]. In the direct decays scenario, this leads to final states 
containing a significant number of jets and missing transverse energy carried away by the 
LSPs; it has to be noticed also that gluino decays yield on average one more jet than 
squarks decays. The cascade decays result in the production of a larger number of jets, 
with a reduced but still significant E/. and occasional production of leptons. The recent 
dilepton based SUSY search by CDF and DO are sensitive to cascade decays. 

No excess of data has be found at the two Tevatron experiments ||138| |- p39| ; as an 
illustration, the CDF experiment for an integrated luminosity of 19 pb _1 , has selected 
18 events in the three jets sample and 6 events in the four jets sample, these numbers 
are consistent with estimates of SM processes, respectively 33ljg (stat) (syst) and 
8l| (stat) t 4 4 (syst). The main background in this analysis comes from W/Z + jets 
processes, top quark production and QCD multi-jet production in which the E/_ originates 
from mismeasurement. 

These searches yield lower limits on gluino and squark masses and the domains ex- 
cluded by both experiments in the (mg,m^) plane are shown in Figs. |31|, where it is 
assumed that all but the top squark masses are degenerate. Assuming that squarks are 
heavier than the gluino, a lower limit on the common squark mass equal to 216 and 260 
GeV/c 2 is set by the CDF and DO collaborations, respectively. These limits do not depend 
on tan/3 but hold only for large \i values \\i > 200 GeV or /x < — 100 GeV]. A tighter 
squark mass limit nig > 267 GeV/c 2 , using the dilepton channel, is set by DO; in this 
case however the result is more model dependent [charginos and heavy neutralinos have 
to decay leptonically, etc.]. 
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Figure 31: DO (top) ||139|| and CDF (bottom) ||138|| exclusion domains in the plane (m 
from the multi-jet plus E/. or dileptons plus jets plus E/_ searches at the Tevatron. 
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6.3. The gaugino sector 
6.3.1 Gluinos 

The gluino mass bounds derived from the qg searches previously described, are respectively 
173 and 187 GeV/c 2 for the CDF and DO collaborations, for any squark mass. When one 
imposes an equal mass for squarks and gluinos, the limits become 216 and 260 GeV/c 2 
respectively, but are more model dependent. 

In 1996, two small windows for very light gluinos [below 1.5 GeV/c 2 and between 3 
and 5 GeV/c 2 ] remained unexcluded by any of the existing searches. Such light gluinos 
would modify the decay pattern of the SUSY particles considered; since with gluinos 
lighter than 5 GeV/c 2 , E/_ becomes small and the usual jets + E/[ signature is no more 
effective. However, resonant production of squarks would have a large cross section and, 
if squarks are not very heavy, broad peaks in the dijet mass distributions are expected. 
Comparison of the observed spectrum with theoretical estimates rules out light gluinos 
if squarks are lighter than about 600 GeV/c 2 ||127|| . Such light gluinos would also affect 



the usual phenomenology of QCD at LEP, by changing the topology of four jet events 
via g* — > gg splitting; they would also contribute to the running of a s as three additional 
flavors in leading order, up to mass effects. A full four-jet analysis carried out at LEP1 
by ALEPH [ |140H excludes gluinos with masses smaller than 6.3 GeV/c 2 . 



6.3.2. Charginos and neutralinos 

We now turn to the search for charginos and neutralinos at LEP2. At e + e~ colliders, 
charginos XiXi (neutralinos XiX^ with i,j — 1, . . .,4 ordered by their masses) are pair- 
produced via s-channel j/Z (Z) boson and t-channel sneutrino v (selectron e) exchange. 
When the masses of the sfermions are very large, the xt (Xj* with j > 2) states decay 
via an exchange of a virtual W (Z) boson as follows: xf ~~ > XiW* — > Xiff (Xj ~ * 
XkZ* ~> Xkff with k < j). If the slepton and sneutrino masses are comparable to Mw 
(Mz) the charginos (next-to-lightest neutralinos) decay via virtual slepton or sneutrino 
exchange and the leptonic branching fraction is enhanced. Finally for slepton or sneutrino 
lighter than the chargino (next to lightest neutralino), the decay modes xt ~~ > or l^v 
(x° — > or Du) become dominant. The radiative decays x°j ~^ Xkl are a l so possible 
via higher-order diagrams. For all LEP2 studies, the gluino is supposed to be heavy, 
otherwise the chargino decay patterns would be different, as reminded in the previous 
section. 



The charginos and heavy neutralinos are searched for at LEP2 [ 141 ] | 146|| , in a first 



time within the large m scenario, i.e large scalar fermion masses. In this case, charginos 
and neutralinos decay into the LSP accompanied by virtual W and Z bosons. The three 
possible final state topologies arising from the W and Z decays are: hadrons plus missing 
energy, acoplanar lepton pairs plus missing energy and mixed final state (hadrons plus 
leptons) plus missing energy. As explained before, signal topologies and the associated 
background sources depend on AM. Therefore, for each topology, selections were opti- 
mized for four different AM ranges: the very low (3-5 GeV), the low (5-10 GeV), the 
medium (20-40 GeV) and the large (> 50 GeV). In addition, the DELPHI experiment 
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searched for topologies arising in the very small AM regime, with decays in flight and 
ISR tags ||144|| . Altogether a total of 27 events was observed by the four experiments 
with 25±5 expected from standard processes with efficiencies ranging from 5% to 65% 
depending on the analysis and the AM ranges as illustrated in Tab. [14. 
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Table 14: Typical efficiencies, number of expected and observed events obtained by the 



DELPHI experiment in the chargino searches at a cm. energy of 183 GeV [|143j . 



Depending on the neutralino-chargino field content, one distinguishes the following 
cases for the determination of the lower limits on neutralino-chargino masses: 

- Higgsino-like x° an d Xi {M 2 3> |/i|): in this case, the production cross sections do 
not depend on the scalar lepton masses, and AM is low and decreases with increasing M 2 . 
Consequently, the limits on the masses of the next-to-lightest neutralino and the lightest 
chargino decrease with M 2 . The present LEP2 combined ||126|| limit on the chargino mass 
is shown in Fig. ^ (top left) turning into a limit on the chargino mass of 63 GeV/c 2 for 
AM > 3 GeV/c 2 . The DELPHI experiment has excluded regions for lower AM ranges 
as depicted in Fig. |32|(bottom) . Masses lower than 80 GeV/c 2 are excluded by all LEP 
experiments for the next-to-lightest neutralino x\ assuming M 2 less than 1500 GeV/c 2 . 
This limit holds only for higgsino-like x 2 > since the coupling to the Z-boson vanishes for 
gaugino-like neutralinos. 

- Gaugino-like chargino (\fi\ ^> M 2 ): the cross section depends strongly on the scalar 
neutrino mass. For 50 < vrty < 80 GeV/c 2 the destructive interference term reduces the 
cross section by one order of magnitude compared to what is expected for rriy > 300 
GeV/c 2 . When the two body decay xf ~^ is dominant, the relevant AM becomes 
equal to m + —mp. Therefore, the limit on the chargino mass will depend clearly on the my 
value. For large niy values (> 300 GeV/c 2 ), the situation is the ideal one: we benefit there 
from the best detection sensitivity and the highest cross section production. This is the 
reason why the kinematical limit is reached with only a few pb _1 . Up to now, including 
the recent 189 GeV data, LEP experiments exclude a chargino mass below 94.3 GeV/c 2 , 
irrespective of tan/?. Extending the study to all values where the chargino leptonic 
decay is enhanced for low slepton masses, dedicated analyses then were developed by the 
LEP experiments. However, it appears that the chargino cannot be directly detected 
when it is degenerate in mass with the sneutrino since the soft leptons produced in the 
final state escape detection. In this case, the LEP1 limit of 45 GeV/c 2 remains, that 
is deduced from the Z-boson total decay width measurement which is insensitive to the 
chargino decay pattern. Recently, this limit has been improved in two independent ways: 
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Figure 32: Top left: LEP combined [|126|| lower limit on the chargino mass for different 
AM values at cm. energies up to 183 GeV. Top right: lower limit on the chargino mass 
for different field content, as a function of tan/?, by the L3 experiment [ 145|| . Bottom: 



excluded domains in the plane 



AM) in the very low AM range obtained by DELPHI 



144j| at a cm. energy up to 183 GeV. 
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(z) the ALEPH experiment, excludes m x ± below 51 GeV/c 2 by measuring the invisible 
VT-boson width, within some model assumptions [|142|| ; and (ii) the L3 experiment ||145 | 



exploits the complementarity of the scalar lepton searches, excluding chargino masses 
below 57.2 GeV/c 2 irrespective of tan/3 and m , as depicted in Fig. [32] (top right). 

At the Tevatron, charginos and neutralinos may be produced directly in the anni- 
hilation process qq — > X2xt or hi the decays of heavier squarks as mentioned before. 
They are searched for in the leptonic decays [mainly assuming the decays xt ~^ ^ ±h> x\ 
and X2 ~^ l + l~Xi\'> see section 5.3. The requirement of three leptons in the final state 
reduces backgrounds to a very small level but is efficient for the signal only in peculiar 
cases. Therefore, charginos and neutralinos may be discovered at the Tevatron above the 
kinematical reach of LEP2, but no robust limits can be derived and the results reported 
today are not competitive with the LEP bounds. 



6.3.3 Neutralino LSP 

Finally, let us discuss the neutralino-LSP mass limits. Since the neutralino Xi is the 
LSP and escapes detection, direct searches for the LSP neutralino cannot be performed 
at LEP2. However, indirect limits have been derived from the constraints on chargino- 
neutralino and slepton searches. As detailed previously, the limits depend on the slepton 
masses, in particular on m^. 

- Large sneutrino mass: For high tan/3 values, the chargino search provides the most 
stringent bound on m x o. For tan/3 values below 2, the neutralino searches including the 

processes e + e" — > X1X3, XiX4> X2X3 but also the radiative decays x\ ~~ > Xi7> exclude 
additional regions of the parameter space (mixed region). The LEP1 exclusion limits still 
play some role, due to the fact that the neutralino search is limited in the region near the 
point where the Xi mass limit is set by the value of the coupling to the Z-boson rather 
than by kinematics. The lower limit obtained for m x o as a function of tan/3 and for heavy 
sleptons is displayed in Fig. ^ for LEP center-of-mass energies increasing from 91.5 up 
to 183 GeV. A summary of all limits [|141|| -[ |l46|j obtained both up to 183 GeV and more 
recently to 189 GeV is listed in Tab. [15]. 

- Any sneutrino mass: For lower slepton masses, as already discussed, the constraints 
in the chargino-neutralino sectors are weaker and the limits on m x o therefore degradate. 
Constrains on m x o have been obtained by all LEP experiments with a systematic scan of 
the MSSM parameter space, as a function of tan/? for all mo values, as depicted in the 
bottom of Fig. |33|. For low tan/3 values (< 1.6) the minimum allowed value for m x o is found 
in the parameter space where the production cross section for charginos is minimal [m^ ~ 
90 GeV/c 2 in the mixed region fi ~ —70 GeV/c 2 ] ||143jl - Hl46|l and the heavy neutralino X3. 



are mostly invisible while for higher tan/3 values (> 2) the lower Xi mass limit is found 
in the parameter space region where the chargino and sneutrino mass difference is small 
and where the next-to-lightest neutralinos X2 decay invisibly. The exclusion limit came 
only from slepton searches; it is localized in the deep gaugino-like region for charginos 
(/i < —400 GeV/c 2 ) [ |141| |. The limits obtained by each experiment are listed in Tab. |TB 



as well as the MSSM parameter ranges used in the scan. The most recent limit on the e 
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Figure 33: Top: ALEPH [141] lower limit on the LSP mass as a function of tan/3, for 
large slepton masses, at cm. energies up to 183 GeV. Bottom: Lower limit on the LSP 



mass as a function of tan/3, for any slepton mass obtained by L3 ||145|| at cm. energies 
up to 183 GeV. 
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m x o > (G 
up to 183 Gev 


eV/c 2 ) 

up to 189 GeV 


Large m , tan/3 > 1 


~ 30.5 (ADLO) 


~ 32.5 (AO) 


Any mo, tan/? >1 
A n up to -2000 GeV/c 2 
D fi up to -200 GeV/c 2 
L ^ up to -500 GeV/c 2 
fi up to -500 GeV/c 2 


26* 
23.4 
25.9 
25.4 


28 



Table 15: Indirect individual mass limits for the LSP Xi U P to y/s = 183, 189 GeV for 
large mo or any mo- In the latter case, the minimum value (*) for ALEPH experiment is 
found in the deep-gaugino region (up to 183 GeV data) while other experiments found 
the minimum in the mixed region. 



mass obtained with 189 GeV data by ALEPH [|141|| allows to set the best absolute limit 
on the LSP mass of 28 GeV/c 2 , which can be turned into an absolute limit on Mi- 

6.4 Summary and Bounds on the MSSM parameters 

The summary of the bounds on the SUSY particles obtained at LEP2 and the Tevatron 



are listed in Tab. [L6[ We also indicate the parameter ranges for which these limits hold. 

The limits can be turned into bounds on the cMSSM parameters; as an example, at the 
Tevatron, the excluded region in the (m^,mg) plane is translated in the plane (mo, mi/2) 
for different tan/3 values, see Fig. |4] (top). Similarly, LEP experiments have derived from 
chargino/neutralino searches excluded contours in the plane (fi, M 2 ) for different mo, tan/3 



values as illustrated on Fig. [35|; in particular the clear complementarity of the searches for 
neutralinos and charginos is visible for the low tan/? regime. Moreover, the combination 
of the gaugino and the slepton searches allowed to exclude regions of parameters in the 
plane (mo,mi/ 2 ) and to derive an absolute limit on M 2 of 47 GeV/c 2 independently of 
mo, fi and tan/3. 

The low values of the parameter tan/? are already constrained, for a given set of 
parameters, by the Higgs boson searches performed at LEP2. The combination of results 
obtained at the Tevatron and at LEP2, allow to set more robust and stringent constraints 



for the low tan/3 regime, as already tested in Ref. ||118|| and illustrated in Fig. Finally 
constraints form the Higgs sector and the other SUSY particle sector, could be combined 
to derive improved bounds on the cMSSM parameters; a first attempt has been done in 



Ref. | |147| j when combining Higgs and SUSY limits obtained at LEP2 for a cm. energy 



up to 172 GeV. It was shown that the interplay between Higgs and stop mass lower limits 
allows mi/2 to be constrained for low values of tan/3 and mo: for small (large) mixing, the 



Higgs (stop) mass limit provides the most effective constraint. On the bottom of Fig. |34], 
one can see clearly how the lower limit on the LPS is strengthened by taking the Higgs 
information into account for small tan /3. 
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Particle 


Assumptions 


Lower limit 
(GeV/c 2 ) 


Experimental 
sources 


h 


tan/3 >0.8 


78.8 


LEP2 Comb. 


A 


tan/5 >0.8 


79.1 


LEP2 Comb. 


H ± 


Br(H -> cs) + Br(H -> cs) = l 
any tanp 

tan/5 <1.1 or tan/3 >40 


68 
100 


LEP2 Comb. 
CDF/DO 


en 


Br(eR -> exi)=l, AM> 20 GeV/c 2 
Br(g R -> ex?)=l, tan/3 >2 


85 
65 


LEP2 Comb. 
LEP2 




BR(// R — > ) = 1, A M > 20 GeV/c 2 


71 


LEP2 Comb. 


tr 


BR(f R -> rxS)=l, A M > 20 GeV/c 2 


75 


LEP2 Comb. 


V 




43 


LEP1 - Z width 


tl 


BR(ti -> cx?)=l, A M > 15 GeV/c 2 
BR(ti -> i/Zxi ) =1 , A M > 15 GeV/c 2 
BR(ti -> cxi )=1, LEP2; m x o < 50 GeV 


83 
85 
122 


LEP2 Comb. 
LEP2 Comb. 
CDF 


Di 


DD/k . L...0 \ 1 A 1\/T ~~> 1 ^ PoAA / n 2 

-Dit^Di — > 0^ J — 1, ZA 1V1 10 VjreV/C 


i 


T T7PO Pnmli 


q 


m q > nig 
mq > m g 


216 
260 


CDF 
DO 


g 


/i < -100 or /i > 200 GeV/c 2 
/i < -100 or /i > 200 GeV/c 2 


173 
187 


CDF 
DO 


xt 


Higgsino AM > 3 GeV/c 2 
Gaugino m > 200 GeV/c 2 
Gaugino any m 


63 
94.3 
57.2 


LEP2 Comb. 

LEP2* 

LEP2 


xS 


Large m 
Any m 


32.5 
28 


LEP2* 
LEP2* 



Table 16: Summary of lower mass limits obtained at the Tevatron and at LEP2 up to \fs 
= 183, 189 GeV [* at which just a small fraction of data has been used]. 
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Figure 34: Top: DO ||139|| exclusion contours in the plane (mo, ^1/2) from the multi-jet 
plus E/_ or dileptons plus jets plus E/. searches. Bottom: Lower limit on the mass of the 
lightest neutralino as a funtion of tan/3 for a series of rrtQ values by the ALEPH experiment 
147] for a cm. energy up to 172 GeV. The dashed curve shows the result coming from 



the Higgs constraints for mo=75 GeV/c 
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Figure 35: Regions in the (/i,M2) plane excluded by DELPHI [ 143 1 for different values 
of tan/?, assuming large mo- The lightly shaded areas are those excluded by lower energy 
LEPl results. The intermediate shading shows regions excluded by the chargino search at 
183 GeV. The darkly shaded areas show the regions excluded by the neutralino searches 
at these energies. 
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